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INTRODUCTION 


The  purpose  of  our  project  was  to  characterize  national  patterns  in  the  treatment  of 
early  invasive  breast  cancer  in  older  women  with  incident  disease.  We  specifically  sought  to 
characterize  disparities  in  care  and  regional  variation  in  treatment  patterns.  Our  study  subject 
was  in  response  to  prior  research  evidence  that  has  suggested  that  regional  variation  and 
socioeconomic  barriers  in  breast  cancer  treatment  remain  substantial  problems  for  patients 
across  the  nation.  In  fact,  though  a  full  decade  ago,  in  1999,  the  Institute  of  Medicine  National 
Cancer  Policy  Board  issued  a  call  to  improve  the  quality  of  cancer  care  nationally1,  still  recent 
data  indicate  that  progress  in  overcoming  disparities  in  cancer  care  has  been  insufficient2.  For 
our  project,  we  applied  a  novel  resource,  comprehensive  national  Medicare  claims  data,  to 
study  disparities  in  care  and  outcomes  in  women  with  breast  cancer.  Our  project  has  spanned  a 
total  of  three  years. 

At  the  culmination  of  Year  1  of  this  project,  our  main  objective  was  to  characterize  the 
scope  of  treatment  disparities  and  the  magnitude  of  regional  variation  in  care,  using  cross- 
sectional  data.  At  the  culmination  of  Year  2,  we  assessed  the  factors  that  impacted  the  choices 
for  treatment— including  non-standard  and  developing  treatment— as  well  as  the  implications 
of  these  treatment  choices  on  variations  on  spending  (cost)  in  breast  cancer  care.  At  the 
culmination  of  Year  3,  we  assessed  the  impact  of  utilization  patterns  on  breast  cancer 
outcomes.  The  following  narrative  will  summarize  results  obtained  over  Years  1  through  3. 

BODY 

Task  Summary  from  Statement  of  Work  (SOW) 

Task  1.  To  assess  standardized  utilization  rates  of  radiation  therapy  and  chemotherapy. 
Objectives 

Analysis  1:  To  present  overall  national  and  state-by-state  absolute  and  standardized  utilization 
rates  of  treatment 

Analysis  2:  To  present  utilization  rates  of  radiotherapy  (RT)  after  conservative  surgery  (CS)  by 
race,  a  significant  modifying  factor,  in  order  to  quantify  disparities  in  breast  cancer  treatment, 
given  that  RT  after  CS  is  generally  considered  standard  therapy 

Analysis  3:  To  present  utilization  rates  of  brachytherapy,  in  order  to  help  quantify  the  uptake  in 
an  emerging  area  of  treatment  across  the  United  States 

Analysis  4:  To  present  a  validation  sample  of  breast  cancer  patients  and  evaluate  and  validate 
the  utility  of  claims-based  covariates  in  predicting  breast  cancer  stage 
*Please  also  see  Annual  Report,  September  2008. 

Task  2.  To  assess  initial  and  continuing  care  costs  of  breast  cancer  care. 

Objectives 

Analysis  1:  To  present  overall  national  and  state-by-state  absolute  and  standardized  utilization 
rates  of  mastectomy  versus  breast-conserving  surgery  (BCS) 

Analysis  2:  To  identify  variations  and  disparities  in  use  of  mastectomy  versus  BCS 
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Analysis  3:  To  quantify  the  costs  associated  with  mastectomy  and  BCS,  and  to  compare  these 
costs  with  other  breast  cancer  treatment  related  spending,  including  chemotherapy, 
radiotherapy,  and  other  surgeries 

Analysis  4:  To  compare  breast  cancer  associated  costs  with  non  cancer  associated  costs  in  the 
year  after  cancer  diagnosis 

Analysis  5:  To  identify  predictors  of  breast  cancer-related  costs,  including  treatment,  disease, 
patient,  and  socioeconomic  factors 

Analysis  6:  To  present  utilization  rates  of  brachytherapy,  an  emerging  but  potentially  costly 
breast  cancer  treatment,  in  order  to  help  quantify  the  diffusion  of  this  technology  in  an 
emerging  area  of  treatment  across  the  United  States  and  discuss  the  implications  of  findings  on 
the  future  incorporation  of  costly  treatments  and  technologies. 

*Please  also  see  Annual  Report,  September  2009. 

Task  3.  To  assess  outcomes  relevant  to  breast  cancer. 

Objectives 

Analysis  1:  To  create  a  longitudinal  cohort  of  breast  cancer  patients  treated  with  BCS,  with  or 
without  RT. 

Analysis  2:  To  examine  overall  survival  rates  in  patients  treated  with  vs.  without  RT. 

Analysis  3:  To  examine  geographic  variation  in  survival. 

Methods 

Study  sample 

Our  main  study  sample  of  patients  diagnosed  with  invasive  breast  cancer  and  treated 
with  either  BCS  or  mastectomy  was  derived  from  the  national  Medicare  dataset.  In  addition,  we 
used  several  subsamples  in  our  analysis  to  accomplish  the  above  mentioned  objectives,  which 
will  also  be  detailed  below. 

Medicare.  The  national  Medicare  dataset  includes  comprehensive  claims  information 
with  beneficiary-specific  data  on  all  Medicare  beneficiaries  in  the  United  States.  Files  contain 
data  collected  by  Medicare  for  reimbursement  of  health  care  services  for  each  beneficiary  and 
include  institutional  (inpatient  and  outpatient)  as  well  as  non-institutional  (physician  or  other 
providers)  final  action  claims3. 

Cross-sectional  cohort.  To  define  a  cohort  of  patients  with  incident  disease  in  2003 
required  claims  data  spanning  2002  to  2004  to  have  complete  information  on  the  claims  history 
the  year  prior  to  diagnosis  and  information  on  claims  up  a  year  after  diagnosis,  as  detailed 
below.  In  summary,  our  initial  study  population  consisted  of  853,273  women  who  had  any 
diagnosis  of  invasive  breast  cancer  in  2003,  based  on  an  International  Classification  of  Diseases, 
Ninth  Revision  (ICD-9)  diagnosis  code  of  174.  As  this  denominator  would  have  included  incident 
and  prevalent  cases  in  2003,  we  used  the  following  algorithm  to  identify  patients  with  incident 
breast  cancer,  treated  with  either  breast-conserving  surgery  (BCS)  or  mastectomy.  This  method 
was  based  on  a  prior  validated  algorithm  for  claims  data.4, 5  We  included  women  (age>65  years) 
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who  had  any  diagnosis  of  invasive  breast  cancer  in  2003  (defined  as  an  International 
Classification  of  Diseases,  Ninth  Revision  (ICD-9)  diagnosis  code  of  174)  who  underwent  BCS  (N= 
83,611)  or  mastectomy  (N=42,504)  between  January  1,  2003,  and  December  31,  2003.  From 
this  sample,  to  increase  the  specificity  of  the  definition,  we  excluded  23,715  patients  who  did 
not  have  at  least  2  claims  (on  different  dates)  indicating  a  diagnosis  of  invasive  breast  cancer 
between  January  1,  2003  and  December  31,  2004  (at  least  1  claim  must  have  occurred  during 
2003),  with  no  more  than  6  months  between  the  date  of  BCS  or  mastectomy  and  the  earliest 
breast  cancer  diagnosis  claim  date.  To  exclude  the  prevalent  cases,  we  excluded  16,471 
patients  who  had  a  breast  cancer-related  diagnosis  or  procedure  claim  between  January  1, 
2002,  and  December  31,  2002.  To  reduce  misclassification  of  the  primary  intended  surgery,  we 
excluded  630  patients  who  underwent  both  types  of  surgery  (either  date  of  mastectomy  claim 
preceding  date  of  BCS  claim  or  mastectomy  occurring  more  than  3  months  after  BCS),  except 
for  patients  who  received  a  mastectomy  within  3  months  of  BCS  for  whom  mastectomy  was 
considered  the  definitive  surgery.  To  limit  our  sample  to  patients  with  non-metastatic  invasive 
breast  cancer,  we  then  excluded  2,122  patients  who  had  two  or  more  claims  specifying 
metastatic  breast  cancer  from  3  months  before  to  3  months  after  the  diagnosis  date.  To 
improve  sample  homogeneity,  we  also  excluded  5,719  patients  who  were  receiving  Medicare 
coverage  because  of  end-stage  renal  disease  or  disability.  Finally,  to  ensure  we  had  complete 
claims  information  to  determine  patients'  cancer  treatment  course  and  comorbidities,  we 
excluded  6,612  patients  who  lacked  Part  A  or  B  coverage  or  who  had  intermittent  health 
maintenance  organization  coverage  in  the  9  months  after  or  in  the  1  year  before  their  breast 
cancer  diagnosis  date  (of  these  patients,  3,561  had  incomplete  information  in  the  year  prior  to 
diagnosis  because  they  were  <66  years  of  age).  For  this  analysis,  the  breast  cancer  diagnosis 
date  was  defined  as  the  date  of  the  earliest  claim  for  a  diagnosis  of  breast  cancer.  This  left  a 
final  sample  size  of  56,725  patients.  Our  algorithm  was  based  on  a  prior  validated  algorithm  for 
identifying  breast  cancer  patients  using  claims  data4. 

Longitudinal  cohort.  We  subsequently  derived  a  longitudinal  cohort,  which  included 
breast  cancer  patients  (as  defined  by  the  algorithm  above)  diagnosed  between  2000  and  2003. 
In  addition,  we  retained  all  denominator  and  claims  data  for  these  patients  through  2004,  thus 
allowing  for  up  to  four-year  follow-up  in  these  patients. 

Surgical  treatment  and  other  covariates 

Covariates  derived  from  Medicare  files  (denominator  and  claims  files)  included  cancer 
treatment  variables,  other  clinical  variables,  and  demographic  data.  Patients  were  classified  as 
treated  with  BCS  or  mastectomy  if  a  claim  for  the  surgery  (Appendix  A)  occurred  within  6 
months  of  the  breast  cancer  diagnosis  date.  Claims  for  chemotherapy  must  have  occurred 
within  6  months  and  RT  claims  within  9  months  of  the  breast  cancer  diagnosis  date.  These 
claims-based  treatment  definitions  have  been  applied  in  prior  studies  of  breast  cancer 
patients.6  11 

Other  disease-  and  treatment-related  variables  included  axillary  lymph  node 
involvement,  axillary  lymph  node  dissection,  sentinel  node  biopsy,  receipt  of  any  RT  (including 
brachytherapy),  receipt  of  any  chemotherapy,  specific  receipt  of  doxorubicin  or  paclitaxel. 
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receipt  of  any  imaging  studies  for  staging,  number  of  hospitalizations  in  the  year  after 
diagnosis,  and  number  of  medical  oncology,  radiation  oncology,  and  surgery  visits  in  the  year 
after  diagnosis.  Of  patients  who  received  RT,  patients  were  further  classified  as  having  received 
external  beam  radiation  therapy  (EBRT),  brachytherapy,  or  both  (EBRT  plus  brachytherapy 
boost),  as  indicated  by  claims  codes.  Patients  treated  with  brachytherapy  were  further 
classified  as  having  received  balloon-based  treatment  if  a  procedure  code  was  found  specifying 
accordingly. 

Variables  indicating  preventive  healthcare  and  interactions  with  the  healthcare  system 
included  mammography  in  the  year  prior  to  diagnosis  and  number  of  physician  visits  in  the  year 
prior  to  diagnosis.  In  addition,  we  calculated  the  severity  of  comorbid  disease  for  each  patient 
based  on  a  modified  Charlson  comorbidity  score  validated  in  a  prior  claims-based  study:  0  (no 
comorbidity),  1  (mild  to  moderate),  or  2  or  more  (severe)12.  This  score  combined  comorbidities 
recorded  in  Medicare  claims  during  the  12  months  prior  to  the  patient's  cancer  diagnosis.  To 
enhance  specificity  of  comorbid  disease  diagnoses,  patients  must  have  had  at  least  1  inpatient 
(Part  A)  claim  or  at  least  2  outpatient  (Part  B)  claims  more  than  30  days  apart.12 

Demographic  data  available  through  Medicare  files  included  patient  age  at  diagnosis, 
race  (categorized  as  white,  black,  and  other),  and  state  and  county  of  residence.  Classification 
of  geographic  regions  was  based  on  Census  Divisions  definitions.13  Socioeconomic  variables, 
obtained  from  the  2003  Area  Resource  File  (ARF)14  and  linked  to  the  Medicare  dataset  by 
patients'  county,  included  (by  county  of  patient's  residence)  median  household  income, 
percent  living  below  poverty  level,  percent  completing  ninth  grade  education,  high  school,  and 
college.  Supply  of  healthcare  providers  (for  breast  cancer  treatment)  was  measured  by  the 
density  of  general  surgeons,  and  density  of  radiation  oncologists  at  county-level,  obtained  from 
the  2003  ARF. 

Predictor  covariates  were  obtained  by  searching  through  inpatient,  outpatient,  and 
carrier  Medicare  claims  or  the  denominator  file  for  SEER-Medicare  linked  data  for  demographic 
variables.  A  comprehensive  list  of  International  Classification  of  Diseases,  Ninth  Revision  (ICD- 
9),  Common  Procedural  Terminology  (CPT),  and  Revenue  Center  codes  for  each  predictor  are 
listed  in  Table  1. 

Cost 

Total  health  care  costs  for  each  patient  were  calculated  based  on  Medicare  spending. 
Any  claim,  and  the  associated  total  payment  amount  reported  by  Medicare,  identified  through 
the  inpatient,  outpatient,  or  carrier  claims  files  was  added,  for  a  sum  total  of  costs  in  the  year 
after  breast  cancer  diagnosis.  Breast-cancer  related  spending  was  summed  for  each  patient 
with  claims  belonging  in  the  following  categories:  surgery  (mastectomy,  BCS,  or  reconstruction 
after  mastectomy),  axillary  nodal  dissection  or  involvement,  chemotherapy,  radiotherapy,  and 
other  breast  cancer  diagnosis-related  claim. 

Survival 

All-cause  mortality  is  tracked  through  Medicare  claims,  specifically  through  the 
denominator  data.  Thus  we  searched  longitudinal  claims  to  obtain  a  date  of  death.  If  date  of 
death  was  missing,  then  the  patient  was  assumed  to  be  alive  at  last  follow-up. 


Statistical  analyses 

Statistical  analyses  were  conducted  using  SAS  version  9.1.3  (SAS  Institute  Inc,  Cary,  NC), 
and  all  statistical  tests  assumed  a  2-tailed  a  of  0.05. 

Treatment  utilization.  (See  details  in  Annual  Report,  September  2008  and  Annual 
Report,  September  2009)  For  characterization  of  treatment  utilization,  bivariate  associations 
between  treatment  (radiotherapy  and  chemotherapy;  percent  BCS  versus  mastectomy)  and 
other  covariates  were  tested  using  the  Pearson  chi-square  test  for  categorical  variables  and  the 
Wilcoxon  rank  sum  test  for  continuous  variables.  Percent  BCS  versus  mastectomy  use  was 
calculated  for  the  entire  sample,  by  state,  and  by  region.  A  multivariate  logistic  model  tested 
the  adjusted  association  between  predictors  and  treatment.  Covariates  were  selected  a  priori 
based  on  significance  in  bivariate  analyses  (P<0.25)  and  significance  in  prior  studies  of  cancer 
patients*6"121.  Logistic  models  also  derived  standardized  treatment  rates,  based  on  unadjusted 
percent  use  of  radiotherapy  and  chemotherapy,  standardized  for  covariates. 

Outcomes.  (See  details  in  Annual  Report,  September  2009)  Mean  breast  cancer-related 
and  overall  healthcare  costs  were  calculated  for  the  entire  sample,  by  state,  and  by  region.  A 
multivariate  logistic  model  tested  adjusted  associations  with  BCS  and  a  multivariate  linear 
regression  model  tested  adjusted  associations  with  breast  cancer-related  costs,  with  covariates 
for  the  final  multivariate  models  selected  a  priori  based  on  significance  in  bivariate  analyses 
(P<0.25)  and  significance  in  prior  studies  of  cancer  patients19"24  The  final  parsimonious  model 
was  selected  based  on  statistical  significance,  goodness  of  fit,  and  minimizing  multicollinearity. 
Finally,  percent  died  over  longitudinal  follow-up  was  calculated  for  the  entire  sample,  by  state, 
and  by  region. 

Early  stage  breast  cancer  subgroup 

Breast  cancer  stage  is  not  directly  available  through  Medicare  claims  data.  To  select  a 
subgroup  of  patients  with  early  stage  breast  cancer,  given  that  surgical  treatment  strategies  are 
typically  dependent  on  disease  stage,  we  applied  a  previously  validated  predictive  algorithm 
that  used  claims-based  covariates  to  identify  patients  with  a  high  probability  of  having  stage  I  or 
II  disease.15  (See  also.  Annual  Report,  September  2008).  Therefore,  in  the  selected  subgroup  of 
43,706  predicted  early  stage  patients,  we  further  examined  the  adjusted  associations  between 
receipt  of  BCS  and  covariates  using  multivariate  logistic  modeling. 

In  a  secondary,  validating  analysis  on  this  selected  group  of  predicted  early  stage 
patients,  we  also  identified  the  subgroup  of  42,499  patients  who  did  not  have  claims  for  axillary 
involvement  and  chemotherapy,  as  these  patients  would  also  be  more  likely  to  have  early  stage 
disease  (kappa=0.73  for  the  two  selected  subgroups).  The  adjusted  associations  were  compared 
in  this  group  to  the  associations  calculated  for  the  predicted  early  stage  patients. 

Brachytherapy  pilot  subgroup 

(See  also.  Annual  Report,  September  2009)  A  limitation  of  the  current  national 
Medicare  data-based  cohort  for  our  analysis  is  that,  for  the  study  of  treatment  patterns,  this  is 
essentially  a  cross-sectional  sample  (as  the  initial  course  of  cancer  treatment  can  occur  months 
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or  even  up  to  a  year  after  the  date  of  diagnosis).  Thus,  for  the  study  of  temporal  patterns,  the 
national  Medicare  data  is  particularly  limited.  Therefore,  we  sought  to  study  a  pilot  sample 
derived  from  the  The  MarketScan  Medicare  Supplemental  database  is  a  large,  nationwide, 
employment-based  claims  database  which  includes  Medicare  beneficiaries  with  private 
supplemental  insurance  obtained  through  their  former  employers.  We  identified  6,882  women 
aged  65  years  and  older  with  a  diagnosis  of  invasive  breast  cancer  and  treated  with  BCS.  Using 
the  same  algorithm  as  listed  above,  patients  who  had  claims  for  RT  were  further  classified  as 
having  received  as  treated  with  external  beam  radiation  therapy  (EBRT),  brachytherapy,  or  both 
(EBRT  plus  brachytherapy  boost),  as  indicated  by  claims  codes.  Patients  treated  with 
brachytherapy  were  further  classified  as  having  received  balloon-based  treatment  if  a 
procedure  code  was  found  specifying  accordingly.  To  address  our  second  objective,  we 
evaluated  for  a  time  trend  using  the  Mantel-Hanszel  chi-square  and  Cochran-Armitage  tests  for 
trend.  We  also  benchmarked  the  time  trends  against  two  major  policy  events  that  occurred 
during  the  study  period:  Food  and  Drug  Administration  (FDA)  approval  of  the  first  balloon- 
based  brachytherapy  device  for  breast  cancer  in  June  2002;  and  Medicare  reimbursement  of 
breast  brachytherapy  in  April  2004. 

Table  1.  Claims  codes  searched  to  calculate  costs  and  define  variables  of  interest 


Variables 


ICD-9  Diagnosis  ICD-9  Procedure  CPT  Revenue  Center 


Radiotherapy 


EBRT 


9221-6,  9228,  9229 


77427,  77431-2,  77401 
9,  77411-14,  77416, 
77418,  77470,  77499, 
77520,  77522-3, 
77525,  77750,  77789, 
77790 


0330,  0333 


Brachytherapy 


19296-8,  77326-8, 
77761-3,  77776-8, 
77781-4,  Q3001 , 
C9714-5 


Balloon-based  brachytherapy 


C9714-5 


Extent  of  disease  at  diagnosis 


Axillary  LN  involvement 


1963 


Metastatic  disease 


1962,  1965-6,  197, 
1970,  1971,  1972, 
1973,  1974,  1975, 
1976,  1977,  1978,  198, 
1980,  1981,  1982, 
1983,  1984,  1985, 
1986,  1987,  1988, 
19881,  19882,  19889 
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Table  1,  continued 


Variables  ICD-9  Diagnosis 

Cancer  diagnosis  and  treatment 

Imaging  (CT,  MRI,  PET,  or  bone  scan) 

Breast  conserving  surgery 
Mastectomy 

Axillary  surgery  (LN  dissection  or  sentinel  node) 

Chemotherapy  (any  agent)  V581 ,  V662,  V672 

Preventive  care  and  interaction  with  healthcare  system 

No.  physician  visits 

Screening  mammography  V7611,  V7612 

Influenza  vaccine  V0481 

General  health  status 

Charlson  comorbidity  score3 

Abbreviations:  CPT  Common  Procedural  Terminology;  ICD  International  Cl, 
3  Klabunde  CN,  et  al.  J  Clin  Epidemiol  2000;53:1258-1267. 


ICD-9  Procedure 


CPT 


Revenue  Center 


70450,  70460,  70470, 

70551-3,  71250, 

8801,  8703,  8741, 

71260,  71270,  72192- 

8891,  8896,  8874, 

4,  74150,  74160, 

9218,  9214 

74170,  76700,  78315, 

78320,  78812-6, 

G0213-5 

8520-3,8525 

19110,  19120,  19125, 

19160,  19162 

8541-8 

19180,  19182,  19200, 

19220,  19240 

38500,  38525,  38740, 

4023,  4051 ,  8543, 

38745,  38792,  19162, 

8547 

19200,  19220,  19240, 

78195 

9925 

96400-96549,  J9000-9, 

Q0083-5  0331,0332,0335 

77055-6,  77058-9, 

8737,  8736  76090-2,  G0202,  0401,  0403 

G0204,  G0206 


90658,  G0008 


li 


Results 


Patient  characteristics  and  treatment  course 

Our  cohort  consisted  of  56,725  women  with  incident,  invasive  breast  cancer  diagnosed 
in  2003  and  treated  with  surgery.  In  our  sample,  median  age  was  76  (interquartile  range  71  to 
81).  Ninety  percent  (N=51,432)  were  white,  7%  (N=3,727)  were  black,  and  3%  (N=l,566)  were 
of  other  race.  As  a  component  of  the  initial  treatment  course,  the  majority  of  patients  were 
treated  with  BCS.  Specifically,  59%  of  patients  (N=33,450)  received  BCS,  while  41%  (N=23,275) 
underwent  mastectomy.  Additionally,  of  the  entire  sample,  50%  received  RT  and  16%  received 
chemotherapy.  Of  BCS  patients,  74%  received  RT  and  13%  received  chemotherapy.  Of 
mastectomy  patients,  14%  received  RT  and  23%  chemotherapy. 


Part  I:  (for  further  details,  please  see  Annual  Report,  September  2008). 

Predictors  of  chemotherapy  and  RT 

Utilization  of  RT  and  chemotherapy  across  the  United  States 

Absolute  percent  use  and  standardized  rates 

Seventy-three  percent  of  women  received  RT  after  CS  and  13%  received  chemotherapy 
as  part  of  their  initial  treatment  course.  However,  there  was  significant  variation  in  the 
utilization  of  these  treatments  across  the  United  States,  with  as  little  as  50%  to  as  much  as  85% 
utilization  of  RT  after  CS  by  state  (P<0.001);  while  utilization  of  chemotherapy  ranged  from  8% 
to  22%  by  state  (P<0.001)  (Figure  la,  2a).  After  standardization  of  utilization  rates  by 
covariates,  significant  geographic  variation  still  existed  (Figure  lb,  2b). 

Figure  la.  RT  in  the  United  States  (all  women).  Darker  shading  represents  higher  percent 
utilization.  Categories:  50%  to  60%,  61%  to  65%,  66%  to  70%,  71%  to  73%,  74%  to  75%,  76%  to 
85%,  86%  to  100% 


m  "p 

■f 
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Figure  lb.  RT  in  the  United  States  (standardized  risk  adjusted  rates) 


Figure  2a.  Chemotherapy  in  the  United  States  (all  women).  Categories:  8%  to  11%,  12%  to 
15%,  16%  to  18%,  19%  to  22% 


Figure  3.  RT  use  in  white  women.  Categories:  12%  to  16%,  17%  to  33%,  34%  to  50%,  51%  to 
66%,  67%  to  72%,  73%  to  75%,  76%  to  77%.  Gray=  data  insufficient. 


Figure  4.  RT  use  in  black  women. 
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Utilization  of  RT  and  by  race 


Significant  disparities  existed  in  RT  treatment  by  race,  with  74%  of  white  women,  65%  of 
black  women,  and  66%  of  other  women  who  received  RT  after  CS  (P<0.001)  (Figures  3-4).  After 
adjusting  for  covariates,  white  women  were  still  significantly  more  likely  than  black  women  to 
have  received  RT  after  CS  (Odds  Ratio  [OR]  1.48,  95%  Confidence  Interval  [Cl]  1.34-1.63, 
P<0.001).  There  was  a  trend  toward  white  women  also  being  more  likely  than  other  (non-white, 
non-black)  women  to  have  received  RT  after  CS,  but  this  difference  was  marginally  significant 
(OR=  1.22,  95%  Cl  1.04-1.42,  P=0.01)  (Table  2). 

Table  2.  Unadjusted  and  adjusted  utilization  of  radiotherapy  (RT)  by  race 

Model:  RT  use  OR  95%  Cl  P 


Unadjusted 


White  vs.  black  patients 

1.52 

1.38 

1.67 

<0.001 

White  vs.  other+  patients 

1.46 

1.27 

1.68 

<0.001 

Adjusted* 

White  vs.  black  patients 

1.48 

1.34 

1.63 

<0.001 

White  vs.  other+  patients 

1.22 

1.04 

1.42 

0.01 

*  Adjusted  for  age,  comorbidity,  chemotherapy  (adriamycin  or  taxol),  axillary  lymph  node  involvement,  staging  imaging,  surgeon  visits, 
mammography,  physician  visits,  and  region. 


Race  and  utilization  in  younger  patients 

In  the  subset  of  patients  age  70  and  younger,  in  whom  RT  utilization  after  CS  would  be 
expected  to  be  most  common,  racial  disparities  persisted,  particularly  between  white  women 
versus  black  women.  A  total  of  83%  of  white  women,  72%  of  black  women,  and  78%  of  other 
women  in  this  younger  group  received  RT  after  CS  (P<0.001).  The  disparity  between  white  and 
back  women  persisted  even  after  adjusting  for  covariates  (OR=  1.81,  95%  Cl  1.46-2.25, 
P<0.001). 

Race  and  geographic  variation  in  utilization 

There  was,  however,  substantial  geographic  variation  in  racial  disparities.  Regions  with 
the  most  marked  racial  disparities  included  the  Pacific  West,  the  East  South  Central  region,  and 
the  Northeast  (Table  3).  Black  women  fared  particularly  poorly  in  these  regions,  with  less  than 
60%  RT  utilization  after  CS.  (Note  that  regional  and  state  data  were  not  presented  for  "other" 
race  due  to  insufficient  numbers  of  patients  in  this  category). 


14 


Table  3.  Racial  disparities  in  RT  utilization  after  CS  by  region 


Region 

States 

%  in  Whites 

%  in  Blacks  P 

Pacific  West 

CA,  OR, 

WA 

72 

55 

<0.001 

Mountain  West 

AZ,  CO, 

ID,  MT,  NV,  NM,  UT,  WY 

76 

74 

0.81 

Midwest,  West  North  Central 

IA,  KS,  MN,  MO,  NE,  ND,  SD 

74 

72 

0.77 

Midwest,  East  North  Central 

IL,  IN,  Ml,  OH,  Wl 

76 

71 

0.04 

Northeast,  New  England 

CT,  ME, 

MA,  NH,  NJ,  NY,  VT 

70 

58 

<0.001 

Northeast,  Mid-Atlantic 

DE,  DC, 

MD,  PA,  Rl 

72 

59 

<0.001 

Southwest,  West  South  Central 

AR,  LA, 

OK,  TX 

73 

64 

0.003 

Southeast,  East  South  Central 

AL,  KY, 

MS,  TN 

72 

57 

<0.001 

Southeast,  South  Atlantic 

FL,  GA, 

NC,  SC,  VA,  WV 

77 

69 

<0.001 

Part  II  (for  further  details,  please  see  Annual  Report,  September  2009). 

Predictors  of  BCS  versus  mastectomy  use 

The  use  of  BCS  was  associated  with  both  clinical  and  non-clinical  factors.  Patients  who 
were  younger,  white,  and  had  fewer  comorbidities,  lymph  node-negative  disease,  and 
predicted  early  stage  disease  were  more  likely  to  undergo  BCS.  Patients  who  did  not  receive 
chemotherapy  or  did  not  undergo  axillary  surgery  were  also  more  likely  to  undergo  BCS  (Table 
4).  In  addition,  neighborhood  socioeconomic  factors  were  also  highly  associated  with  BCS. 
Specifically,  patients  living  in  metropolitan  areas  and  in  counties  with  higher  median  household 
income,  lower  percent  living  below  poverty  level,  and  higher  percent  with  college  education 
were  more  likely  to  undergo  BCS.  Of  patients  living  in  non-metropolitan  areas,  only  51% 
underwent  BCS.  Supply  of  healthcare  providers  also  influenced  treatment,  with  BCS  use  more 
likely  in  patients  residing  in  counties  with  a  higher  density  of  surgeons  and  radiation  oncologists 
(Table  4).  Finally,  significant  geographic  variation  existed  (P<0.001),  with  patients  in  the 
Northeast  and  Pacific  West  most  likely  to  undergo  BCS.  In  contrast,  patients  in  the  South  were 
least  likely  to  undergo  BCS,  with  half  or  fewer  of  all  patients  in  these  regions  treated  with  BCS 
(Figure  5,  Table  5).  On  adjusted  analysis,  higher  density  of  surgeons  was  no  longer  a  significant 
predictor  of  BCS  use  (P=0.13),  specifically  once  the  multivariate  model  accounted  for 
geographic  region.  However,  a  higher  density  of  radiation  oncologists  remained  a  significant 
predictor  of  BCS  use  (P=0.01).  (Table  6). 

Early  stage  breast  cancer  subgroup 

In  the  selected  subgroup  of  43,706  patients  with  predicted  early  stage  (Stage  I  or  II) 
disease  (77%  of  the  entire  sample),  a  total  of  68%  (29,828  of  43,706)  of  this  selected  group 
received  BCS.  This  was  consistent  with  a  total  of  65%  (27,544  of  42,499)  received  BCS  in  the 
validation  subgroup  of  patients  who  did  not  have  axillary  involvement  and  did  not  receive 
chemotherapy.  Geographic  variation  persisted  in  the  use  of  BCS  for  patients  with  predicted 
early  stage.  Patients  in  the  Northeast  (78-79%)  and  Pacific  West  (71%)  were  still  the  most  likely 
to  undergo  BCS,  while  patients  in  the  South  (57-59%)  and  portions  of  the  Midwest  (58%)  were 
the  least  likely  (Table  5).  The  validation  subgroup  was  similar,  with  BCS  ranging  from  54% 
(South)  to  75%  (Northeast).  On  adjusted  analysis,  significant  predictors  of  BCS  use  included 
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similar  demographic,  clinical,  and  socioeconomic  factors  as  predictors  for  the  entire  sample 
(Table  6).  Significant  predictors  of  BCS  use  identified  in  the  validation  group  for  early  stage 
breast  cancer  were  also  consistent  with  this  analysis. 

Table  4.  Predictors  of  treatment  utilization 


Patient  Characteristic 

%  Treated  with  BCS  (N) 

%  Treated  with  Mastectomy  (N) 

P 

Demographic 

Age,  mean  (standard  deviation) 

76  (7) 

77  (7) 

<0.001 

66  to  <70  years 

63  (7,112) 

37  (4,260) 

<0.001 

>70  years 

58  (26,338) 

42  (19,015) 

Race 

White 

59  (30,572) 

41  (20,860) 

<0.001 

Black 

54  (2,026) 

46  (1,701) 

Other 

54  (852) 

46  (714) 

Clinical 

Charlson  comorbidity  score 

0  comorbid  conditions 

61  (22,735) 

39  (14,578) 

<0.001 

1  comorbid  condition 

57  (6,414) 

43  (4,875) 

2  or  more  comorbid  conditions 

54  (2,673) 

46  (2,314) 

Unknown 

52  (1,628) 

48  (1,508) 

Disease  Stage  and  Treatment 

Predicted  early  stage  (stage  1  or  II)  disease 

68  (29,828) 

32  (3,622) 

<0.001 

Axillary  lymph  node  positive  disease 

37  (17,587) 

63  (30,050) 

<0.001 

Lymph  node-negative  disease 

63  (5,688) 

37  (3,400) 

Axillary  lymph  node  dissection 

42  (11,995) 

58  (16,715) 

<0.001 

No  axillary  dissection 

77  (21,455) 

23  (6,560) 

Radiation  therapy 

88  (24,823) 

12  (3,358) 

<0.001 

No  radiation  therapy 

30  (8,627) 

70  (19,917) 

Chemotherapy 

44  (4,299) 

56  (5,447) 

<0.001 

No  chemotherapy 

62  (29,151) 

38  (17,828) 

Healthcare  access* 

Median  surgeon  density*  (IQR) 

11 (8-16) 

10  (6-15) 

<0.001 

Median  radiation  oncologist  density*  (IQR) 

13  (4-20) 

11  (0-19) 

<0.001 

Socioeconomic  status* 

Living  in  metropolitan  area 

62  (25,979) 

38  (16,132) 

<0.001 

Living  in  non-metropolitan  area 

51  (7,193) 

49  (6,971) 

Median  income  (IQR) 

41,691  (36,221-48,059) 

39,879  (34,267-45,922) 

<0.001 

Median  percent  living  in  poverty  (IQR) 

10.7  (8-14) 

11.4  (9-14) 

<0.001 

Median  percent  with  college  education  (IQR) 

24  (17-29) 

22  (15-28) 

<0.001 

Abbreviations:  BCS  breast-conserving  surgery;  IQR  interquartile  range 

*  By  patient  county  of  residence 
+  Per  100,000  persons 
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Figure  5.  Percent  use  of  BCS  by  state  (darker  shading  represents  higher  frequency  of  use). 
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Table  5.  Unadjusted  and  adjusted  percent  use  of  breast-conserving  surgery  (BCS)  versus 


mastectomy  by  geographic  region. 

%  BCS  (N=23,275)  %  Mastectomy  (N=33,450) 


Region 

States 

Overall 

Adjusted* 

Early  Stage 

Overall 

Adjusted* 

Early  Stage 

West,  Pacific  West 

AK,  CA,  HI,  OR,  WA 

62 

62 

71 

38 

38 

29 

West,  Mountain  West 

AZ,  CO,  ID,  MT,  NV,  NM,  UT, 

WY 

57 

58 

66 

43 

42 

34 

Midwest,  West  North  Central 

IA,  KS,  MN,  MO,  NE,  ND,  SD 

50 

51 

58 

50 

49 

41 

Midwest,  East  North  Central 

IL,  IN,  Ml,  OH,  Wl 

61 

60 

70 

39 

40 

30 

Northeast,  New  England 

CT,  MA,  NH,  ME,  Rl,  VT 

70 

63 

79 

30 

36 

21 

Northeast,  Mid-Atlantic 

NJ,  NY,  PA 

67 

64 

78 

33 

35 

22 

South,  South  Atlantic 

DE,  DC,  FL,  GA,  MD,  NC,  SC, 

VA,  WV 

59 

59 

68 

41 

41 

32 

South,  West  South  Central 

AK,  LA,  OK,  TX 

50 

54 

59 

50 

45 

41 

South,  East  South  Central 

AL,  KY,  MS,  TN 

48 

51 

57 

52 

48 

43 

Abbreviations:  BCS  breast-conserving  surgery 

*  Adjusted  for  covariates  including  age,  race,  comorbidity  score,  axillary  lymph  node  involvement,  axillary  dissection,  chemotherapy,  screening 
mammography,  physician  visits,  surgeon  density,  radiation  oncologist  density,  metropolitan  area,  poverty,  education. 
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Table  6.  Multivariate  logistic  model:  Predictors  of  utilization  of  breast-conserving  surgery, 


Entire  sample  Early  stage  subgroup 


Covariate 

OR 

95%  Cl 

P 

OR 

95%  Cl 

P 

Age  66  to  <70  years  vs.  >70  years 

1.37 

1.31 

1.44 

<0.001 

1.15 

1.09 

1.22 

<0.001 

Race 

White  vs.  black  race 

1.14 

1.05 

1.23 

<0.001 

0.95 

0.86 

1.04 

0.26 

White  vs.  other  race 

1.29 

1.15 

1.45 

<0.001 

1.17 

1.08 

1.27 

<0.001 

Charlson  comorbidity  score 

0  vs.  1  comorbid  conditions 

1.18 

1.13 

1.24 

<0.001 

1.21 

1.15 

1.28 

<0.001 

0  vs.  2  or  more  comorbid  conditions 

1.38 

1.29 

1.47 

<0.001 

1.17 

1.08 

1.27 

<0.001 

Ovs.  unknown  comorbid  conditions 

1.12 

1.01 

1.25 

0.03 

0.91 

0.79 

1.05 

0.20 

Lymph  node-negative  disease 

1.60 

1.52 

1.68 

<0.001 

- 

- 

- 

- 

No  axillary  lymph  node  dissection* 

4.00 

3.85 

4.17 

<0.001 

4.65 

4.46 

4.88 

<0.001 

No  chemotherapy* 

1.32 

1.25 

1.39 

<0.001 

0.63 

0.57 

0.69 

<0.001 

Screening  mammography§ 

2.02 

1.87 

2.17 

<0.001 

0.89 

0.80 

1.00 

0.04 

1  or  more  visits  to  physician§ 

1.43 

1.22 

1.68 

<0.001 

1.97 

1.60 

2.43 

<0.001 

Geographic  Region11 

West,  Pacific  West 

0.86 

0.78 

0.95 

0.008 

0.81 

0.72 

0.92 

0.001 

West,  Mountain  West 

0.71 

0.63 

0.80 

<0.001 

0.69 

0.60 

0.80 

<0.001 

Midwest,  West  North  Central 

0.50 

0.45 

0.55 

<0.001 

0.45 

0.40 

0.51 

<0.001 

Midwest,  East  North  Central 

0.80 

0.73 

0.88 

<0.001 

0.75 

0.67 

0.84 

<0.001 

Northeast,  Mid-Atlantic 

1.05 

0.95 

1.16 

0.34 

1.06 

0.94 

1.19 

0.36 

South,  South  Atlantic 

0.76 

0.69 

0.83 

<0.001 

0.69 

0.62 

0.78 

<0.001 

South,  West  South  Central 

0.58 

0.53 

0.65 

<0.001 

0.54 

0.48 

0.61 

<0.001 

South,  East  South  Central 

0.51 

0.45 

0.57 

<0.001 

0.45 

0.39 

0.52 

<0.001 

Surgeon  density^ 

1.00 

0.99 

1.01 

0.13 

1.00 

0.99 

1.01 

0.14 

Radiation  oncologist  density^ 

1.30 

1.06 

1.59 

0.01 

1.37 

1.07 

1.75 

0.01 

Living  in  metropolitan  area 

1.20 

1.14 

1.26 

<0.001 

1.20 

1.13 

1.27 

<0.001 

Percent  living  in  poverty  <  11% 1 

1.05 

1.00 

1.09 

0.03 

1.05 

1.00 

1.11 

0.06 

Percent  with  college  education  >  23%' 

1.13 

1.08 

1.19 

<0.001 

1.15 

1.09 

1.21 

<0.001 

Abbreviations:  Cl  confidence  interval,  OR  odds  ratio 

*  A  model  excluding  axillary  lymph  node  dissection  and  chemotherapy,  which  are  treatments  likely  to  occur  concurrently  or  after  surgery, 
did  not  affect  risk  estimates  for  other  covariates, 
t  Increased  odds  per  1  surgeon  per  100,000  persons 
§  Increased  odds  per  1  radiation  oncologist  per  10,000  persons 

I  Continuous  variables  dichotomized  at  the  median  value 

I I  Compared  with  reference  category  Northeast  New  England.  The  Likelihood  ratio  test  for  all  strata  of  the  variable  for  region  was 
statistically  significant  (PcO.001). 
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Part  III 

Costs  of  cancer  and  non-cancer  care 

In  the  year  following  diagnosis,  the  median  breast  cancer  care-related  costs  (as 
reimbursed  by  Medicare)  was  $6,101  (interquartile  range  [IQR]  $2,900  to  $13,058).  The  most 
costly  contributor  to  these  costs  was  surgery-related  costs  (Table  7).  These  cancer-related  costs 
compare  with  the  median  overall  health  care  costs  of  $12,274  (IQR  $7,623  to  $19,041)  and 
median  non-cancer-related  costs  of  $4,376  (IQR  $668  to  $9042).  There  was  a  significant 
correlation  between  cancer-related  costs  and  overall  health  care  costs  in  our  sample  (Pearson's 
R=  0.56,  P<0.001). 

Table  7.  Median  costs  of  breast  cancer  care 


Category 

Median  costs  ($) 

IQR 

Mean(  SD) 

Total  breast  cancer-related  costs 

Surgery 

6,101 

2,900  - 13,058 

9,973  (11,726) 

Breast  conserving  surgery 

411 

163  -  865 

727  (1,400) 

Mastectomy 

626 

163  -  1,017 

1,098  (1,783) 

Axillary  treatment 

643 

212  -  1,086 

1,097  (2,007) 

Reconstruction 

1,045 

210-2,171 

1,511  (1,729) 

Chemotherapy 

895 

205  -  3,401 

2,524  (4,055) 

Radiotherapy 

2,042 

1,105  -  4,489 

3,136  (3,329) 

Other/  Unspecified 

3,344 

1,597  -  6,859 

5,770  (7,471) 

Abbreviations:  IQR  interquartile  range,  SD  standard  deviation 


Our  data  further  demonstrated  that  there  was  significant  geographic  variation  in  breast 
cancer-related  costs  (P<0.001)  (Table  8,  Figure  6).  In  particular,  there  appeared  to  be  lower 
cancer-related  spending  in  the  Southeastern  region  of  the  US  and  parts  of  the  Midwestern  US 
(Table  8)  (P<0.001).  This  geographic  variation  appeared  to  correlate  with  variation  in  overall 
healthcare  costs,  with  overall  spending  also  lowest  in  the  Southeastern  US. 

Other  predictors  were  also  found  to  have  a  significant  correlation  with  breast  cancer- 
related  costs  on  univariate  analysis.  Not  surprisingly,  treatment  utilization  with  any  treatment 
modality  (surgery,  chemotherapy,  or  radiotherapy)  was  associated  with  increased  costs.  In 
addition,  other  clinical  variables  showed  strong  associations,  including  indicators  of  greater 
disease  involvement  such  as  axillary  disease  (nodal  involvement).  Patient  factors  also  were 
associated,  with  younger  women  with  fewer  comorbidities  more  likely  to  incur  breast  cancer 
related  costs.  Interestingly,  however,  socioeconomic  variables  also  demonstrated  associations, 
with  patients  who  were  white,  had  higher  income  and  higher  education,  and  lived  in  a 
metropolitan  area  more  likely  to  incur  breast  cancer  related  costs  (Table  9).  These  variables 
remained  significant  on  multivariate  analysis. 
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Table  8.  Geographic  variation  in  breast  cancer-related  costs  and  overall  costs  (in  $). 

Breast  Cancer-Related  Overall 


States 

Median  ($) 

IQR  25th 

IQR  75th 

Median  ($) 

IQR  25th 

IQR  75th 

AK,  CA,  HI,  OR,  WA 

6,367 

2,993 

13,166 

13,079 

8,224 

19,452 

AZ,  CO,  ID,  MT,  NV,  NM,  UT, 

WY 

6,480 

3,108 

13,409 

11,875 

7,561 

17,971 

IA,  KS,  MN,  MO,  NE,  ND,  SD 

5,076 

2,636 

11,368 

10,885 

6,501 

16,094 

IL,  IN,  Ml,  OH,  Wl 

5,993 

2,971 

12,817 

12,255 

7,770 

18,719 

CT,  MA,  NH,  ME,  Rl,  VT 

6,025 

2,993 

12,933 

12,596 

8,616 

17,650 

NJ,  NY,  PA 

6,568 

2,962 

13,677 

13,474 

8,662 

20,748 

DE,  DC,  FL,  GA,  MD,  NC,  SC, 

VA,  WV 

6,300 

2,997 

13,797 

12,656 

7,837 

20,219 

AK,  LA,  OK,  TX 

6,541 

2,919 

13,112 

11,612 

7,069 

19,207 

AL,  KY,  MS,  TN 

5,018 

2,424 

11,502 

10,804 

6,304 

16,785 

|  8,800  to  10,000 
]  7,600  to  8,799 
]  6,400  to  7,599 
□  5,200  to  6,399 
]  4,000  to  5, 199 


Figure  6.  Breast-cancer  related  costs  across  the  United  States 
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Table  9.  Predictors  of  breast  cancer-related  spending 

Patient  Characteristic _ Median  Cost  ($) _ P _ 

Demographic 

Age 

66  to  <70  years  8,120  <0.001 

>70  years  5,670 

Race 

White  6,123  0.03 

Black  5,952 

Other  5,766 

Clinical 

Charlson  comorbidity  score 

0  comorbid  conditions  6,077  <0.001 

1  comorbid  condition  5,675 

2  or  more  comorbid  conditions  5,887 

Unknown  9,652 

Disease  Stage  and  Treatment 

Breast-conserving  surgery  7,132  <0.001 

Mastectomy  4,530 

Axillary  lymph  node  positive  disease  12,022  <0.001 

Lymph  node-negative  disease  5,514 

Axillary  lymph  node  dissection  7,564  <0.001 

No  axillary  dissection  4,885 

Post-mastectomy  reconstruction  8,890  <0.001 

No  post-mastectomy  reconstruction  6,037 

Radiation  therapy  10,087  <0.001 

No  radiation  therapy  3,208 

Chemotherapy  18,512  <0.001 

No  chemotherapy  5,012 

Healthcare  access* 

Surgeon  density*  (Pearson's  R)  0.00036  0.93 

Radiation  oncologist  density*  (Pearson's  R)  0.016  0.0002 

Socioeconomic  status* 

Living  in  metropolitan  area  6,326  <0.001 

Living  in  non-metropolitan  area  5,371 

Income  (Pearson's  R)  0.034  <0.001 

Percent  with  college  education  (Pearson's  R) _ 0.021 _ <0.001 

Abbreviations:  BCS  breast-conserving  surgery 
*  Defined  by  patient  county  of  residence 
+  Per  100,000  persons 
X  Per  10,000  persons 
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Brochytheropy  utilization  and  pilot  longitudinal  subgroup 

Of  all  patients  from  the  national  Medicare  dataset  treated  BCS,  97%  were  treated  with 
external  beam  radiotherapy  (EBRT)  alone,  3%  with  brachytherapy  alone,  and  <1%  with  EBRT 
plus  brachytherapy  boost.  For  patients  treated  with  brachytherapy,  98%  received  interstitial 
therapy  and  2%  intracavitary  therapy.  Though  percent  utilization  of  brachytherapy  modalities 
ranged  from  1%  to  4%  across  different  locations  in  the  US,  no  statistically  significant  variation 
was  detected  by  state  (P=0.62)  or  by  region  (P=0.32).  In  addition,  brachytherapy  use  did  not 
differ  by  race  (P=0.63)  or  age  (P=0.59). 

Of  the  entire  sample,  5%  (333  of  6,882)  received  brachytherapy  alone  (multi-catheter  or 
balloon-based),  95%  (6,521  of  6,882)  received  EBRT,  and  <1%  (28  of  6,865)  received  EBRT  plus 
brachytherapy  boost.  Treatments  with  brachytherapy  alone  significantly  increased  over  time, 
from  <1%  in  2001,  2%  in  2002,  3%  in  2003,  5%  in  2004,  8%  in  2005,  and  10%  in  2006  (P<0.001) 
(Table  10).  The  most  notable  increases  could  be  benchmarked  against  two  major  policy  events: 
First,  an  increase  in  use  was  noted  after  July  2002,  correlating  with  FDA  approval  of  the  balloon- 
based  breast  brachytherapy  device  (June  2002);  and  also,  a  further  increase  was  noted  after 
July  2004,  correlating  with  Medicare  reimbursement  of  treatment  (April  2004)  (Figure  7).  Of 
patients  treated  with  any  form  of  brachytherapy  alone,  the  proportion  who  received  balloon- 
based  treatment  also  increased  dramatically  over  time,  with  89%  receiving  balloon-based 
treatment  by  2006  (Figure  8).  In  multivariate  analysis,  the  temporal  trend  indicating  a  steady 
increase  in  the  use  of  brachytherapy  remained  significant  (P<0.001). 

Table  10.  Temporal  trends  in  brachytherapy  utilization 


Time  period  Total 

%  with  EBRT  only 

(N) 

%  with  EBRT  +  Boost 

(N) 

%  Brachytherapy  only 

(N) 

1/1/2001-6/30/2001  183 

98.91 

181 

0.00 

0 

1.09 

2 

7/1/2001  - 12/30/2001  363 

98.90 

359 

0.55 

2 

0.55 

2 

1/1/2002-6/30/2002  469 

99.36 

466 

0.21 

1 

0.43 

2 

7/1/2002  - 12/30/2002  494 

96.76 

478 

0.20 

1 

3.04 

15 

1/1/2003-6/30/2003  711 

97.19 

691 

0.14 

1 

2.67 

19 

7/1/2003  - 12/30/2003  828 

96.14 

796 

0.48 

4 

3.38 

28 

1/1/2004-6/30/2004  874 

94.62 

827 

0.34 

3 

5.03 

44 

7/1/2004  - 12/30/2004  851 

94.48 

804 

0.24 

2 

5.29 

45 

1/1/2005-6/30/2005  942 

91.83 

865 

0.85 

8 

7.32 

69 

7/1/2005  - 12/30/2005  534 

91.57 

489 

0.37 

2 

8.05 

43 

1/1/2006-6/30/2006  633 

89.26 

565 

0.63 

4 

10.11 

64 

Total  6882 

94.75 

6521 

0.41 

28 

4.83 

333 
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%  Brachytherapy  Alone 


Figure  7.  Temporal  trends:  Percent  of  patients  treated  with  brachytherapy  (multi-catheter  or 
balloon-based)  as  the  sole  modality  of  radiotherapy  after  breast-conserving  therapy.  Year  'a' 
refers  to  January  through  June  and  'b'  July  through  December.  (P<0.001) 


2001a  2001b  2002a  2002b  2003a  2003b  2004a  2004b  2005a  2005b  2006a 


Time  Period 
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Figure  8.  Temporal  trends:  Percent  of  patients  treated  with  brachytherapy  as  the  sole 
modality  of  radiotherapy  after  breast-conserving  therapy  that  received  balloon-based 
brachytherapy.  Year  'a'  refers  to  January  through  June  and  'b'  July  through  December. 
(P<0.001) 
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Part  IV 


Survival.  In  the  longitudinal  cohort  of  women  treated  with  BCS,  of  62,625  patients,  by 
the  end  of  the  follow-up  period,  a  total  of  11%  of  patients  had  died  (n=6,705).  For  those  who 
died,  the  median  time  to  death  was  709  days  (IQR  445,  1049).  Survival  was  significantly 
associated  with  the  use  of  RT  after  BCS.  Sixty-four  percent  of  all  patients  received  RT  over  the 
study  period.  Of  patients  who  received  RT,  17%  died,  compared  with  7%  of  those  who  did  not 
receive  RT  (P<0.001).  Survival  rates  varied  significantly  by  state  (P<0.001)  and  geographic 
region  (P<0.001,  Figure  9). 

Figure  9.  Percent  died.  Darker  shading  represents  a  higher  percentage  of  patients  who  died. 
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Notably,  high  death  rates  (and  thus  low  survival  rates)  were  still  found  in  areas  of  the  South  and 
Midwest,  areas  which  had  also  previously  demonstrated  lower  rates  of  appropriate  treatment 
utilization.  The  Northeast  also  showed  higher  death  rates  on  unadjusted  analysis. 

On  multivariate  analysis,  younger  age  and  white  race  were  associated  with  better 
overall  survival.  In  addition,  use  of  RT  was  associated  with  better  survival  (OR=  0.50,  95%  Cl 
0.48,  0.53).  Interestingly,  even  after  adjusting  for  these  covariate  factors,  geographic  region 
remained  a  significant  predictor  (Table  11). 

Table  11.  Adjusted  and  unadjusted  survival  by  region 


Region 

States 

%  Died 

OR 

95%  Cl 

West,  Pacific  West 

AK,  CA,  HI,  OR,  WA 

11 

- 

West,  Mountain  West 

AZ,  CO,  ID,  MT,  NV,  NM,  UT,  WY 

10 

1.23 

1.06 

1.41 

Midwest,  West  North  Central 

IA,  KS,  MN,  MO,  NE,  ND,  SD 

10 

1.18 

1.04 

1.34 

Midwest,  East  North  Central 

IL,  IN,  Ml,  OH,  Wl 

12 

1.41 

1.28 

1.55 

Northeast,  New  England 

CT,  MA,  NH,  ME,  Rl,  VT 

11 

1.30 

1.15 

1.47 

Northeast,  Mid-Atlantic 

NJ,  NY,  PA 

12 

1.28 

1.16 

1.41 

South,  South  Atlantic 

DE,  DC,  FL,  GA,  MD,  NC,  SC,  VA,  WV 

11 

1.26 

1.15 

1.39 

South,  West  South  Central 

AK,  LA,  OK,  TX 

11 

1.28 

1.13 

1.45 

South,  East  South  Central 

AL,  KY,  MS,  TN 

12 

1.43 

1.24 

1.64 
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Discussion 


Summary,  Comments,  and  Future  Work 

In  this  contemporary,  national  cohort  of  older  breast  cancer  patients,  we  demonstrated 
the  following  novel  findings: 

•  Significant  geographic  (regional)  variation  existed  for  the  use  of  chemotherapy  and 
RT.  Furthermore,  even  though  multivariate  adjustment  attempted  to  adjust  for 
variations  in  case-mix  across  the  United  States,  geographic  variation  was  still  a 
significant  predictor  of  treatment  choices. 

•  For  the  use  of  RT  after  BCS,  a  significant  amount  of  variation  also  existed.  The 
absolute  frequency  of  underutilization  of  RT  was  also  surprising,  as  high  as  27% 
across  the  country  and  even  as  high  as  50%  in  some  states.  Racial  differences  in 
utilization  were  also  significant. 

•  The  use  of  breast  conservation  for  early  invasive  disease  also  significantly  varied. 

Use  of  BCS  was  as  low  as  approximately  50%  in  some  regions  of  the  United  States, 
and  was  associated  with  socioeconomic  factors,  suggesting  that  socioeconomic 
barriers  have  a  significant  effect  on  the  use  of  breast  conservation. 

•  Clinical  and  non-clinical  variations  in  breast  cancer-related  costs  were  also  of 
interest.  Increased  spending  appeared  to  reflect  adherence  to  standard  treatment 
(for  example,  the  use  of  RT  after  BCS)  and  higher  socioeconomic  status. 

•  Use  of  newer  therapies  was  also  associated  with  both  clinical  and  non-clinical 
factors,  with  both  patient  and  physician  factors  contributing.  Interestingly,  some 
socioeconomic  and  policy  factors  may  overshadow  clinical  efficacy  data  in 
influencing  physician  or  patient  behavior  in  selecting  treatments  for  breast  cancer. 

•  Survival  was  associated  with  use  of  RT  after  BCS.  Flowever,  survival  varied 
significantly  by  region.  While  regional  patterns  resembled  utilization  patterns,  still 
utilization  (or  underutilization)  of  standard  care  could  not  fully  explain  differences  in 
survival  across  the  United  States. 

Please  see  Annual  Reports,  September  2008  and  September  2009,  for  any  further  detailed 
discussion  of  these  points. 

The  body  of  work  that  we  have  produced  from  our  analyses  has  provided  a  unique 
contribution  to  the  literature.  Strengths  of  our  project  include  the  use  of  a  novel  dataset,  the 
national  Medicare  dataset.  This  is  one  of  the  largest  cohorts  assembled  to  date.  Our  long-term 
goal  of  obtaining  near  and  long-term  follow-up  allows  our  dataset  to  be  one  of  the  most  unique 
for  studying  the  course  of  both  treatment  and  outcomes  in  breast  cancer  patients.  Finally,  we 
contributed  to  gaps  in  the  literature  on  disparities  in  breast  cancer  care.  We  demonstrated  and 
quantified  persistent  disparities,  particularly  for  the  use  of  surgery  and  radiation  therapy  in 
breast  cancer  patients.  In  addition,  however,  using  one  of  the  unique  strengths  of  the  Medicare 
dataset— the  comprehensive  information  provided  on  all  beneficiaries  across  the  United 
States— we  were  able  to  study  in  detail  the  geographic  variations  of  the  disparities.  The 
significance  and  the  absolute  magnitude  of  these  variations  were  staggering. 
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Our  results  prompt  several  areas  of  future  directions  and  needs  in  the  area  of 
investigating  disparities  in  breast  cancer  care.  While  we  were  able  to  detail  the  magnitude  of 
disparities  and  regional  variation  in  treatment,  cost,  and  outcomes  using  our  Medicare  dataset, 
future  investigations  may  focus  on  what  are  the  underlying  causes  of  regional  differences.  We 
hypothesize,  in  part  prompted  by  our  hypothesis-generating  results,  that  both  physician  and 
patient  influences- and  consideration  of  the  patient-physician  interaction-  may  be 
contributory.  Of  importance  is  understanding  these  potential  factors  and  being  able  to  take  into 
account  individual,  regional,  or  cultural  preferences  for  care  when  determining  treatment  goals. 
Another  important  question  that  our  results  prompt  is,  what  are  the  best  strategies  that  may 
be  applied  to  overcome  real  socioeconomic  barriers  that  clearly  continue  to  affect  treatment 
utilization  across  the  United  States?  Our  finding  of  socioeconomic  barriers  to  breast 
conservation  in  our  contemporary  cohort  suggest  that  little  progress  has  been  made  in 
overcoming  disparities  in  the  use  of  breast  conserving  therapy,  particularly  for  patients  living  in 
the  most  disadvantaged  neighborhoods  and  regions.  Future  research  may  explore  whether 
interventions  at  the  individual  patient  level,  combined  with  changes  at  the  policy  level,  may 
help  to  overcome  barriers  to  quality  breast  cancer  care.  Finally,  our  investigations  have  focused 
on  older  patients,  and  these  data  should  be  validated  in  populations  of  younger  breast  cancer 
patients  as  well.  The  implications  of  answering  these  questions  are  on  the  ever-present 
demands  to  simultaneously  optimize  treatment  quality,  effectiveness  and  costs.  It  is  hoped  that 
the  insights  gained  in  this  area  of  research  may  help  to  improve  the  rationale  by  which 
therapies  are  promoted  and  adopted  into  standard  care. 

KEY  RESEARCH  ACCOMPLISHMENTS 

•  Assembled  a  unique,  contemporary  cross-sectional  Medicare  breast  cancer  cohort 

•  Assembed  a  unique,  longitudinal  Medicare  breast  cancer  cohort  with  follow-up 

•  Obtained  and  analyzed  treatment  utilization  in  several  validation  populations,  including 
SEER-Medicare  and  MarketScan  databases. 

•  Developed,  applied,  and  validated  unique  algorithms  for  identifying  breast  cancer  patients 
and  early  stage  subgroup  of  patients  based  on  claims  data. 

•  Conducted  univariate  and  multivariate  data  analyses  on  treatment  utilization,  cost,  and 
outcomes  in  breast  cancer  patients. 

•  Further  conducted  stratified  and  subsidiary  analyses,  including  an  analysis  of  the 
disparities  in  treatment  and  costs  by  clinical  and  non-clinical  factors 

•  Publication  of  multiple  manuscripts  and  preparation  of  manuscripts  in  progress,  at  peer- 
reviewed  journals 

•  Presentation  of  results  at  national  scientific  meetings 
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CONCLUSION 


At  the  conclusion  of  Year  3  of  our  research,  we  have  made  significant  progress  toward 
accomplishing  our  project  goals.  Specifically,  we  have  worked  to  accomplish  the  objectives 
stated  for  Year  3  in  the  Statement  of  Work.  Using  our  novel,  comprehensive  national  Medicare 
dataset,  we  conducted  several  retrospective  analyses  on  a  cohort  of  older  women  diagnosed 
with  early  invasive  breast  cancer.  Results  from  our  analyses  provided  novel  insights  that 
contribute  to  the  existing  scientific  literature.  The  most  striking  results  from  our  analyses 
suggest  that  variation  in  breast  cancer  care  is  significant,  and  these  variations  appear  to 
contribute  to  variations  in  costs  and  outcomes  for  breast  cancer  care  across  the  United  States. 
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Original  Article 


Racial  Disparities  in  the  Use  of  Radiotherapy 
After  Breast-Conserving  Surgery 

A  National  Medicare  Study 

Grace  L.  Smith,  MD,  PhD,  MPH1;  Ya-Chen  T.  Shih,  PhD2;  Ying  Xu,  MD,  MS2; 

Sharon  H.  Giordano,  MD,  MPH3;  Benjamin  D.  Smith,  MD1,4;  George  H.  Perkins,  MD1; 

Welela  Tereffe,  MD1;  Wendy  A.  Woodward,  MD,  PhD1;  and  Thomas  A.  Buchholz,  MD1 


BACKGROUND:  In  prior  studies,  the  use  of  standard  breast  cancer  treatments  has  varied  by  race,  but  previous  analy¬ 
ses  were  not  nationally  representative.  Therefore,  in  a  comprehensive,  national  cohort  of  Medicare  patients,  racial  dis¬ 
parities  in  the  use  of  radiotherapy  (RT)  after  breast-conserving  surgery  (BCS)  for  invasive  breast  cancer  were 
quantified.  METHODS:  A  national  Medicare  database  was  used  to  identify  all  beneficiaries  (age  >65  years)  treated 
with  BCS  for  incident  invasive  breast  cancer  in  2003.  Claims  codes  identified  RT  use,  and  Medicare  demographic 
data  indicated  race.  Logistic  regression  modeled  RT  use  in  white,  black,  and  other-race  patients,  adjusted  for  demo¬ 
graphic,  clinical,  and  socioeconomic  covariates.  RESULTS:  Of  34,080  women,  91%  were  white,  6%  were  black,  and  3% 
were  another  race.  The  mean  age  of  the  patients  was  76  ±  7  years.  Approximately  74%  of  whites,  65%  of  blacks,  and 
66%  of  other-race  patients  received  RT  (P  <  .001).  After  covariate  adjustment,  whites  were  found  to  be  significantly 
more  likely  to  receive  RT  than  blacks  (odds  ratio,  1.48;  95%  confidence  interval,  1.34-1.63  [P  <  .001]).  Disparities 
between  white  and  black  patients  varied  by  geographic  region,  with  blacks  in  areas  of  the  northeastern  and  southern 
United  States  demonstrating  the  lowest  rates  of  RT  use  (57%  in  these  regions).  In  patients  age  <70  years,  racial  dis¬ 
parities  persisted.  Specifically,  83%  of  whites,  73%  of  blacks,  and  78%  of  other  races  in  this  younger  group  received 
RT  (P  <  .001).  CONCLUSIONS:  In  this  comprehensive  national  sample  of  older  breast  cancer  patients,  substantial 
racial  disparities  were  identified  in  RT  use  after  BCS  across  much  of  the  United  States.  Efforts  to  improve  breast 
cancer  care  require  overcoming  these  disparities,  which  exist  on  a  national  scale.  Cancer  2010:116:734-41.  ©  2009 
American  Cancer  Society. 

KEYWORDS:  breast  cancer,  Medicare,  radiotherapy,  race,  disparities. 


Postoperative  radiotherapy  (RT)  is  generally  considered  a  standard  component  of  therapy  for  patients  diagnosed 
with  early  invasive  breast  cancer  who  undergo  breast-conserving  surgery  (BCS).  Cumulative  results  from  randomized  trials 
have  demonstrated  local  control  and  survival  benefits  with  the  addition  of  RT  after  BCS. 11  However,  despite  the  evidence 
supporting  the  use  of  RT  after  BCS,  a  recent  analysis  of  treatment  patterns  in  breast  cancer  patients  enrolled  in  the  North 
American  Fareston  versus  Tamoxifen  Adjuvant  (NAFTA)  trial  suggested  that  racial  disparities  may  exist  in  RT  use  after 
BCS.  Specifically,  this  study  found  that  postoperative  RT  was  administered  in  82%  of  white  patients  but  only  70%  of 
black  patients,  although  the  difference  was  not  found  to  be  statistically  significant.5 

Other  studies  of  breast  cancer  patients  have  also  suggested  racial  disparities  in  the  receipt  of  standard  treatments  for 
locally  advanced  and  regional  disease,  including  differences  among  white  and  nonwhite  patients  in  the  use  of  RT,  surgery, 
and  chemotherapy.5"15  Results  from  many  of  these  studies  suggest  that  nonwhite  patients  may  more  frequently  encounter 
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barriers  to  receiving  standard  care.  However,  prior  analy¬ 
ses  were  conducted  on  limited  samples,  including  single¬ 
state,  single-center,  or  clinical  trial  patient  cohorts,  and 
studies  of  population-based  cohorts  have  not  to  our 
knowledge  been  nationally  comprehensive.  Conse¬ 
quently,  the  true  magnitude  of  and  regional  variation  in 
racial  disparities  in  breast  cancer  treatment  across  the 
United  States  remain  uncertain. 

Quantifying  the  magnitude  of  racial  disparities  in 
the  use  of  standard  treatments  across  the  nation  is  an  im¬ 
portant  step  toward  identifying  and  ultimately  reducing 
barriers  to  high-quality  breast  cancer  care.  Therefore, 
using  data  from  a  comprehensive,  national  sample  of 
Medicare  patients,  we  sought  to  quantify  racial  differences 
in  the  use  of  RT  after  BCS  for  women  with  invasive  breast 
cancer.  In  addition,  we  sought  to  determine  whether  there 
was  geographic  variation  across  the  United  States  in  the 
association  between  race  and  RT  use. 

MATERIALS  AND  METHODS 

National  Medicare  Data  Set  and  Study 
Sample  Derivation 

The  national  Medicare  data  set  includes  comprehensive 
claims  data  for  all  Medicare  beneficiaries  in  the  United 
States.  Files  contain  data  collected  by  Medicare  for  reim¬ 
bursement  of  healthcare  services  for  each  beneficiary  and 
include  institutional  (inpatient  and  outpatient)  and  non- 
institutional  (physicians  or  other  providers)  final  action 
claims.10  The  University  of  Texas  M.  D.  Anderson  Can¬ 
cer  Center’s  institutional  review  board  approved  use  of 
the  Medicare  database  for  this  study. 

We  used  the  following  algorithm  to  identify  patients 
with  incident  breast  cancer,  which  we  based  on  a  prior 
validated  algorithm  using  claims  data.17’18  Our  initial 
study  population  was  comprised  of  853,273  women  who 
had  any  diagnosis  of  invasive  breast  cancer  in  2003, 
defined  as  an  International  Classification  of  Diseases- 
Ninth  Revision  (ICD-9)  diagnosis  code  of  174.  Because 
this  denominator  would  have  included  incident  and  prev¬ 
alent  cases  in  2003,  as  well  as  potentially  disease-free 
patients  undergoing  breast  cancer  screening,  we  applied 
the  following  algorithm  to  refine  our  study  sample  to 
identify  incident  cases  only.  From  our  initial  study  popu¬ 
lation,  we  included  83,611  patients  aged  >65  years  who 
underwent  BCS  between  January  1,  2003,  and  December 
31,  2003.  From  this  sample,  we  excluded  17,872  patients 
who  did  not  have  at  least  2  claims  (on  different  dates) 
specifying  a  diagnosis  of  invasive  breast  cancer  between 
January  1,  2003,  and  December  31,  2004  (at  least  1  claim 


must  have  occurred  during  2003),  with  no  more  than  6 
months  between  the  date  of  BCS  and  the  earliest  breast 
cancer  diagnosis  claim  date.  To  exclude  the  prevalent 
cases,  we  excluded  10,362  patients  who  had  a  breast  can¬ 
cer-related  diagnosis  or  procedure  claim  between  January 
1,  2002,  and  December  31,  2002.  To  ensure  that  BCS 
was  intended  to  be  the  primary  cancer-directed  surgery, 
we  excluded  13,470  patients  who  underwent  mastectomy 
within  3  months  of  BCS.  To  limit  our  sample  to  patients 
with  early  stage  invasive  breast  cancer,  we  then  excluded 
1082  patients  who  had  >2  claims  specifying  metastatic 
breast  cancer  from  3  months  before  to  3  months  after  the 
diagnosis  date.  To  improve  sample  homogeneity,  we  also 
excluded  2964  patients  who  were  receiving  Medicare  cov¬ 
erage  because  of  end-stage  renal  disease  or  disability. 
Finally,  to  ensure  we  had  complete  claims  information  to 
determine  patients’  cancer  treatment  course  and  comor¬ 
bidities,  we  excluded  3781  patients  who  lacked  Medicare 
Part  A  or  B  coverage  or  who  had  intermittent  health  main¬ 
tenance  organization  coverage  in  the  9  months  after  or  in 
the  1  year  before  their  breast  cancer  diagnosis  date  (of 
these  patients,  2143  had  incomplete  information  in  the 
year  before  diagnosis  because  they  were  aged  <66  years). 
For  our  study,  the  breast  cancer  diagnosis  date  was  consid¬ 
ered  the  date  of  the  earliest  claim  for  a  diagnosis  of  breast 
cancer.  This  left  a  final  sample  size  of  34,080  patients. 
The  complete  list  of  claims  codes  applied  in  our  algorithm 
can  be  found  in  Table  1 . 

RT  and  Race 

Patients  were  considered  to  have  received  breast  RT  if  a  claim 
for  RT  (Appendix  A)  occurred  within  9  months  of  the  breast 
cancer  diagnosis  date.  This  definition  has  been  validated  in 
prior  studies  of  breast  cancer  patients.1 9-24  In  addition,  for 
patients  classified  as  receiving  RT,  completion  of  the  thera¬ 
peutic  course  was  defined  as  at  least  3  complete  weeks  of 
therapy  (5  treatments  per  week).  Patient  race  was  determined 
using  the  Medicare  denominator  file,  which  contains  demo¬ 
graphic  information  on  each  beneficiary.  In  this  file,  race  and 
ethnicity  data  are  based  on  the  patients’  self-report.  For  our 
analysis,  we  categorized  patients  as  white,  black,  and  other 
(nonwhite,  nonblack)  race  based  on  a  prior  study  of  the  sen¬ 
sitivity  and  specificity  of  Medicare  race  and  ethnicity 
25 

groupings. 

Other  Covariates 

Demographic  and  clinical  covariates  were  derived  from 
Medicare  files,  including  the  denominator  file  and  claims 
files.  Demographic  data  included  age  at  diagnosis  and 
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Table  1.  Claims  Codes 

Variable 

Breast-conserving  surgery 

Breast  cancer-related  diagnosis  and  procedure  claims 
Biopsy 

Breast-conserving  surgery 
Mastectomy 
Lymph  node  dissection 
Radiotherapy 
History  of  breast  cancer 


ICD-9  Code 

85.20,  85.21,  85.22,  85.23 


85.1- 85.19 
As  above 
85.33-85.48 
40.3 

92.2- 92.29,  0330* *,  0333* 
VI  0.3 


CPT  Code 

19120,  19125,  19126,  19160,  19162 


19000,  19001,  19100,  19101,  19110,  19112 

As  above 

19180-19255 

38525,  38740,  38745 

77400-77499,  77520-77525,  77750-77799 


ICD-9  indicates  International  Classification  of  Diseases,  Ninth  Revision;  CPT,  Current  Procedural  Terminology. 

*  Revenue  Center  code. 


state  of  residence,  with  geographic  regions  based  on  the 
Census  Divisions.26  Disease-related  and  treatment-related 
variables  included  axillary  lymph  node  involvement;  axil¬ 
lary  lymph  node  dissection;  sentinel  lymph  node  biopsy; 
receipt  of  any  chemotherapy;  receipt  of  doxorubicin  or 
paclitaxel;  any  staging  imaging;  number  of  hospitaliza¬ 
tions  in  the  year  after  diagnosis;  and  number  of  medical 
oncology,  radiation  oncology,  and  surgery  visits  in  the 
year  after  diagnosis.  Variables  indicating  preventive 
healthcare  and  interactions  with  the  healthcare  system 
included  mammography  in  the  year  before  diagnosis  and 
number  of  physician  visits  in  the  year  before  diagnosis.  In 
addition,  we  calculated  the  severity  of  comorbid  disease 
for  each  patient  based  on  a  modified  Charlson  comorbid¬ 
ity  score  validated  in  a  prior  claims-based  study:  0  (no 
comorbidity),  1  (mild  to  moderate),  or  >2  (severe).27 
This  score  combined  comorbidities  recorded  in  Medicare 
claims  during  the  12  months  before  the  patient’s  cancer 
diagnosis.  To  enhance  the  specificity  of  comorbid  disease 
diagnoses,  patients  must  have  had  at  least  1  inpatient 
(Part  A)  claim  or  at  least  2  outpatient  (Part  B)  claims 
more  than  30  days  apart.27  Socioeconomic  covariates 
were  derived  from  the  2003  Area  Resource  File28  linked 
to  the  Medicare  files  by  county  and  state,  and  included  ru¬ 
ral/ urban  status,  percentage  of  population  (by  county)  liv¬ 
ing  in  poverty,  median  income,  education  level,  density  of 
surgeons,  and  density  of  radiation  oncologists. 

Statistical  Analysis 

All  analyses  were  conducted  using  SAS  software  (version 
9.1.3;  SAS  Institute  Inc,  Cary,  NC),  and  all  statistical  tests 
assumed  a  2-tailed  a  of  0.05.  We  calculated  percentage  RT 
use  for  the  entire  sample,  by  state  and  by  region,  and  tested 
the  unadjusted  association  between  receipt  of  RT  and  race 
using  the  Pearson  chi-square  test.  Bivariate  associations 
between  receipt  of  RT  and  other  covariates  were  tested 


using  the  Pearson  chi-square  test  for  categoric  variables  and 
the  Wilcoxon  rank  sum  test  for  continuous  variables.  A 
multivariate  logistic  model  tested  the  adjusted  association 
between  receipt  of  RT  and  race.  Covariates  were  selected  a 
priori  based  on  significance  in  bivariate  analyses  (P  <  .25) 
and  significance  in  prior  studies  of  cancer  patients.17"24 
Goodness  of  fit  was  assessed  using  the  Hosmer  and  Leme- 
show  test.  The  final  parsimonious  multivariate  logistic 
model  was  then  also  used  to  calculate  adjusted  rates  of  RT 
use.27  A  subsidiary  model  included  socioeconomic  covari¬ 
ates  and  tested  the  association  in  the  subset  of  33,172 
patients  who  had  socioeconomic  covariate  data  available. 

We  further  examined  the  association  between 
receipt  of  RT  and  race  in  the  subset  of  patients  age  <70 
years  (n  —  7270),  because  prior  evidence  has  suggested 
that  the  use  of  RT  after  BCS  is  more  common  in  younger 
patients.5  In  sensitivity  analyses,  we  tested  the  association 
between  race  and  receipt  of  RT  in  the  subset  of  patients 
with  a  more  restrictive  (specific)  definition  of  invasive 
breast  cancer,  which  required  the  2  diagnosis  codes  for 
invasive  breast  cancer  to  be  >30  days  apart  (n  =  33,1 14), 
and  in  a  more  restrictive  definition  of  breast  cancer,  which 
excluded  patients  with  any  diagnosis  code  of  ductal  carci¬ 
noma  in  situ  (DCIS)  (ICD-9  diagnosis  code  233.0) 
within  6  months  of  the  invasive  breast  cancer  diagnosis 
date  (n  =  22,301).  We  also  tested  a  less  restrictive  treat¬ 
ment  window  of  RT  within  12  months  of  breast  cancer 
diagnosis  to  take  into  account  possible  racial  differences  in 
treatment  delays.  Finally,  in  patients  (of  any  age)  who 
received  RT,  we  also  examined  the  association  between 
completion  of  RT  course  and  race. 

RESULTS 

Patient  Characteristics 

Of  the  34,080  women  with  invasive  breast  cancer  who 
were  treated  with  BCS,  91%  (N  =  31,127)  were  white, 
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6%  (N  =  2077)  were  black,  and  3%  (N  =  876)  were 
another  race.  The  mean  age  of  the  sample  was  76  (stand¬ 
ard  deviation,  7)  years.  Approximately  73%  of  patients 
received  RT  after  BCS,  and  13%  received  chemotherapy 
as  part  of  their  initial  treatment  course.  Approxmately 
39%  of  patients  had  a  claim  indicating  axillary  lymph 
node  dissection  or  axillary  lymph  node  involvement,  and 
55%  had  a  claim  indicating  sentinel  lymph  node  biopsy 
(Table  2). 

Race  and  RT  Use 

Significant  racial  differences  were  found  to  exist  in  the 
receipt  of  RT;  74%  of  white  women,  65%  of  black 
women,  and  66%  of  other-race  women  received  RT  after 
BCS  (P  <  .001).  After  adjustment  for  demographic,  clin¬ 
ical,  and  socioeconomic  covariates,  white  women  were 
still  significantly  more  likely  than  black  women  to  have 
received  RT  after  BCS  (odds  ratio  [OR],  1.48;  95%  confi¬ 
dence  interval  [95%  Cl],  1.34-1.63  [ P  <  .001]).  In  addi¬ 
tion,  there  was  a  trend  toward  white  women  also  being 
more  likely  than  other-race  women  to  have  received  RT 
after  BCS  (OR,  1.22;  95%  Cl,  1.04-1.42  [P  =  .01]) 
(Table  3).  Other  variables  found  to  be  independently 
associated  with  higher  rates  of  RT  use  included  geo¬ 
graphic  region,  younger  age,  lower  comorbidity  score, 
treatment  with  doxorubicin,  axillary  lymph  node  dissec¬ 
tion,  use  of  staging  imaging,  use  of  mammography  before 
cancer  diagnosis,  and  greater  number  of  physician  visits 
(interactions  with  the  healthcare  system)  in  the  year  before 
cancer  diagnosis.  In  the  subset  of  patients  with  complete 
socioeconomic  covariate  data  available,  the  association 
between  receipt  of  RT  and  race  was  unchanged.  In  this 
model,  lower  poverty  level,  higher  median  income,  higher 
education  level,  higher  surgeon  density,  and  higher  radia¬ 
tion  oncologist  density  were  also  found  to  be  associated 
with  receipt  of  RT. 

Sensitivity  analyses  using  more  restrictive  definitions 
of  invasive  breast  cancer  did  not  appear  to  alter  the  magni¬ 
tude  or  significance  of  the  associations  between  race  and 
RT  use  significantly.  Sensitivity  analyses  using  a  less  re¬ 
strictive  treatment  window  also  did  not  appear  to  alter  the 
magnitude  or  significance  of  these  associations. 

In  the  subset  of  patients  age  <70  years,  82% 
received  RT  after  BCS.  Specifically,  83%  of  white 
women,  73%  of  black  women,  and  78%  of  other-race 
women  received  RT  after  BCS  (P  <  .001).  On  multivari¬ 
ate  analysis  of  this  younger  group,  white  women  were  still 
found  to  be  significantly  more  likely  than  black  women  to 
have  received  RT  after  BCS  (OR,  1.73;  95%  Cl,  1.40- 


Table  2.  Sample  Characteristics 

Characteristic  Percentage 


Demographic 

Mean  age  (SD),  y  76  (7) 

Race 

White  91 

Black  6 

Other  3 

Region 

West,  Pacific  West  1 1 

West,  Mountain  West  5 

Midwest,  West  North  Central  7 

Midwest,  East  North  Central  19 

Northeast,  New  England  7 

Northeast,  Mid-Atlantic  16 

South,  South  Atlantic  21 

South,  West  South  Central  9 

South,  East  South  Central  5 

Not  known  1 

Clinical 

No  comorbid  diseases  32 

Axillary  lymph  node  involvement  39 

Sentinel  lymph  node  biopsy  55 

Radiotherapy  73 

Chemotherapy  1 3 

Socioeconomic 

Mean  income  (SD),  $US  43,429  (10,649) 

Education 

Mean  percent  completed  less  7  (4) 

than  ninth  grade  (SD) 

Mean  percent  completed  24  (9) 

college  or  higher  (SD) 

Mean  surgeon  density  (per  100,000)  (SD)  13  (10) 

Mean  radiation  oncologist  density  1.5  (1.4) 

(per  100,000)  (SD) 


SD  indicates  standard  deviation. 


2.14  [P=  .003]),  although  there  was  no  significant  differ¬ 
ence  noted  between  white  women  and  other-race  women 
(OR,  1.28;  95%  Cl,  0.86-1.90  [P  —  .22]).  The  interac¬ 
tion  term  between  race  (black  vs  white  race)  and  age  (age 
<70  years  vs  >70  years)  was  not  found  to  be  statistically 
significant  ( P=  .19). 

In  the  subset  of  patients  who  received  RT  and  had 
the  number  of  treatments  documented,  85%  were  found 
to  have  completed  a  course  of  therapy,  whereas  15%  had 
an  incomplete  course.  Approximately  85%  of  white 
women,  85%  of  black  women,  and  82%  of  other-race 
women  completed  their  RT  course  (P  =  .02).  After 
adjusting  for  covariates,  there  remained  no  significant  dif¬ 
ference  in  rates  of  completion  by  race  (black  vs  white 
women:  OR,  1.00;  95%  Cl,  0.83-1.20  [P  =  .99]  and 
other-race  vs  white  women:  OR,  1.02;  95%  Cl,  0.73-1.43 
[P—  .89]). 
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Table  3.  Racial  Disparities  in  the  Use  of  Radiotherapy  After 
Breast-Conserving  Surgery  for  Breast  Cancer 


Model:  RT  use 

OR 

95%  Cl 

P 

Unadjusted 

White  vs  black  patients 

1.52 

1.38-167 

<.001 

White  vs  other*  patients 

1.46 

1.27-1.68 

<.001 

Adjusted! 

White  vs  black  patients 

1.48 

1.34-1.63 

<.001 

White  vs  other*  patients 

1.22 

1.04-1.42 

.01 

*RT  indicates  radiotherapy;  OR, 

odds  ratio; 

95%  Cl,  95% 

confidence 

interval. 

f  Adjusted  for  age,  comorbidity,  chemotherapy  (doxorubicin  or  paclitaxel), 
axillary  lymph  node  involvement,  staging,  imaging,  surgeon  visits,  mam¬ 
mography,  physician  visits,  and  region. 

Geographic  Variation  in  Racial  Disparities 

We  observed  significant  geographic  variations  in  racial 
disparities.  The  absolute  difference  in  the  percentage  of 
RT  use  between  white  and  black  women  by  state  ranged 
from  —28%  to  74%.  However,  for  most  states,  the  per¬ 
centage  of  white  women  who  received  RT  was  greater 
than  the  percentage  of  black  women  who  received  RT, 
and  in  no  state  was  the  rate  of  RT  use  found  to  be  signifi¬ 
cantly  greater  in  black  women  compared  with  white 
women  (Fig.  1).  In  areas  of  the  northeastern  and  southern 
United  States,  the  racial  disparity  between  white  and  black 
women  with  regard  to  RT  use  was  particularly  pro¬ 
nounced  (Fig.  2)  (Table  4).  Our  data  were  insufficient  to 
assess  variations  by  state  or  region  for  other-race  women. 

DISCUSSION 

In  our  analysis  of  a  national  sample  of  US  women  aged 
>65  years  who  were  diagnosed  with  invasive  breast  can¬ 
cer,  white  women  (74%)  were  significantly  more  likely 
than  black  women  (65%)  to  receive  RT  after  BCS.  The 
results  of  the  current  study  also  suggested  that  white 
patients  were  also  more  likely  than  patients  of  other  races 
(66%)  to  receive  RT  after  BCS.  Furthermore,  even  after 
we  took  into  consideration  other  covariates  such  as  patient 
comorbidity,  white  women  were  still  nearly  50%  more 
likely  than  black  women  to  receive  RT.  The  racial  dispar¬ 
ity  between  whites  and  blacks  persisted  even  among  the 
younger  patients  in  the  cohort,  those  age  <70  years,  in 
whom  RT  after  BCS  was  more  common.  Our  analysis 
also  identified  geographic  variations  in  the  magnitude  of 
treatment  disparities  between  white  and  black  patients 
across  the  United  States,  although  some  degree  of  treat¬ 
ment  disparity  existed  in  most  states.  Furthermore,  the 
disparities  across  the  nation  persisted  even  after  adjusting 
for  variations  in  socioeconomic  factors  such  as  income, 


education,  and  access  to  healthcare  resources.  The  large 
magnitude  of  racial  disparities  in  RT  use  identified  in  the 
current  study  is  concerning  given  that  the  administration 
of  RT  after  BCS  in  patients  with  early  invasive  breast 
cancer  reduces  the  risk  of  disease  recurrence  and  breast 

cancer-related  death  and  is  considered  the  standard  of 

1-4,30 

care. 

Other  studies  have  also  described  disparities  in 
breast  cancer  care  between  black  patients  and  white 
patients.  5  In  a  recent  secondary  analysis  of  the  NAFTA 
trial,  in  which  all  patients  in  the  study  cohort  were  treated 
with  BCS  and  hormonal  therapy,  80%  of  all  enrollees 
received  RT.  Similar  to  the  9%  to  10%  absolute  differ¬ 
ence  noted  in  the  current  study,  there  was  a  12%  absolute 
difference  found  between  white  and  black  patients  in  the 
rate  of  RT  use  in  the  NAFTA  study.5  A  recent  analysis  of 
the  Surveillance,  Epidemiology,  and  End  Results  (SEER)- 
Medicare  cohort  reported  an  8%  absolute  difference 
between  white  and  black  patients  with  regard  to  the 
receipt  of  RT  after  BCS,  with  86%  of  white  women  and 
78%  of  black  women  between  the  ages  of  66  and  85  years 
with  invasive  breast  cancer  receiving  the  treatment.31  A  li¬ 
mitation  of  this  study  was  that  SEER-Medicare  sampled 
<  26%  of  the  US  Medicare  population,  compared  with 
the  current  study,  which  sampled  the  Medicare  popula¬ 
tion  in  its  entirety. 

A  previous  study  of  breast  cancer  patients  who  were 
treated  in  the  state  of  Florida  found  that  Hispanic  women 
were  significantly  less  likely  to  receive  local  therapy  for 
breast  cancer  than  non-Hispanic  white  women,14  whereas 
another  study  of  patients  in  Hawaii  found  that  Asian 
women  may  be  less  likely  to  receive  therapy  than  white 
women.  3-  In  our  cross-sectional  study,  we  were  unable  to 
explore  the  use  of  RT  in  specific  nonwhite,  nonblack  eth¬ 
nic  groups  because  of  the  limited  sample  size  of  each  eth¬ 
nic  group  in  our  sample.  Future  analyses  may  seek  to 
focus  on  exploring  national  patterns  in  potential  treat¬ 
ment  disparities  affecting  women  of  other  racial/ethnic 
groups. 

The  underlying  reasons  for  the  racial  disparities 
observed  in  the  current  study  cohort  remain  to  be  deter¬ 
mined.33  Prior  studies  have  identified  access  to  care  and 
socioeconomic  factors  (such  as  the  patients’  income  and 
education,  the  cost  of  RT,  or  the  availability  of  supple¬ 
mental  insurance)  as  factors  that  may  influence  racial  dis¬ 
parities  in  cancer  care,34’35  but  even  after  adjustment  for 
markers  of  healthcare  access  in  the  current  analysis,  such 
as  number  of  physician  visits  and  use  of  mammography 
before  cancer  diagnosis,  racial  disparities  persisted.  Other 
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unmeasured  potential  explanatory  factors  include  the 
impact  of  the  physician-patient  interaction.36  For  exam¬ 
ple,  it  is  not  known  whether  physicians  offer  treatment 
less  frequently  to  nonwhite  patients,  whether  substandard 
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Figure  1.  Percentage  radiotherapy  use  in  (a)  white  patients 
versus  (b)  black  patients  is  shown.  Gray  shading  indicates 
that  the  sample  size  was  too  small  to  provide  meaningful 
data. 


care  occurs  more  frequently  in  predominantly  nonwhite 
communities,  or  whether  nonwhite  patients  are  more 
likely  to  decline  treatment.  Additional  social  factors,  such 
as  culturally  specific  health  beliefs,  the  presence  of  social 
support,  and  marital  status,37  also  could  affect  racial  dis¬ 
parities  in  care  and  may  be  important  variables  to  explore 
in  future  studies  seeking  to  identify  specific  barriers  in 
breast  cancer  care. 

The  underlying  causes  of  the  geographic  variations 
in  racial  disparities  also  require  further  study.  Several  pub¬ 
lished  analyses  have  indicated  that  geographic  variations 
in  breast  cancer  care  exist,7’8’38  but  to  the  best  of  our 
knowledge  no  prior  study  has  described  how  geographic 
patterns  in  care  may  be  modified  by  patient  race.  Our 
analysis  found  that  disparities  between  white  and  black 
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Figure  2.  Absolute  difference  in  the  rate  of  radiotherapy  (RT) 
use  between  white  and  black  patients  is  shown.  An  absolute 
difference  <0%  indicates  that  the  percentage  of  black 
women  receiving  RT  was  greater  than  the  percentage  of 
white  women.  Gray  shading  indicates  that  the  sample  size 
was  too  small  to  provide  meaningful  data. 


Table  4.  Unadjusted  and  Adjusted  Regional  Variations  in  the  Rate  of  Radiotherapy  Usage  Between  White  and  Black  Patients 


With  Breast  Cancer 

Unadjusted 

Adjusted* 

Region 

States 

%  Whites 

%  Blacks 

P 

%  Whites 

%  Blacks 

P 

West,  Pacific  West 

AK,  CA,  HI,  OR,  WA 

72 

54 

<.001 

74 

61 

.07 

West,  Mountain  West 

AZ,  CO,  ID,  MT,  NV, 

77 

78 

.88 

72 

67 

.59 

NM,  UT,  WY 

Midwest,  West  North  Central 

IA,  KS,  MN,  MO,  NE, 

74 

73 

.94 

74 

70 

.67 

ND,  SD 

Midwest,  East  North  Central 

IL,  IN,  Ml,  OH,  Wl 

76 

72 

.04 

79 

73 

.15 

Northeast,  New  England 

CT,  MA,  NH,  ME,  Rl,  VT 

72 

69 

.53 

78 

71 

<.001 

Northeast,  Mid-Atlantic 

NJ,  NY,  PA 

71 

57 

<.001 

70 

62 

.01 

South,  South  Atlantic 

DE,  DC,  FL,  GA,  MD, 

76 

68 

<.001 

73 

64 

.009 

NC,  SC,  VA,  WV 

South,  West  South  Central 

AR,  LA,  OK,  TX 

74 

63 

.001 

69 

59 

<.001 

South,  East  South  Central 

AL,  KY,  MS,  TN 

72 

57 

<.001 

78 

60 

.001 

•Adjusted  for  all  covariates  except  state/region. 
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patients  in  RT  after  BCS  appeared  most  pronounced  in 
the  northeastern  and  southern  regions  of  the  United 
States.  Given  that  racial  disparities  do  not  appear  to  be 
concentrated  in  a  single  area  of  the  country,  the  underly¬ 
ing  reasons  for  disparities  in  the  different  regions  of  the 
United  States  are  likely  to  be  varied,  complex,  and  multi¬ 
factorial.  For  example,  rural/urban  differences  may  play  a 
role  in  the  geographic  variations  in  treatment  disparities. 
Prior  studies  have  identified  patients’  distance  from  RT 
facilities  as  a  potential  barrier  to  treatment,  ’  and  sug¬ 
gested  that  patients  living  in  rural  areas  may  have  less 
access  to  RT.39,40  However,  a  recent  study  reported  simi¬ 
lar  rates  of  RT  use  after  mastectomy  and  after  lumpec¬ 
tomy  in  urban  versus  rural  locales.41  The  study’s  results 
suggest  that  disadvantaged  patients  in  either  setting  may 
experience  distinct  but  persistent  treatment  barriers.  Pro¬ 
spective  efforts  to  overcome  disparities  in  cancer  treat¬ 
ment,  such  as  the  Patient  Navigation  program,42  are 
currently  being  evaluated  and  may  serve  as  an  important 
source  for  identifying  the  key  barriers  to  cancer  care  in  dif¬ 
ferent  geographic  regions  of  the  United  States. 

The  current  study  has  limitations  to  consider.  Our 
sample  focused  on  older  breast  cancer  patients  with  con¬ 
tinuous  fee-for-service  Medicare  Parts  A  and  B  insurance 
coverage.  Therefore,  future  studies  will  be  required  to  vali¬ 
date  the  magnitude  and  significance  of  racial  disparities  in 
the  treatment  of  younger  patients.  In  addition,  studies  are 
also  needed  to  determine  the  magnitude  of  racial  dispar¬ 
ities  among  patients  with  other  insurance  status,  such  as 
those  with  Medicaid,  given  that  prior  evidence  suggests 
that  patients  with  Medicare  insurance  coverage  are 
actually  more  likely  to  receive  appropriate  breast  cancer 
treatment.43  Furthermore,  our  definitions  of  invasive 
breast  cancer  and  RT  were  claims-based  and,  therefore, 
may  be  subject  to  misclassification  bias.  However,  the 
results  of  the  current  study  demonstrated  reproducibility, 
even  after  applying  the  more  restrictive  criteria  to  exclude 
potential  metastatic  disease  and  DCIS.  Prior  studies  indi¬ 
cate  a  high  degree  of  validity  for  claims-based  treatment 
variables  and  suggest  that,  although  our  treatment  defini¬ 
tion  was  unlikely  to  be  100%  sensitive,  our  percentage  use 
may  underestimate  true  use  rates  by  only  approximately 
2%. 37  Finally,  the  racial  distribution  of  patients  in  the 
current  study  sample  differs  from  the  general  population 
distribution.  Notably,  6%  of  our  sample  was  black, 
whereas  9%  of  US  women  (age  >65  years)  are  black.26 
This  difference  is  likely  due  to  the  tendency  of  nonwhite 
patients  to  receive  mastectomy  over  BCS  even  for  the 
treatment  of  early  stage  disease.44"46 


Conclusions 

To  the  best  of  our  knowledge,  the  current  study  provides 
1  of  the  first  available  data  sets  describing  a  quality  indica¬ 
tor  for  breast  cancer  care  across  a  comprehensive  national 
sample  of  older  breast  cancer  patients  and  was  able  to 
identify  substantial  racial  disparities  in  care  across  the 
United  States.  Our  analysis  helps  define  the  scope  of  the 
treatment  disparities  in  RT  after  BCS  and  underscores 
the  concern  that  this  treatment  disparity  occurs  not 
merely  in  isolation  but  is  instead  a  problem  that  exists  on 
a  national  scale.  Future  efforts  to  improve  breast  cancer 
care  will  require  identifying  and  overcoming  the  underly¬ 
ing  causes  of  these  racial  disparities.  As  additional  data 
become  available,  future  studies  may  also  explore  changes 
in  the  magnitude  of  disparities  over  time  and  the  effect  of 
these  disparities  in  care  on  breast  cancer  outcomes. 
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Abstract 

Background:  In  epidemiologic  studies,  cancer  stage  is  an  important  predictor  of  outcomes.  However,  cancer  stage 
is  typically  unavailable  in  medical  insurance  claims  datasets,  thus  limiting  the  usefulness  of  such  data  for 
epidemiologic  studies.  Therefore,  we  sought  to  develop  an  algorithm  to  predict  cancer  stage  based  on  covariates 
available  from  claims-based  data. 

Methods:  We  identified  a  cohort  of  77,306  women  age  >  66  years  with  stage  l-IV  breast  cancer,  using  the 
Surveillence  Epidemiology  and  End  Results  (SEER)-Medicare  database.  We  formulated  an  algorithm  to  predict 
cancer  stage  using  covariates  (demographic,  tumor,  and  treatment  characteristics)  obtained  from  claims.  Logistic 
regression  models  derived  prediction  equations  in  a  training  set,  and  equations'  test  characteristics  (sensitivity, 
specificity,  positive  predictive  value  (PPV),  and  negative  predictive  value  [NPV])  were  calculated  in  a  validation  set. 
Results:  Of  the  entire  sample  of  women  diagnosed  with  invasive  breast  cancer,  51%  had  stage  I;  26%  stage  II;  11% 
stage  III;  and  4%  stage  IV  disease.  The  equation  predicting  stage  IV  disease  achieved  sensitivity  of  81%,  specificity 
89%,  positive  predictive  value  (PPV)  24%,  and  negative  predictive  value  (NPV)  99%,  while  the  equation 
distinguishing  stage  I/ll  from  stage  III  disease  achieved  sensitivity  83%,  specificity  78%,  PPV  98%,  and  NPV  31%. 
Combined,  the  equations  most  accurately  identified  early  stage  disease  and  ascertained  a  sample  in  which  98%  of 
patients  were  stage  I  or  II. 

Conclusions:  A  claims-based  algorithm  was  utilized  to  predict  breast  cancer  stage,  and  was  particularly  successful 
when  used  to  identify  early  stage  disease.  These  prediction  equations  may  be  applied  in  future  studies  of  breast 
cancer  patients,  substantially  improving  the  utility  of  claims-based  studies  in  this  group.  This  method  may  similarly 
be  employed  to  develop  algorithms  permitting  claims-based  epidemiologic  studies  of  patients  with  other  cancers. 


Background 

Administrative  medical  insurance  claims  are  an  impor¬ 
tant  source  of  population-based  data  used  in  epidemio¬ 
logic  studies  of  various  diseases.  Specifically,  in  older 
patients,  national  Medicare  data  have  been  useful  for 
the  study  of  many  conditions,  including  myocardial 
infarction,  heart  failure,  chronic  kidney  disease,  Parkin¬ 
son’s  disease,  and  venous  thromboembolism  [1-5].  How¬ 
ever,  for  studying  cancer,  the  use  of  national  Medicare 
data  has,  to  date,  been  limited.  Medicare  claims  data  are 
clearly  recognized  as  potentially  a  rich  source  for  cancer 
epidemiology  and  outcomes  research,  and  in  fact 
demonstrate  acceptable  validity  for  identifying  cancer 
diagnoses  and  treatment  patterns  [6-13].  Unfortunately, 

*  Correspondence:  tbuchhol@mdanderson.org 

department  of  Radiation  Oncology,  The  University  of  Texas  MD  Anderson 
Cancer  Center,  1515  Holcombe  Blvd  Houston,  Texas  77030,  USA 

(44  BioMed  Central 


the  lack  of  cancer  stage  data  in  Medicare  claims  remains 
a  major  limiting  factor  in  maximizing  the  utility  of  these 
datasets  for  retrospective,  outcomes-based  research  in 
cancer  patients[14,15].  In  particular,  cancer  stage  is  a 
crucial  predictor  of  disease  outcome  and  a  key  factor  in 
determining  the  appropriateness  of  treatment.  For 
example,  in  breast  cancer  patients,  stage  is  associated 
with  overall  and  disease-free  survival  and,  furthermore, 
stage  influences  treatment  decisions  such  as  selection 
and  timing  of  surgery,  radiotherapy,  and  chemotherapy 
[16].  Epidemiologic  studies  of  cancer  patients  typically 
employ  stage  variables  as  covariates  or  as  inclusion  and 
exclusion  criteria,  and  thus  it  is  essential  to  develop 
accurate  algorithms  to  account  for  cancer  stage  in  stu¬ 
dies  using  claims  data. 

Surprisingly,  the  need  to  derive  such  algorithms  has 
largely  been  ignored  in  the  literature.  Only  one  prior 
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study  has  developed  a  claims-based  algorithm  to  predict 
stage  in  breast  cancer  patients.  Cooper  et  al.  used  the 
Surveillance  Epidemiology  and  End  Results  (SEER)-Med- 
icare  database,  and  authors  reported  that  their  claims- 
based,  single-predictor  models  were  insufficient  for 
identifying  patients  with  local,  regional,  and  distant 
stage  disease.  Sensitivity  of  these  models  for  distinguish¬ 
ing  local  from  regional  and  distant  disease  was  low-for 
example,  in  breast  cancer  patient  samples,  only  approxi¬ 
mately  60%[17].  In  the  decade  since  the  prior  algorithm 
was  derived,  no  other  algorithm  has  been  presented  in 
the  literature  attempting  to  improve  cancer  stage  classi¬ 
fication  using  claims  data.  Thus  studies  have  continued 
to  apply  the  algorithm  by  Cooper  and  colleagues  to 
derive  cancer  stage  variables,  despite  the  recognized  lim¬ 
itations  of  this  algorithm  [18,19]  and  the  introduction  of 
measurement  errors  into  such  analyses. 

Accordingly,  we  sought  to  derive  an  expanded  predic¬ 
tive  algorithm  based  on  multivariate  modeling  and  to 
improve  the  sensitivity  and  specificity  for  identifying 
cancer  stage,  using  our  study  sample  of  breast  cancer 
patients  as  an  illustrative  case.  Using  available  Medicare 
claims  for  breast  cancer  patients  found  in  the  SEER- 
Medicare  database,  we  developed  a  prediction  algorithm 
to  identify  patients  with  distant  (stage  IV)  disease  at 
diagnosis  and,  among  patients  without  distant  disease,  a 
prediction  algorithm  to  classify  the  extent  of  locoregio- 
nal  (stages  I-III)  disease. 

Methods 

Algorithm 
Study  sample 

The  SEER-Medicare  database  is  comprised  of  a  popula¬ 
tion-based  cohort  of  Medicare  beneficiaries  with  incident 
cancer  identified  through  SEER  registries,  which  account 
for  up  to  26%  of  the  United  States’  population  [20,21]. 
Our  initial  study  population  consisted  of  150,764  women 
(age  >  65  years)  diagnosed  with  breast  cancer  between 
1992  and  2002  identified  through  SEER-Medicare.  From 
this  population,  we  excluded  5,217  patients  with 
unknown  SEER  historic  stage  (as  this  variable  indicated 
the  presence  or  absence  of  metastases),  and  19,816  with 
in  situ  disease  (as  we  intended  to  focus  only  on  invasive 
disease).  We  further  excluded  47,114  patients  who  did 
not  have  continuous  Medicare  Fee-for-Service  coverage 
or  had  any  HMO  coverage  from  12  months  prior  to  9 
months  after  their  diagnosis  date  (as  claims  information 
might  be  incomplete  during  these  periods),  and  the  308 
patients  age  <66,  since  these  patients  potentially  would 
not  have  had  comprehensive  claims  information  to  define 
the  independent  predictor  covariates.  We  finally  excluded 
1,003  patients  who  died  or  were  lost  to  follow-up  within 
9  months  of  their  diagnosis  date.  This  yielded  a  final 
sample  size  of  77,306  patients  in  our  study. 


Dependent  variable:  Cancer  stage 

The  “gold  standard”  for  identifying  cancer  stage  at  diag¬ 
nosis  was  determined  using  a  combination  of  tumor  vari¬ 
ables  available  through  SEER.  Distant  disease  was 
determined  using  the  American  Joint  Committee  on 
Cancer  (AJCC)  [22]  historic  stage  as  reported  to  SEER, 
which  indicates  tumor  present  in  any  distant  site  at  can¬ 
cer  diagnosis  (compared  with  tumor  limited  only  to  local 
or  regional  sites  at  diagnosis).  For  our  analysis,  patients 
with  any  distant  disease  were  considered  stage  IV. 

Patients  without  distant  (stage  IV)  disease  had  local  or 
regional  AJCC  historic  stage.  T  and  N  classification  in 
these  patients  were  assigned  based  on  SEER  variables 
for  tumor  size  and  extent  of  disease.  Tumor  size  and 
extent  were  categorized  as  <  2  cm  (Tl);  >2  to  5  cm 
(T2);  >5  cm  (T3);  or  invading  into  the  chest  wall,  ribs, 
intercostals  or  serratus  anterior  muscles,  extensive  inva¬ 
sion  into  the  skin,  inflammatory  carcinoma,  or  further 
contiguous  extension  into  the  skin  (T4).  Nodal  disease 
was  categorized  as  0  positive  lymph  nodes  (NO);  1-3 
positive  lymph  nodes  (Nl);  4-9  positive  lymph  nodes 
(N2);  or  10  or  more  positive  lymph  nodes  (N3)[21].  Due 
to  the  extent  of  missing  data  in  the  SEER  database,  loca¬ 
tion  of  positive  lymph  nodes  was  not  included  in  N  clas¬ 
sification.  Stage  I  included  T1N0  disease,  stage  II 
included  TONI,  T1N1,  T2N0,  T2N1,  and  T3N0  disease, 
and  stage  III  included  T0-2N2,  T3N1-2,  T4N0-2,  and 
T0-4N3  disease.  These  classifications  are  based  on 
AJCC  2003  staging  criteria  [22]. 

Independent  predictors 

Candidate  independent  predictors  were  selected  a  priori 
based  on  statistical  significance  in  bivariate  analyses  (P  < 
0.25)  and  clinical  significance  in  prior  studies  of  cancer 
patients[20, 23-28].  Variables  were  defined  by  searching 
through  inpatient,  outpatient,  and  carrier  Medicare 
claims  or  the  denominator  file  for  SEER-Medicare 
linked  data  for  demographic  variables.  A  comprehensive 
list  of  International  Classification  of  Diseases,  Ninth 
Revision  (ICD-9),  Common  Procedural  Terminology 
(CPT),  and  Revenue  Center  codes  for  each  predictor  are 
listed  in  Table  SI,  Additional  file  1. 

Statistical  analysis 

All  statistical  analyses  were  conducted  using  SAS  ver¬ 
sion  9.1.3  (SAS  Institute  Inc,  Cary,  NC),  and  all  statisti¬ 
cal  tests  assumed  a  2-tailed  a  of  0.05.  The  University  of 
Texas  M.  D.  Anderson  Cancer  Center  institutional 
review  board  deemed  this  study  exempt  from  review, 
since  the  data  were  without  identifiers. 

We  derived  two  separate  logistic  models  and  imple¬ 
mented  the  models  sequentially.  The  first  model  tested 
the  associations  between  predictor  covariates  and  the 
dichotomous  outcome  of  stage  IV  versus  non-stage  IV 
disease.  Among  the  subset  of  patients  who  were  not 
categorized  as  having  stage  IV  disease,  the  second 
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model  tested  the  associations  between  predictor  covari¬ 
ates  (excluding  the  covariate  for  metastatic  disease  at 
diagnosis)  and  the  dichotomous  outcome  of  stage  I/II 
(early)  versus  stage  III  disease.  Outcomes  were  dichoto¬ 
mized  based  on  clinical  rationale,  given  that  treatment 
of  metastatic  disease  is  palliative;  and  that  treatment  of 
early  stage  disease  is  distinct  in  that  breast  conserving 
therapy  is  a  treatment  option. 

We  used  a  split  sample  approach  to  develop  and  vali¬ 
date  our  logistic  models.  Each  model  was  derived  from 
the  “training  set,”  selected  using  simple  random  sam¬ 
pling  without  replacement  (38,653  of  77,306  patients). 
Parsimonious  models  were  then  selected  based  on  statis¬ 
tical  significance  ( P  <  0.25),  clinical  significance  of  cov¬ 
ariates  in  prior  studies, [20,23-28]  and  goodness-of-fit. 
Prior  studies  were  used  as  an  initial  guide  for  the  selec¬ 
tion  of  covariates  to  consider.  The  significance  cutoff  (P 
<  0.25)  was  used  to  rule  in  covariates  to  keep.  Examin¬ 
ing  the  goodness-of-fit  of  the  overall  model  was  used  to 
rule  out  covariates  to  exclude.  In  combination,  these 
three  criteria  were  used  to  select  the  final  model. 

Testing 

Patients  not  included  in  the  training  set  constituted  the 
“validation  set”.  In  the  validation  set,  the  parameter  asso¬ 
ciated  with  each  covariate  estimated  from  the  derivation 


set  was  applied  to  each  patient  in  the  validation  set  to  cal¬ 
culate  each  patient’s  predictive  probability  (calculated 
probability  =  ^0+^Xl''+^^  +  ''-)/+exp(j3o+Jgi^+J82;Caj  +  ...)  ) 

of  having  stage  IV  disease  in  the  first  model  and  stage  I/II 
or  stage  III  disease  in  the  second  model.  Test  characteris¬ 
tics  were  calculated  for  probability  cutpoints  between  0.05 
and  0.90,  using  two-by-two  tables.  The  “gold  standard”  for 
stage  was  considered  the  SEER  stage;  the  test  stage  was 
based  on  the  calculated  probability  (for  example,  for  a 
probability  cutpoint  of  0.05,  patients  were  predicted  to 
have  stage  I/II  disease  if  their  calculated  probability  was  > 
0.05,  and  not  to  have  stage  I/II  disease  if  their  calculated 
probability  was  <0.05). 

Combining  equations 

The  prediction  equations  were  then  applied  to  isolate  a 
sample  of  patients  with  early  stage  disease.  Specifically, 
the  two  prediction  equations  were  applied  sequentially 
to  the  validation  sample  in  order  to  identify  a  subset  of 
patients  with  stage  I/II  disease.  The  first  step  used  a 
probability  cutpoint  of  0.05  to  exclude  patients  with  pre¬ 
dicted  stage  IV  disease.  The  second  step  was  applied  to 
the  subset  identified  in  the  first  step  and  used  a  prob¬ 
ability  cutpoint  of  0.90  to  include  patients  with  pre¬ 
dicted  stage  I/II  disease.  These  cutpoints  were  chosen 
based  on  their  test  characteristics  (high  sensitivity,  speci¬ 
ficity,  and  positive  predictive  value  [PPV]  or  negative 
predictive  value  [NPV]).  Finally,  we  also  compared  the 
test  characteristics  derived  from  our  multivariate  predic¬ 
tion  equations  to  test  characteristics  derived  from 


single-predictor  equations  for  distant  and  regional  dis¬ 
ease  used  in  a  prior  study, [17]  to  determine  whether 
multivariate  equations  improved  test  characteristics 
compared  with  the  single-predictor  equations. 

Implementation 

Example  in  Practice:  Medicare  test  sample 

Finally,  we  present  an  example  that  applies  the  predic¬ 
tion  equation.  We  used  a  test  sample  based  on  a  claims- 
only  dataset,  the  national  Medicare  dataset.  The  national 
Medicare  dataset  includes  claims  data  for  all  Medicare 
beneficiaries  in  the  United  States.  Files  contain  data  col¬ 
lected  by  Medicare  for  reimbursement  of  health  care 
services  for  each  beneficiary  and  include  institutional 
(inpatient  and  outpatient)  and  non-institutional  (physi¬ 
cians  or  other  providers)  final  action  claims  [29].  We 
initially  included  127,607  women  (age  >  65)  with  a  diag¬ 
nosis  claim  indicating  invasive  breast  cancer  in  2003 
(International  Classification  of  Diseases,  Ninth  Revision 
(ICD-9)  diagnosis  code  of  174)  who  underwent  a  breast- 
cancer  related  procedure.  We  then  excluded  23,715 
patients  who  did  not  have  at  least  2  claims  (on  different 
dates)  specifying  a  diagnosis  of  invasive  breast  cancer 
between  January  1,  2003  and  December  31,  2004  (at 
least  1  claim  must  have  occurred  during  2003);  16,471 
patients  who  had  a  breast  cancer-related  diagnosis  or 
procedure  claim  between  January  1,  2002,  and  Decem¬ 
ber  31,  2002;  5,719  patients  who  were  receiving  Medi¬ 
care  coverage  because  of  end-stage  renal  disease  or 
disability;  and  6,612  patients  who  lacked  Part  A  or  B 
coverage  or  who  had  intermittent  health  maintenance 
organization  coverage  in  the  9  months  after  or  in  the  1 
year  before  their  breast  cancer  diagnosis  date  (of  these 
patients,  3,561  had  incomplete  information  in  the  year 
prior  to  diagnosis  because  they  were  <66  years  of  age); 
for  a  total  sample  size  of  56,725  patients.  This  method 
for  sample  selection  of  incident  breast  cancer  has  been 
validated  in  a  prior  study  [30]. 

In  this  test  sample,  we  applied  our  derived  algorithm 
and  calculated  the  frequency  of  patients  classified  as 
predicted  stage  IV  and  predicted  stage  I/II  disease. 
Again,  for  this  sample,  the  first  step  used  a  probability 
cutpoint  of  0.05  and  the  second  step  a  probability  cut- 
point  of  0.90.  As  a  test  of  our  algorithm  for  consistency, 
we  compared  the  predicted  frequencies  to  the  actual 
stage  distribution  in  two  populations:  1)  the  SEER-Medi- 
care  population  (age  >65  years)  and  2)  the  National 
Cancer  Data  Base  population  (age  >  70  years)  [31]. 

Results 

Patient  characteristics 

In  77,306  women,  mean  age  was  76  years  (standard 
deviation  7  years),  and  94%  were  white.  Fifty-one  per¬ 
cent  were  stage  I  (39,147),  26%  stage  II  (19,967),  11% 
stage  III  (8,174),  4%  stage  IV  (3,220)  and  9%  with  non- 
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distant  disease  but  T  or  N  classification  unknown 
(6,798).  Forty-five  percent  were  treated  with  breast  con¬ 
serving  surgery,  49%  with  mastectomy,  44%  with  radio¬ 
therapy,  and  18%  with  chemotherapy  (Table  1). 

Prediction  algorithm  equations  and  test  characteristics  for 
probability  cutpoints 

Candidate  covariates  tested  in  prediction  equations  are 
listed  in  Table  SI,  Additional  file  1.  Parameter  estimates 
for  the  covariates  included  in  each  final  prediction  equa¬ 
tion  are  listed  in  Table  S2,  Additional  file  1. 

Stage  IV 

Fourteen  percent  of  all  patients  and  73%  of  patients  with 
stage  IV  disease  had  a  claims  code  indicating  possible 


Table  1  Study  Sample  Patient  Characteristics,  N  =  77,306 


Predictor  Variable 

%  of  All  Patients 

Demographic 

Age,  mean  (SD) 

76  (7) 

White  race 

94 

Stage 

Stage  1 

51 

Stage  II 

26 

Stage  III 

10 

Stage  IV 

4 

Stages  I-III  but  T  or  N  unknown 

9 

Extent  of  disease3 

Axillary  LN  involvement 

19 

Metastatic  disease 

14 

Cancer  treatment3 

No.  visits  to  surgeon,  mean  (SD) 

4(3) 

No.  visits  to  medical  oncologist,  mean  (SD) 

4(9) 

No.  visits  to  radiation  oncologist,  mean  (SD) 

3  (5) 

Imaging  (CT,  MRI,  PET,  or  bone  scan) 

25 

Radiation  therapy 

44 

Breast  conserving  surgery 

45 

Mastectomy 

49 

Axillary  LN  dissection 

72 

Chemotherapy  (any  agent) 

18 

Doxorubicin 

7 

Paclitaxel 

3 

No.  physician  visits,  mean  (SD) 

14(12) 

Screening  mammography 

78 

Influenza  vaccine 

34 

General  health  status3 

No.  hospital  admission  for  any  cause,  mean  (SD) 

1  (D 

Charlson  comorbidity  score 

0  comorbid  conditions 

69 

1  comorbid  condition 

18 

2  or  more  comorbid  conditions 

8 

Unknown 

5 

Abbreviations:  CT  computed  tomography;  LN  lymph  nodes;  MRI  magnetic 
resonance  imaging  PET  positron  emission  tomography;  SD  standard  deviation 
a  As  indicated  by  Medicare  claims  codes;  thus  percentage  of  patients  with 
code  for  metastatic  disease  not  equal  to  patients  with  Stage  IV  disease. 


metastatic  disease.  Accordingly,  the  single-predictor 
model  including  only  this  covariate  had  sensitivity  of 
73%;  specificity  89%;  PPV  22%;  and  NPV  99%  for  identi¬ 
fying  stage  IV  disease.  After  including  covariates  (Table 
S2,  Additional  file  1)  in  the  multivariate  model,  the  sen¬ 
sitivity  was  81%  (95%  Confidence  Interval  80%  -  84%) 
for  identifying  stage  IV  disease  at  a  probability  cutpoint 
of  0.05.  At  this  cutpoint,  specificity  was  89%  (86%  - 
89%);  PPV  24%;  (22%  -  25%)  and  NPV  99%  (99%  -  99%), 
yielding  a  c-statistic  of  0.93.  (Table  S3,  Additional  file  1) 
The  distribution  of  calculated  predicted  probabilities  in 
the  validation  set  for  patients  with  stage  IV  versus  stage 
I-III  disease  is  presented  in  Figure  la. 

Stages  l-lll 

In  patients  with  stages  I-III  disease,  19%  had  a  claims 
code  indicating  axillary  lymph  node  involvement.  Speci¬ 
fically,  2%  of  patients  with  stage  I  disease,  36%  with 
stage  II  disease,  and  64%  of  patients  with  stage  III  dis¬ 
ease  had  this  claims  code.  The  single-predictor  model 
including  only  this  covariate  yielded  sensitivity  of  87% 
(specificity  61%;  PPV  94%;  NPV  39%)  for  identifying 
stage  I/II  disease  and  sensitivity  of  61%  (and  specificity 
87%;  PPV  39%;  NPV  94%)  for  identifying  stage  III  dis¬ 
ease.  After  including  covariates  (Table  S2,  Additional 
file  1)  in  the  multivariate  model,  the  sensitivity  was  91% 
(90%  -  92%)  for  identifying  stage  I/II  disease  at  a  prob¬ 
ability  cutpoint  of  0.80;  and  83%  (83%  -  85%)  at  a  cut- 
point  of  0.90.  At  a  cutpoint  of  0.90,  specificity  was  78% 
(75%  -  79%);  PPV  98%  (97%  -  98%);  and  NPV  31%  (30% 
-  34%).  For  identifying  stage  III  disease,  the  sensitivity 
was  78%  (75%  -  79%)  at  a  cutpoint  of  0.10.  At  this  cut- 
point,  specificity  was  83%  (83%  -  85%);  PPV  30%  (30%  - 
34%);  and  NPV  98%  (97%  -  98%)  (Table  S3,  Additional 
file  1).  These  models  yielded  a  c-statistic  of  0.88.  The 
distribution  of  calculated  predicted  probabilities  in  the 
validation  set  for  patients  with  stage  I/II  versus  stage  III 
disease  is  presented  in  Figure  lb. 

Comparison  with  other  single  predictors 
For  comparison’s  sake,  for  identifying  stage  IV  disease, 
the  second  most  important  predictor  was  axillary  lymph 
node  dissection.  This  predictor  alone  would  yield  the 
following  test  characteristics:  sensitivity  67%;  specificity 
74%;  PPV  10%;  and  NPV  98%.  For  identifying  stage  I/II 
disease,  the  second  most  important  predictor  was  breast 
conserving  surgery  vs.  mastectomy,  yielding  the  follow¬ 
ing  test  characteristics:  sensitivity  49%;  specificity  82%; 
PPV  95%;  and  NPV  18%. 

Combining  equations 

The  prediction  equations  were  most  accurate  for  isolat¬ 
ing  patients  with  early  stage  disease.  Specifically,  after 
applying  the  two  prediction  equations  sequentially  to 
the  validation  sample  to  identify  patients  with  predicted 
stage  I/II  disease,  a  subset  of  23,285  patients  were 
selected  (of  38,653  patients,  36,417  were  predicted  to 
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Figure  1  Distribution  of  calculated  probabilities  (%)  for  patients  with  predicted  stages  l-lll  and  stage  IV  disease  (a)  Distribution  of 
calculated  probabilities  for  patients  with  predicted  stages  l-ll  and  stage  III  disease  (b).  A  calculated  probability  of  5%  corresponds  to  a  cutpoint 
of  0.05.  The  histograms  are  normalized  to  100%. 
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have  non-stage  IV  disease,  and  of  these  patients,  23,285 
were  predicted  to  have  stage  I/II  disease).  The  predictive 
sample  actually  consisted  of  98%  gold  standard  stage  I/II 
disease  (22,706  of  23,285),  2%  stage  III  disease  (549  of 
23,285),  and  <1%  stage  IV  disease  (110  of  23,285).  Of  all 
patients  with  gold  standard  stage  I/II  disease  (29,546  of 
38,653  validation  patients),  23%  (6,840  of  29,546)  were 
excluded  (classified  as  other  than  stage  I/II)  as  a  result 
of  the  algorithm  (4,604  from  the  first  model  and  2,236 
from  the  second  model).  (Figure  2,  Figure  3). 

Example  in  Practice:  Medicare  test  sample  and 
comparison  for  consistency 

In  our  Medicare  test  sample,  after  the  first  predictor 
equation  was  applied,  a  total  of  4%  (2,333  of  56,725)  of 
women  were  predicted  as  having  stage  IV  disease.  This 
compared  favorably  with  the  SEER-Medicare  population, 
which  included  4%  (3,220  of  77,306)  of  women  with 
confirmed  stage  IV  disease;  as  well  as  the  NCDB  popu¬ 
lation,  which  included  5%  (1,913  of  41,071)  confirmed 
stage  IV  disease.  After  the  second  predictor  equation 
was  applied  to  the  remainder  of  the  test  sample,  a  total 
of  79%  (43,169  of  54,392)  of  women  were  predicted  as 
having  stage  I/II  disease.  This  compared  favorably  with 
the  SEER-Medicare  population,  which  included  80% 
(59,114  of  74,086)  of  women  with  confirmed  Stage  I/II 
disease;  as  well  as  the  NCDB  population,  which  included 
84%  (33,036  of  39,158)  confirmed  Stage  I/II  disease. 

Discussion 

In  this  cohort  of  older  breast  cancer  patients,  Medicare 
claims  data  assisted  the  prediction  of  cancer  stage. 


Predictor  equations  using  claims  data  alone  were  able  to 
achieve  approximately  80%  sensitivity  and  specificity  for 
identifying  stage  IV  disease  as  well  as  distinguishing 
stage  I/II  from  stage  III  disease.  Prediction  models  maxi¬ 
mized  NPV  when  distinguishing  stage  IV  from  stage  I- 
III  disease  but  maximized  PPV  when  distinguishing 
stage  I/II  from  III  disease.  With  a  resulting  tradeoff  of 
lower  PPV  in  the  first  model  and  lower  NPV  in  the  sec¬ 
ond  model,  the  algorithm  was  therefore  found  to  be 
best  suited  to  most  accurately  identify  early  stage  dis¬ 
ease.  Specifically,  an  algorithm  combining  the  two  equa¬ 
tions  seeking  to  identify  patients  with  early  stage  disease 
was  able  to  achieve  a  sample  in  which  98%  of  patients 
had  stage  I  or  II  disease. 

Our  prediction  models  represent  an  improvement  over 
the  single  other  previously  published  model.  In  this 
prior  study,  Cooper  et  al.  used  single-predictor  equa¬ 
tions  to  identify  cancer  stage.  To  identify  patients  with 
distant  disease,  authors  tested  a  single  variable  based  on 
claims  codes  for  metastatic  disease.  This  single-predictor 
model  demonstrated  60%  sensitivity  and  58%  PPV.  To 
distinguish  patients  with  local  versus  regional  disease, 
authors  tested  a  single  variable  based  on  the  claim  code 
for  axillary  lymph  node  involvement.  This  single-predic¬ 
tor  model  demonstrated  62%  sensitivity  and  85%  PPV 
[17].  The  relatively  poor  test  characteristics  from  this 
prior  study  demonstrated  that  these  single-predictor 
models  would  be  insufficient  for  predicting  stage  in 
patients  with  breast  cancer  and  suggested  that  claims 
data  alone  would  be  inadequate  for  epidemiologic  stu¬ 
dies  of  cancer  patients. 


1-Specificity 

Figure  2  Receiver  Operating  Curve  (ROC)  for  equation  to  predict  stage  IV  disease 
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1-Specificity 

Figure  3  Receiver  Operating  Curve  (ROC)  for  equation  to  predict  stage  l-lll  disease 


In  contrast,  our  prediction  models  have  improved 
upon  these  test  characteristics.  Our  single  predictor 
model  demonstrated  improvement,  likely  in  part  due  to 
a  more  extensive  list  of  claims  codes,  with  multiple  cov¬ 
ariates  providing  added  value.  Moreover,  our  algorithm 
demonstrated  consistency  when  results  were  compared 
with  population-based  data  from  SEER-Medicare  and 
the  NCDB.  There  are  two  important  future  research 
applications  of  our  prediction  models.  First,  our  multi¬ 
variate  logistic  modeling  method  for  developing  a  stage 
predictor  algorithm  may  similarly  be  applied  to  test 
models  and  potentially  develop  stage  prediction  equa¬ 
tions  for  patients  diagnosed  with  cancers  of  other  sites. 
Second,  our  prediction  equations  may  also  be  applied 
directly  to  claims-based  databases  of  breast  cancer 
patients  who  have  unknown  stage.  Using  a  combination 
of  multiple  predictors  along  with  claims  codes  for  meta¬ 
static  disease  and  axillary  lymph  node  dissection,  para¬ 
meter  estimates  and  calculated  probabilities  can  be 
applied  to  the  prediction  of  patient  breast  cancer  stage. 

Our  algorithm  can  therefore  serve  as  a  tool  to  assist  in 
the  investigation  of  a  variety  of  epidemiologic  research 
questions  in  breast  cancer  patients  by  allowing  a  sample 
selection  of  those  patients  with  early  stage  disease.  In 
addition,  predicted  early  stage  disease  can  be  applied  as 
a  covariate.  Accordingly,  since  disease  stage  may  be  bet¬ 
ter  accounted  for,  claims  databases  of  breast  cancer 
patients  may  also  be  better  applied  to  address  such 
questions  as  the  assessment  of  treatment  utilization, 
geographic  variation,  or  outcomes  in  patients  diagnosed 
with  early  stage  breast  cancer.  Specifically,  for  the  stage 
IV  prediction  model,  a  probability  cutpoint  between  > 
0.05  and  >  0.10  would  be  highly  specific  and  sensitive 


for  identifying  patients  with  stage  IV  disease.  For  the 
stage  I/II  prediction  model,  a  cutpoint  between  >  0.80 
and  >  0.90  would  be  highly  specific  and  sensitive  for  dis¬ 
tinguishing  patients  with  stage  I/II  disease  from  patients 
with  stage  III  disease. 

For  the  identification  of  patients  with  stage  IV  disease, 
a  selected  probability  cutpoint  criterion  could  be  trans¬ 
lated  into  a  dichotomous  variable,  and  used  either  to 
select  a  sample  of  patients  with  stage  IV  disease  or  used 
in  a  “rule  out”  context,  as  an  exclusion  criterion.  The 
high  NPV  in  our  proposed  model  suggests  that  when 
using  these  outpoints  to  identify  a  sample  limited  to 
patients  with  stage  I-III  disease,  the  likelihood  of  mis- 
classification  bias  (bias  due  to  the  inappropriate  inclu¬ 
sion  of  patients  with  stage  IV  disease  in  the  sample) 
would  be  low  in  a  “rule-out”  setting. 

For  distinguishing  patients  with  stage  I/II  versus  III  dis¬ 
ease,  the  probability  cutpoint  criterion,  translated  into  a 
dichotomous  variable,  could  be  useful  in  various  contexts, 
such  as  excluding  patients  with  stage  III  disease  in  order 
to  refine  a  study  population  of  patients  with  early  stage 
breast  cancer,  or  creating  a  dichotomous  covariate  to 
adjust  for  potential  confounding  associated  with  stage  I-III 
disease.  The  test  characteristics  in  our  analysis  suggest 
that  the  combination  of  these  prediction  equations  may  be 
particularly  useful  in  the  context  of  identifying  breast  can¬ 
cer  patients  with  early  stage  (stage  I  and  II)  disease. 

Our  study  has  limitations  to  consider.  First,  our  cohort 
was  limited  only  to  older  patients  with  breast  cancer. 
Although  the  variables  associated  with  stage  are  likely  to 
be  similar  in  younger  patients,  exact  parameter  estimates 
may  differ,  and  the  application  of  these  models  in 
younger  patients  requires  further  validation.  Additionally, 
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as  our  predictor  variables  were  derived  from  Medicare 
claims,  these  models  will  also  require  validation  in  other 
claims  based  data.  If  not  all  the  proposed  variables  in  our 
models  are  available,  however,  at  a  minimum,  adjustment 
in  multivariate  analysis  for  as  many  possible  available 
candidate  predictors  proposed  in  our  study  could  be  use¬ 
ful  to  improve  modeling  of  breast  cancer  outcomes  in 
future  studies.  Although  we  excluded  from  our  parsimo¬ 
nious  model  covariates  that  required  long-term  follow-up 
(specifically,  overall  survival  and  mastectomy  9  or  more 
months  after  diagnosis),  our  models  still  required  both 
retrospective  and  prospective  data  for  up  to  1  year  prior 
to  and  1  year  after  the  date  of  diagnosis.  Thus,  studies 
applying  our  models  would  be  limited  to  patients  with 
continuous  coverage  and  complete  claims  information 
over  this  time  period.  The  gold  standard  for  our  out¬ 
come,  cancer  stage,  was  based  on  pathologic  variables  in 
SEER-Medicare,  though  given  a  lack  of  central  pathology 
review  by  the  SEER  program,  unmeasured  error  may 
have  affected  the  gold  standard,  yielding  potentially  less 
than  100%  accuracy.  Finally,  if  a  sample  was  selected 
based  on  the  algorithm,  sample  characteristics  derived 
from  algorithm  predictor  variables  (for  example,  che¬ 
motherapy,  radiotherapy,  and  surgery  utilization)  may  be 
under-  or  overestimated  compared  with  the  complete 
patient  population,  depending  on  the  direction  and  sig¬ 
nificance  of  their  association  with  disease  stage  in  the 
prediction  equations. 

Conclusions 

Medicare  claims  data  can  be  utilized  to  derive  a  useful 
algorithm  to  predict  stage  in  breast  cancer  patients.  In 
particular,  the  predicted  probability  of  early  stage  disease 
can  be  easily  generated  when  applying  the  prediction 
algorithm  to  this  patient  population,  thus  substantially 
improving  the  utility  of  Medicare  claims  data  for  study¬ 
ing  breast  cancer. 


Additional  file  1:  Table  SI.  Candidate  Covariates  and  Claims  Codes. 
Table  S2.  Prediction  Equations.  Parameter  estimates  for  stage  IV  versus 
stages  l-lll  disease  and  for  stage.  I/ll  versus  stage  III  disease.  Table  S3.  Test 
Characteristics  After  Applying  Prediction  Equations  on  Validation  Set 
Samples. 
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BACKGROUND: 


STUDY  DESIGN: 


RESULTS: 


CONCLUSIONS: 


Breast-conserving  surgery  (BCS)  followed  by  radiotherapy  is  as  effective  as  mastectomy  for 
treatment  of  early  invasive  breast  cancer.  But  earlier  studies  report  low  BCS  use  rates  of  12%  to 
43%  nationally,  especially  in  older  patients.  We  sought  to  determine  current  patterns  and 
predictors  of  BCS  use. 

In  a  national  Medicare  database  of  all  beneficiaries  (age  greater  than  65  years)  with  incident  invasive 
breast  cancer  treated  with  operation  in  2003,  claims  codes  identified  BCS  versus  mastectomy  and 
demographic,  treatment,  and  geographic  region  covariates.  The  2003  Area  Resource  File  provided 
socioeconomic  covariates.  Logistic  regression  modeled  predictors  of  BCS. 

In  56,725  women,  59%  were  treated  with  BCS  versus  41%  with  mastectomy.  BCS  was  more 
likely  in  women  who  were  younger  than  70  years  (odds  ratio  [OR],  1.37;  95%  Cl,  1.31  to  1.44; 
p  <  0.001)  and  had  lymph  node-negative  disease  (OR,  1.60;  95%  Cl,  1.52  to  1.68;  p  <  0.001). 
Socioeconomic  factors  influenced  use,  with  BCS  more  likely  in  areas  with  low  poverty  (OR, 
1.05;  95%  Cl,  1.00  to  1.09;  p  =  0.03),  high  education  (OR,  1.13;  95%  Cl,  1.08  to  1.19),  high 
density  of  radiation  oncologists  (OR,  1.30;  95%  Cl,  1.06  to  1.59),  and  in  metropolitan  areas 
(OR,  1.20;  95%  Cl,  1.14  to  1.26).  Significant  geographic  variation  existed:  70%  of  women 
were  treated  with  BCS  in  northeastern  New  England  compared  with  only  48%  to  50%  in  the 
South  (p  <  0.001). 

Currently,  more  than  half  of  older  women  across  the  US  diagnosed  with  nonmetastatic  invasive 
breast  cancer  treated  surgically  receive  BCS,  representing  a  substantial  increased  use  compared 
with  historical  data.  Lack  of  BCS  use  appears  in  part  associated  with  socioeconomic  disadvan¬ 
tage,  suggesting  that  persistent  barriers  to  breast  conservation  exist.  (J  Am  Coll  Surg  2009;209: 
425-433.  ©  2009  by  the  American  College  of  Surgeons) 


In  patients  diagnosed  with  invasive  breast  cancer,  breast- 
conserving  surgery  (BCS)  followed  by  radiotherapy  (RT)  is 
considered  as  effective  as  mastectomy  for  treatment  of 
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early-stage  disease.1  Cumulative  data  from  randomized  tri¬ 
als  demonstrate  that  patients  treated  with  breast- 
conserving  therapy  have  overall  survival  rates  comparable 
to  those  in  patients  treated  with  mastectomy.2,3  A  20-year 
update  of  the  National  Surgical  Adjuvant  Breast  and  Bowel 
Project  (NSABP)-B06  trial  affirmed  that  these  treatment 
options  were  still  comparable  even  on  longterm  followup.2 
Given  the  strong  evidence  supporting  the  effectiveness 
of  breast-conserving  therapy  for  treating  early-stage  dis¬ 
ease,  many  experts  have  recommended  that  this  option 
be  selected,  when  possible,  in  order  to  promote  breast 
conservation.3'5 

Despite  the  recommendations  supporting  breast-conserving 
therapy  studies  of  breast  cancer  patients  treated  in  the 
1980s  and  1990s  reported  low  rates  of  BCS,  with  less  than 
half  of  all  surgically  treated  patients  with  nonmetastatic 
invasive  disease  receiving  BCS.6,7  These  studies  noted  sig- 
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Abbreviations  and  Acronyms 

BCS  =  breast-conserving  surgery 

NSABP  =  National  Surgical  Adjuvant  Breast  and  Bowel 
Project 

OR  =  odds  ratio 

RT  =  radiotherapy 

SEER  =  Surveillance,  Epidemiology,  and  End  Results 


nificant  geographic  variation  in  treatment,  but  even  in  re¬ 
gions  with  the  highest  use,  only  about  40%  of  patients  were 
treated  with  BCS.6  Results  further  suggested  that  during 
these  earlier  eras,  significant  treatment  disparities  existed, 
with  non-Caucasian  patients  and  patients  facing  health¬ 
care  access  barriers  less  likely  to  receive  BCS.5'9 

Contemporary  data  are  lacking.  It  is  unclear  whether 
more  recently  the  frequency  of  BCS  has  surpassed  mastec¬ 
tomy  use  and  whether  rates  of  BCS  have  increased  in  re¬ 
gions  of  the  US  previously  demonstrating  low  rates.  In 
addition,  it  is  not  known  whether  socioeconomic  factors 
continue  to  be  strong  influential  factors  affecting  treat¬ 
ment.  Characterizing  contemporary  use  patterns  is  impor¬ 
tant  for  identifying  persistent  barriers  to  treatment  with 
BCS.  So,  in  a  contemporary  national  sample  of  Medicare  pa¬ 
tients  diagnosed  with  invasive  breast  cancer,  we  sought  to 
quantify  the  frequency  of  treatment  with  BCS  versus  mastec¬ 
tomy.  In  addition,  we  sought  to  determine  whether  socioeco¬ 
nomic  and  demographic  factors — such  as  geographic  region, 
income,  education,  and  the  supply  of  health-care  providers — 
influenced  surgical  treatment  in  patients  across  the  US. 

METHODS 

National  Medicare  dataset  and  study 
sample  derivation 

The  national  Medicare  dataset  includes  claims  data  for  all 
Medicare  beneficiaries  in  the  US.  Files  contain  data  col¬ 
lected  by  Medicare  for  reimbursement  of  health-care  ser¬ 
vices  for  each  beneficiary  and  include  institutional  (inpa¬ 
tient  and  outpatient)  and  noninstitutional  (physicians  or 
other  providers)  final  action  claims.7  Of  note,  the  national 
Medicare  dataset  is  distinct  from  the  Surveillance,  Epide¬ 
miology,  and  End  Results  (SEER)-Medicare  linked  data¬ 
base;10  our  Medicare  dataset  is  nationally  comprehensive 
and  SEER-Medicare  focuses  only  on  select  SEER  registry 
sites.  The  University  of  Texas  MD  Anderson  Cancer  Cen¬ 
ter’s  institutional  review  board  exempted  use  of  the  Medi¬ 
care  database  for  this  study. 

We  used  the  following  algorithm  to  identify  patients 
with  incident  breast  cancer,  treated  with  either  BCS  or 


mastectomy.  This  method  was  based  on  a  previous  vali¬ 
dated  algorithm  for  claims  data.11'12  We  included  women 
(age  65  years  or  more)  who  had  any  diagnosis  of  invasive 
breast  cancer  in  2003  (defined  as  an  ICD-9  diagnosis  code 
of  174)  who  underwent  BCS  (n  =  83,61 1)  or  mastectomy 
(n  =  42,504)  between  January  1,2003,  and  December  31, 
2003.  From  this  sample,  to  increase  the  specificity  of  the 
definition,  we  excluded  23,715  patients  who  did  not  have 
at  least  2  claims  (on  different  dates)  indicating  a  diagnosis 
of  invasive  breast  cancer  between  January  1,  2003,  and 
December  31,  2004  (at  least  1  claim  must  have  occurred 
during  2003),  with  no  more  than  6  months  between  the 
date  of  BCS  or  mastectomy  and  the  earliest  breast  cancer 
diagnosis  claim  date.  To  exclude  the  prevalent  cases,  we 
excluded  16,471  patients  who  had  a  breast  cancer-related 
diagnosis  or  procedure  claim  between  January  1,  2002,  and 
December  31,  2002.  To  reduce  misclassification  of  the  pri¬ 
mary  intended  operation,  we  excluded  630  patients  who 
underwent  both  types  of  operation  (either  date  of  mastec¬ 
tomy  claim  preceding  date  of  BCS  claim  or  mastectomy 
occurring  more  than  3  months  after  BCS),  except  for  pa¬ 
tients  who  received  a  mastectomy  within  3  months  of  BCS 
and  for  whom  mastectomy  was  considered  the  definitive 
operation.  To  limit  our  sample  to  patients  with  nonmeta¬ 
static  invasive  breast  cancer,  we  then  excluded  2,122  pa¬ 
tients  who  had  2  or  more  claims  specifying  metastatic 
breast  cancer  from  3  months  before  to  3  months  after  the 
diagnosis  date.  To  improve  sample  homogeneity,  we  also 
excluded  5,719  patients  who  were  receiving  Medicare  cov¬ 
erage  because  of  end-stage  renal  disease  or  disability.  Fi¬ 
nally,  to  ensure  we  had  complete  claims  information  to 
determine  patients’  cancer  treatment  course  and  comor¬ 
bidities,  we  excluded  6,612  patients  who  lacked  part  A  or  B 
coverage  or  who  had  intermittent  health  maintenance  or¬ 
ganization  coverage  in  the  9  months  after  or  in  the  1  year 
before  their  breast  cancer  diagnosis  date  (of  these  patients, 
3,561  had  incomplete  information  in  the  year  before  diag¬ 
nosis,  because  they  were  less  than  66  years  of  age).  For  this 
analysis,  the  breast  cancer  diagnosis  date  was  defined  as  the 
date  of  the  earliest  claim  for  a  diagnosis  of  breast  cancer. 
This  left  a  final  sample  size  of  56,725  patients. 

Surgical  treatment  and  other  covariates 

Covariates  derived  from  Medicare  files  (denominator  and 
claims  files)  included  cancer  treatment  variables,  other 
clinical  variables,  and  demographic  data.  Patients  were  clas¬ 
sified  as  treated  with  BCS  or  mastectomy  if  a  claim  for  the 
operation  (Appendix,  online)  occurred  within  6  months  of 
the  breast  cancer  diagnosis  date.  Claims  for  chemotherapy 
must  have  occurred  within  6  months  and  RT  claims  within 


Vol.  209,  No.  4,  October  2009 


Smith  et  al 


Breast-Conserving  Surgery  in  the  United  States  427 


9  months  of  the  breast  cancer  diagnosis  date.  These  claims- 
based  treatment  definitions  have  been  applied  in  previous 
studies  of  breast  cancer  patients.13'18 

Other  disease-  and  treatment-related  variables  included 
axillary  lymph  node  involvement,  axillary  lymph  node  dis¬ 
section,  sentinel  node  biopsy,  receipt  of  any  chemotherapy, 
specific  receipt  of  doxorubicin  or  paclitaxel,  receipt  of  any 
imaging  studies  for  staging,  number  of  hospitalizations  in 
the  year  after  diagnosis,  and  number  of  medical  oncology, 
radiation  oncology,  and  surgery  visits  in  the  year  after  di¬ 
agnosis.  Variables  indicating  preventive  health  care  and  in¬ 
teractions  with  the  health-care  system  included  mammog¬ 
raphy  in  the  year  before  diagnosis  and  number  of  physician 
visits  in  the  year  before  diagnosis.  In  addition,  we  calcu¬ 
lated  the  severity  of  comorbid  disease  for  each  patient  based 
on  a  modified  Charlson  comorbidity  score  validated  in  an 
earlier  claims-based  study:  0  (no  comorbidity),  1  (mild  to 
moderate),  or  2  or  more  (severe).19  This  score  combined 
comorbidities  recorded  in  Medicare  claims  during  the  12 
months  before  the  patient’s  cancer  diagnosis.  To  enhance 
specificity  of  comorbid  disease  diagnoses,  patients  must 
have  had  at  least  1  inpatient  (part  A)  claim  or  at  least  2 
outpatient  (part  B)  claims  more  than  30  days  apart.19 

Demographic  data  available  through  Medicare  files  in¬ 
cluded  patient  age  at  diagnosis,  race  (categorized  as  Cauca¬ 
sian,  African  American,  and  other),  and  state  and  county  of 
residence.  Classification  of  geographic  regions  was  based 
on  US  Census  Division  definitions.20  Socioeconomic  vari¬ 
ables,  obtained  from  the  2003  Area  Resource  File21  and 
linked  to  the  Medicare  dataset  by  county  of  patient’s  resi¬ 
dence,  included  median  household  income,  percent  living 
below  poverty  level,  and  percent  completing  ninth  grade 
education,  high  school,  or  college.  Supply  of  health-care 
providers  (for  breast  cancer  treatment)  was  measured  by  the 
density  of  general  surgeons  and  of  radiation  oncologists  at 
county  level,  obtained  from  the  2003  Area  Resource  File. 

Statistical  analysis 

All  analyses  were  conducted  using  SAS  software  version 
9.1.3  (SAS  Institute),  and  all  statistical  tests  assumed  a 
2-tailed  a  of  0.05.  We  calculated  percent  BCS  versus  mas¬ 
tectomy  use  for  the  entire  sample,  by  state  and  by  region. 
We  also  tested  the  unadjusted  bivariate  associations  be¬ 
tween  receipt  of  BCS  and  treatment,  clinical,  demo¬ 
graphic,  and  socioeconomic  covariates  using  the  Pearson 
chi-square  test  for  categorical  variables  and  the  Wilcoxon 
rank-sum  test  for  continuous  variables. 

A  multivariate  logistic  model  tested  adjusted  associa¬ 
tions  between  BCS  and  covariates.  Covariates  for  the  final 
multivariate  model  were  initially  selected  a  priori  based  on 
significance  in  bivariate  analyses  (p  <  0.25)  and  signifi¬ 


cance  in  earlier  studies  of  cancer  patients.19'24  The  final 
parsimonious  model  was  selected  based  on  statistical  sig¬ 
nificance,  goodness  of  fit,  and  minimizing  multicollinear- 
ity.  This  model  was  also  used  to  calculate  adjusted  rates  of 
BCS  use  by  state.  Goodness  of  fit  was  assessed  based  on  the 
model  R-square  value  and  the  Hosmer-Lemeshow  test. 

Early-stage  breast  cancer  subgroup 

Breast  cancer  stage  is  not  directly  available  through  Medi¬ 
care  claims  data.  To  select  a  subgroup  of  patients  with  early- 
stage  breast  cancer,  given  that  surgical  treatment  strategies 
are  typically  dependent  on  disease  stage,  we  applied  a  pre¬ 
viously  validated  predictive  algorithm  that  used  claims- 
based  covariates  to  identify  patients  with  a  high  probability 
of  having  stage  I  or  II  disease  (unpublished  data;  Appendix, 
online).  So,  in  the  selected  subgroup  of  43,706  predicted 
early-stage  patients,  we  further  examined  the  adjusted  as¬ 
sociations  between  receipt  of  BCS  and  covariates  using 
multivariate  logistic  modeling. 

In  a  secondary,  validating  analysis  on  this  selected  group 
of  predicted  early-stage  patients,  we  also  identified  the  sub¬ 
group  of 42,499  patients  who  did  not  have  claims  for  axil¬ 
lary  involvement  and  chemotherapy;  these  patients  would 
also  be  more  likely  to  have  early-stage  disease  (kappa  = 
0.73  for  the  2  selected  subgroups).  The  adjusted  associa¬ 
tions  were  compared  in  this  group  with  the  associations 
calculated  for  the  predicted  early-stage  patients. 

RESULTS 

Patient  characteristics  and  treatment  course 

Our  cohort  consisted  of  56,725  women  with  incident,  in¬ 
vasive  breast  cancer  diagnosed  in  2003  and  treated  with 
surgery.  In  our  sample,  median  age  was  76  years  (interquar¬ 
tile  range,  71  to  81  years).  Ninety  percent  (n  =  51,432) 
were  Caucasian,  7%  (n  =  3,727)  were  African  American, 
and  3%  (n  =  1,566)  were  of  other  race.  As  a  component  of 
the  initial  treatment  course,  the  majority  of  patients  were 
treated  with  BCS.  Specifically,  59%  of  patients  (n  = 
33,450)  received  BCS,  and  41%  (n  =  23,275)  underwent 
mastectomy.  Additionally,  of  the  entire  sample,  50%  re¬ 
ceived  RT  and  16%  received  chemotherapy.  Of  BCS  pa¬ 
tients,  74%  received  RT  and  13%  received  chemotherapy. 
Of  mastectomy  patients,  14%  received  RT  and  23% 
chemotherapy. 

Predictors  of  breast-conserving  surgery  versus 
mastectomy  use 

The  use  of  BCS  was  associated  with  both  clinical  and  non- 
clinical  factors.  In  the  entire  sample,  on  unadjusted  analy¬ 
sis,  patients  who  were  younger,  Caucasian,  had  fewer  co¬ 
morbidities,  lymph  node-negative  disease,  and  predicted 
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Table  1.  Patient  Characteristics  and  Use  of  Breast-Conserving  Surgery  Versus  Mastectomy  in  Women  (n  =  56,725) 

Patient  characteristic 

Treated  with  BCS 

Treated  with  mastectomy 

n 

% 

n 

% 

Demographic 

Age,  y,  mean  (SD) 

76  (7) 

77(7) 

66  to  <  70 

63 

7,112 

37 

4,260 

IV 

\] 

o 

58 

26,338 

42 

19,015 

Race 

Caucasian 

59 

30,572 

41 

20,860 

African  American 

54 

2,026 

46 

1,701 

Other 

54 

OO 

Wl 

to 

46 

714 

Clinical 

Charlson  comorbidity  score 

0  comorbid  conditions 

61 

22,735 

39 

14,578 

1  comorbid  condition 

57 

6,414 

43 

4,875 

2  or  more  comorbid  conditions 

54 

2,673 

46 

2,314 

Unknown 

52 

1,628 

48 

1,508 

Disease  stage  and  treatment 

Predicted  early-stage  (stage  I  or  II)  disease 

68 

29,828 

32 

3,622 

Axillary  lymph  node-positive  disease 

37 

17,587 

63 

30,050 

Lymph  node-negative  disease 

63 

5,688 

37 

3,400 

Axillary  lymph  node  dissection 

42 

11,995 

58 

16,715 

No  axillary  dissection 

77 

21,455 

23 

6,560 

Radiation  therapy 

88 

24,823 

12 

3,358 

No  radiation  therapy 

30 

8,627 

70 

19,917 

Chemotherapy 

44 

4,299 

56 

5,447 

No  chemotherapy 

62 

29,151 

38 

17,828 

Health  care  access* 

Median  surgeon  density*  (IQR) 

11  (8-16) 

10  (6-15) 

Median  radiation  oncologist  density*  (IQR) 

13  (4-20) 

11  (0-19) 

Socioeconomic  status* 

Living  in  metropolitan  area 

62 

25,979 

38 

16,132 

Living  in  nonmetropolitan  area 

51 

7,193 

49 

6,971 

Median  income,  $  (IQR) 

41,691  (36,221- 

-48,059) 

39,879  (34,267-45,922) 

Median  %  living  in  poverty  (IQR) 

10.7  (8-14) 

11.4  (9-14) 

Median  %  with  college  education  (IQR) 

24  (17-29) 

22  (15-28) 

p  values  are  <  0.001  for  all  data. 

*By  patient  county  of  residence. 

^Per  100,000  persons. 

*Per  10,000  persons. 

BCS,  breast-conserving  surgery;  IQR,  interquartile  range. 


early-stage  disease  were  more  likely  to  undergo  BCS.  Pa¬ 
tients  who  did  not  receive  chemotherapy  or  did  not  un¬ 
dergo  axillary  surgery  were  also  more  likely  to  undergo 
BCS  (Table  1).  In  addition,  neighborhood  socioeconomic 
factors  were  also  highly  associated  with  BCS.  Specifically, 
patients  living  in  metropolitan  areas  and  in  counties  with 
higher  median  household  income,  lower  percentage  living 
below  poverty  level,  and  higher  percentage  with  college 
education  were  more  likely  to  undergo  BCS.  Of  patients 
living  in  nonmetropolitan  areas,  only  51%  underwent 
BCS.  Supply  of  health-care  providers  also  influenced  treat¬ 


ment,  with  BCS  use  more  likely  in  patients  residing  in 
counties  with  a  higher  density  of  surgeons  and  radiation 
oncologists  (Table  1). 

Finally,  significant  geographic  variation  existed  (p  < 
0.001),  with  patients  in  the  Northeast  and  Pacific  West  most 
likely  to  undergo  BCS.  In  contrast,  patients  in  the  South  were 
least  likely  to  undergo  BCS,  with  half  or  fewer  of  all  patients  in 
these  regions  treated  with  BCS  (Fig.  1;  Table  2).  On  adjusted 
analysis,  higher  density  of  surgeons  was  no  longer  a  significant 
predictor  of  BCS  use  (p  =  0.13),  specifically,  once  the  multi¬ 
variate  model  accounted  for  geographic  region.  But  a  higher 
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Figure  1.  Percent  use  of  breast-conserving  surgery  (BCS)  by  state.  Darker  shading  represents  higher  fre¬ 
quency  of  use.  Not  shown:  Alaska  (55%)  and  Hawaii  (67%). 


density  of  radiation  oncologists  remained  a  significant  predic¬ 
tor  of  BCS  use  (p  =  0.01)  (Table  3). 

Early-stage  breast  cancer  subgroup 

In  the  selected  subgroup  of 43,706  patients  with  predicted 
early-stage  (stage  I  or  II)  disease  (77%  of  the  entire  sample), 
a  total  of  68%  (29,828  of  43,706)  of  this  selected  group 
received  BCS.  This  was  consistent  with  a  total  of  65% 
(27,544  of  42,499)  who  received  BCS  in  the  validation 
subgroup  of  patients  who  did  not  have  axillary  involvement 
and  did  not  receive  chemotherapy.  Geographic  variation 
persisted  in  the  use  of  BCS  for  patients  with  predicted 
early-stage  disease.  Patients  in  the  Northeast  (78%  to 
79%)  and  Pacific  West  (71%)  were  still  the  most  likely 


to  undergo  BCS;  patients  in  the  South  (57%  to  59%) 
and  portions  of  the  Midwest  (58%)  were  the  least  likely 
(Table  2).  The  validation  subgroup  was  similar,  with 
BCS  ranging  from  54%  (South)  to  75%  (Northeast). 
On  adjusted  analysis,  significant  predictors  of  BCS  use 
included  similar  demographic,  clinical,  and  socioeco¬ 
nomic  factors  as  predictors  for  the  entire  sample  (Table 
3).  Significant  predictors  of  BCS  use  identified  in  the 
validation  group  for  early-stage  breast  cancer  were  also 
consistent  with  this  analysis. 

DISCUSSION 

In  this  contemporary  national  cohort  of  older  patients,  we 
found  that  overall,  a  majority  of  surgically  treated  patients 


Table  2.  Unadjusted  and  Adjusted  Percent  Use  of  Breast-Conserving  Surgery  Versus  Mastectomy  by  Geographic  Region 


BCS,  %  (n  = 

23,275) 

Mastectomy,  %  (n 

=  33,450) 

Region 

States 

Overall 

Adjusted 

*  Early  stage 

Overall 

Adjusted* 

Early  stage 

West,  Pacific  West 

AK,  CA,  HI,  OR,  WA 

62 

62 

71 

38 

38 

29 

West,  Mountain  West 

AZ,  CO,  ID,  MT,  NV,  NM,  UT,  WY 

57 

58 

66 

43 

42 

34 

Midwest,  West  North  Central 

IA,  KS,  MN,  MO,  NE,  ND,  SD 

50 

51 

58 

50 

49 

41 

Midwest,  East  North  Central 

IL,  IN,  MI,  OH,  WI 

61 

60 

70 

39 

40 

30 

Northeast,  New  England 

CT,  MA,  NH,  ME,  RI,  VT 

70 

63 

79 

30 

36 

21 

Northeast,  Mid-Atlantic 

N],  NY,  PA 

67 

64 

78 

33 

35 

22 

South,  South  Atlantic 

DE,  DC,  FL,  GA,  MD,  NC,  SC,  VA, 
WV 

59 

59 

68 

41 

41 

32 

South,  West  South  Central 

AK,  LA,  OK,  TX 

50 

54 

59 

50 

45 

41 

South,  East  South  Central 

AL,  KY,  MS,  TN 

48 

51 

57 

52 

48 

43 

*Adj listed  for  covariates,  including  age,  race,  comorbidity  score,  axillary  lymph  node  involvement,  axillary  dissection,  chemotherapy,  screening  mammography, 
physician  visits,  surgeon  density,  radiation  oncologist  density,  metropolitan  area,  poverty,  and  education. 

BCS,  breast-conserving  surgery. 
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Table  3.  Multivariate  Logistic  Model: 

Predictors  of  Use  of  Breast-Conserving  Surgery 

Entire  sample 

Early-stage  subgroup 

Covariate 

Odds  ratio 

95%  Cl 

p  Value  Odds  ratio 

95%  Cl 

p  Value 

Age  66  to  <  70  y  vs  £  70  y 

1.37 

1.31 

1.44 

<  0.001 

1.15 

1.09 

1.22 

<  0.001 

Race 

Caucasian  versus  African  American 

1.14 

1.05 

1.23 

<  0.001 

0.95 

0.86 

1.04 

0.26 

Caucasian  versus  other  race 

1.29 

1.15 

1.45 

<  0.001 

1.17 

1.08 

1.27 

<  0.001 

Charlson  comorbidity  score 

0  versus  1  comorbid  condition 

1.18 

1.13 

1.24 

<  0.001 

1.21 

1.15 

1.28 

<  0.001 

0  versus  2  or  more  comorbid  conditions 

1.38 

1.29 

1.47 

<  0.001 

1.17 

1.08 

1.27 

<  0.001 

0  versus  unknown  comorbid  conditions 

1.12 

1.01 

1.25 

0.03 

0.91 

0.79 

1.05 

0.20 

Lymph  node-negative  disease 

1.60 

1.52 

1.68 

<  0.001 

— 

— 

— 

— 

No  axillary  lymph  node  dissection* 

4.00 

3.85 

4.17 

<  0.001 

4.65 

4.46 

4.88 

<  0.001 

No  chemotherapy* 

1.32 

1.25 

1.39 

<  0.001 

0.63 

0.57 

0.69 

<  0.001 

Screening  mammography 

2.02 

1.87 

2.17 

<  0.001 

0.89 

0.80 

1.00 

0.04 

One  or  more  visits  to  physician 

1.43 

1.22 

1.68 

<  0.001 

1.97 

1.60 

2.43 

<  0.001 

Geographic  region4 

West,  Pacific  West 

0.86 

0.78 

0.95 

0.008 

0.81 

0.72 

0.92 

0.001 

West,  Mountain  West 

0.71 

0.63 

0.80 

<  0.001 

0.69 

0.60 

0.80 

<  0.001 

Midwest,  West  North  Central 

0.50 

0.45 

0.55 

<  0.001 

0.45 

0.40 

0.51 

<  0.001 

Midwest,  East  North  Central 

0.80 

0.73 

0.88 

<  0.001 

0.75 

0.67 

0.84 

<  0.001 

Northeast,  Mid-Atlantic 

1.05 

0.95 

1.16 

0.34 

1.06 

0.94 

1.19 

0.36 

South,  South  Atlantic 

0.76 

0.69 

0.83 

<  0.001 

0.69 

0.62 

0.78 

<  0.001 

South,  West  South  Central 

0.58 

0.53 

0.65 

<  0.001 

0.54 

0.48 

0.61 

<  0.001 

South,  East  South  Central 

0.51 

0.45 

0.57 

<  0.001 

0.45 

0.39 

0.52 

<  0.001 

Surgeon  density* 

1.00 

0.99 

1.01 

0.13 

1.00 

0.99 

1.01 

0.14 

Radiation  oncologist  density^ 

1.30 

1.06 

1.59 

0.01 

1.37 

1.07 

1.75 

0.01 

Living  in  metropolitan  area 

1.20 

1.14 

1.26 

<  0.001 

1.20 

1.13 

1.27 

<  0.001 

Percent  living  in  poverty  <  11%^ 

1.05 

1.00 

1.09 

0.03 

1.05 

1.00 

1.11 

0.06 

Percent  with  college  education  >  23%^ 

1.13 

1.08 

1.19 

<  0.001 

1.15 

1.09 

1.21 

<  0.001 

*A  model  excluding  axillary  lymph  node  dissection  and  chemotherapy,  which  are  treatments  likely  to  occur  concurrently  or  after  surgery,  did  not  affect  risk 
estimates  for  other  covariates. 

tCompared  with  reference  category  Northeast  New  England.  The  likelihood  ratio  test  for  all  strata  of  the  variable  for  region  was  statistically  significant  (p  <  0.001). 
^Increased  odds  per  1  surgeon  per  100,000  persons  (continuous  variable). 

^Increased  odds  per  1  radiation  oncologist  per  10,000  persons  (continuous  variable). 

^Continuous  variables  dichotomized  at  the  median  value. 


with  invasive  breast  cancer  received  BCS  rather  than  mas¬ 
tectomy:  59%  of  the  entire  sample  and  up  to  68%  of  those 
with  predicted  early-stage  disease.  But  we  found  consider¬ 
able  variations  in  BCS  use  by  several  important  factors. 
First,  clear  regional  differences  existed.  In  some  regions  an 
exceptionally  high  percentage  of  patients  received  BCS  (up 
to  79%  in  predicted  early-stage  patients  in  the  Northeast), 
but  only  about  half  of  patients  in  areas  of  the  South  re¬ 
ceived  BCS.  Second,  neighborhood  socioeconomic  charac¬ 
teristics  also  appeared  to  influence  BCS  use,  with  the  pres¬ 
ence  of  county-level  poverty  and  lower  education  levels 
both  associated  with  lower  use  of  BCS.  Furthermore, 
health-care  access  also  appeared  to  influence  treatment  pat¬ 
terns,  with  patients  living  in  metropolitan  areas  and  areas 
with  a  higher  density  of  radiation  oncologists  more  likely  to 


receive  BCS.  Surgeon  density  did  not  appear  to  influence 
BCS  versus  mastectomy  use. 

Previous  studies  of  patients  treated  in  the  1980s  and 
1 990s  reported  lower  rates  of  BCS  use  compared  with  these 
more  contemporary  results.  For  example,  in  a  similar  co¬ 
hort  of  Medicare  patients  diagnosed  with  invasive  breast 
cancer  in  1986,  Nattinger  and  colleagues7  reported  only  a 
12%  use  of  BCS  just  after  the  initial  publication  of  the 
NSABP-B06  trial  in  1985. 1,4,7  Subsequently,  in  1990,  a 
National  Institutes  of  Health  Consensus  Development 
Conference  statement  declared  that  breast  conservation 
was  considered  the  preferred  treatment  choice  for  the  ma¬ 
jority  of  women  with  early-stage  disease.4  Since  that  time, 
steady  increases  have  occurred  in  BCS  use,8,22  but  studies 
have  consistently  found  that  the  majority  of  patients  with 
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invasive  breast  cancer  in  the  US  have  been  treated  with 
mastectomy.  For  example,  a  more  recent  study  by  Morrow 
and  associates6  found  that  in  1994,  BCS  use  was  still  only 
43%  in  patients  with  stage  I  or  stage  II  breast  cancer  across 
the  US.  Continued  evidence  has  accumulated,  providing 
increasingly  convincing  data  that  BCS  plus  radiotherapy  is 
comparable  to  mastectomy  in  early-stage  patients,  includ¬ 
ing  the  20-year  followup  of  the  NSABP-B06  trial  pub¬ 
lished  in  2002, 2,3  demonstrating  low  local  failure  rates 
after  breast  conservation.  With  this  published  data  it  is 
not  surprising  that  the  overall  frequency  of  BCS  in  our 
more  recent  cohort  demonstrates  continued  increases  in 
use.  Future  studies  may  seek  to  determine  how  trends  in 
mastectomy  versus  BCS  continue  to  evolve.  Recent  data 
suggest  a  reversal  in  use  trends  may  be  occurring,  with  an 
increasing  number  of  patients  receiving  mastectomy,  po¬ 
tentially  influenced  by  changing  technologies  in  breast 
cancer  care,  such  as  the  use  of  breast  MRI.23  Alterna¬ 
tively,  the  convenience  of  other  emerging  technologies, 
such  as  accelerated  partial  breast  irradiation  and  hypo- 
fractionated  whole  breast  irradiation,  may  increase  the 
frequency  of  breast  conservation. 

Substantial  geographic  variation  in  BCS  use  also  existed 
in  earlier  treatment  eras,  with  a  nearly  identical  pattern  to 
that  found  in  our  study.  Patients  in  the  northeastern  re¬ 
gions  of  the  US  have  consistently  been  the  most  likely  to 
receive  BCS  and  those  in  southern  regions  the  least  likely, 
and  this  pattern  has  persisted  over  approximately  the  last 
two  decades. s’7,24'25  Several  potential  causes  for  this  marked 
variation  have  been  studied.  First,  physician  influences  on 
surgical  treatment  choices  have  been  considered.  Various 
surgeon  characteristics  have  been  found  to  influence  the 
use  of  BCS  versus  mastectomy,  including  surgeon  gender, 
practice  setting,  geographic  setting  (rural  versus  urban), 
and  practice  volume.26'28  Other  surgeon  characteristics,  in¬ 
cluding  surgical  or  educational  background  and  surgical 
experience,  have  been  found  to  influence  other  aspects  of 
treatment,  such  as  receipt  of  radiotherapy  and  use  of  senti¬ 
nel  lymph  node  biopsy.29,30  In  addition,  patient  preferences 
and  the  patient-surgeon  interaction  influence  surgical 
treatment  choice.  Several  earlier  studies  have  noted  that 
even  when  offered  or  recommended  the  option  of  breast- 
conserving  therapy,  some  patients  still  choose  mastecto¬ 
my.31  Accordingly,  treatment  goals  should  account  for  in¬ 
dividual,  regional,  or  cultural  preferences  for  care. 

Finally,  previous  studies  of  patients  treated  in  earlier  eras 
demonstrated  that  socioeconomic  factors,  such  as  non- 
Caucasian  race,  lower  income,  higher  poverty  levels,  and 
lower  education,  were  significant  barriers  to  receipt  of 
BCS.32  Our  analysis  of  the  early-stage  disease  subgroup 
suggests  that  racial  differences  in  treatment  for  African- 


American  patients  and  Caucasian  patients  may  have  de¬ 
creased  compared  with  those  reported  in  earlier  studies.8 
Unfortunately,  the  persistent  finding  of  the  other  socioeco¬ 
nomic  barriers  in  our  contemporary  cohort  suggest  that 
little  progress  has  been  made  in  overcoming  disparities  in 
the  use  of  breast-conserving  therapy,  particularly  for  pa¬ 
tients  living  in  the  most  disadvantaged  neighborhoods  and 
regions.  Indeed,  initiation  and  completion  of  radiotherapy 
after  BCS  itself  has  been  found  in  earlier  studies  to  vary 
with  socioeconomic  factors.33,34  So,  strategies  to  improve 
access  to  radiotherapy35  may  also  serve  as  important  factors 
in  promoting  BCS,  especially  in  patients  who  would  prefer 
this  surgical  treatment  modality.  Notably,  a  unique  finding 
of  our  study  was  that  the  supply  of  radiation  oncologists 
(as  measured  by  the  density  of  radiation  oncologist  by 
county)  was  a  major  factor  affecting  the  choice  to  pursue 
BCS.  Other  than  geographic  region,  this  socioeconomic 
variable  had  the  largest  effect  size  for  association  with 
surgical  treatment  choice.  Previous  studies  have  noted 
that  increased  distance  from  radiation  oncologists  is  a 
barrier  to  treatment  with  BCS.35,36  Similarly,  our  result 
supports  the  hypothesis  that  access  to  radiation  oncolo¬ 
gists  continues  to  be  an  important  factor,  which  may 
affect  decision  making  by  both  patients  and  physicians 
regarding  surgical  treatment. 

Limitations 

Our  study  has  several  limitations  to  consider.  Our  sample 
focused  on  older  breast  cancer  patients  with  continuous 
fee-for-service  Medicare  Part  A  and  B  insurance  coverage. 
So,  future  studies  will  be  required  to  validate  the  magnitude 
and  significance  of  predictors  of  BCS  use  in  younger  pa¬ 
tients.  In  addition,  studies  are  also  needed  to  determine  the 
magnitude  of  potential  disparities  in  use  among  patients 
with  other  insurance  status,  given  that  earlier  evidence  sug¬ 
gests  that  insurance  status  may  influence  surgical  treatment 
strategies.32  Misclassification  of  the  outcome,  surgical 
treatment,  was  possible,  given  that  this  variable  was  claims 
based.  We  excluded  the  very  small  percentage  (1%)  of  pa¬ 
tients  who  had  both  claims  for  BCS  and  mastectomy,  but 
mastectomy  did  not  appear  to  be  the  definitive  operation 
(mastectomy  preceded  BCS  or  occurred  more  than  3 
months  after  BCS)  in  order  to  reduce  this  potential  source 
of  misclassification.  Additionally,  in  our  subsidiary  analy¬ 
sis,  claims-based  algorithms  were  used  as  a  proxy  for  early- 
stage  disease,  so  sensitivity  and  specificity  of  the  definitions 
were  not  100%.  Because  early-stage  disease  is  an  important 
determinant  of  whether  BCS  is  a  surgical  treatment  choice, 
it  is  likely  that  the  true  frequency  of  BCS  use  lies  between 
39%  and  68%.  So  future  studies  may  seek  to  validate  and 
refine  the  estimated  frequency  of  BCS  use  in  patients  with 
confirmed  early-stage  disease.  Finally,  although  our  multi- 
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variate  model  adjusted  for  a  large  number  of  covariates, 
residual  confounding  may  still  exist,  so  our  model  may  not 
be  able  to  distinguish  exactly  the  contribution  of  physician, 
patient,  geographic,  and  socioeconomic  factors  to  treat¬ 
ment  decisions.  Our  large,  national,  retrospective  cohort- 
based  dataset  serves  as  an  important  complement  to  other 
studies  that  have  more  detailed  data  on  the  interactions 
between  patients  and  physicians  and  patients  and  socioeco¬ 
nomic  factors.31 

In  conclusion,  in  this  contemporary  study  we  demon¬ 
strated  that  currently  more  than  half  of  older  women  across 
the  US  diagnosed  with  nonmetastatic  invasive  breast  can¬ 
cer  are  treated  with  BCS.  This  represents  substantially  in¬ 
creased  use  rates  compared  with  historical  data.  Yet  marked 
variation  across  the  US  persists  in  BCS  use.  Lack  of  BCS 
use  appears  to  be  associated  with  socioeconomic  disadvan¬ 
tage,  suggesting  that  persistent  barriers  to  breast  conserva¬ 
tion  exist. 
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Appendix.  Covariates  and  Associated  Claims  Codes  Derived  from  Medicare  Files 

Predictor  variable 

Time  period 
searched 

Source  file 

ICD-9  diagnosis 

ICD-9  procedure 

CPT 

Revenue 

center 

Demographic 

Age 

At  diagnosis 

Denominator 

Race 

At  diagnosis 

Denominator 

Extent  of  disease  at  diagnosis 

Axillary  LN  involvement 

3  mo  before  to  3 
mo  after 

1963 

Metastatic  disease 

3  mo  before  to  3 
mo  after 

1962,  1965-6,  197,  1970,  1971,  1972, 
1973,  1974,  1975,  1976,  1977, 

1978,  198,  1980,  1981,  1982,  1983, 
1984,  1985,  1986,  1987,  1988, 

19881,  19882,  19889 

Cancer  treatment 

No.  visits  to  surgeon 

In  the  year  after 
diagnosis 

Carrier  claims 

No.  visits  to  medical 
oncologist 

In  the  year  after 
diagnosis 

Carrier  claims 

No.  visits  to  radiation 
oncologist 

In  the  year  after 
diagnosis 

Carrier  claims 

Imaging  (CT,  MRI,  PET, 
or  bone  scan) 

In  the  3  mo  after 
diagnosis 

8801,  8703,  8741, 
8891,  8896, 

8874,  9218, 

9214 

70450,  70460,  70470, 
70551-3,71250, 
71260,  71270, 
72192-4,  74150, 
74160,  74170,  76700, 
78315,  78320, 
78812-6,  G0213-5 

Radiation  therapy 

In  the  9  mo  after 
diagnosis 

V580,  V661,  V671 

9221-7,  9229 

77401-77525, 

77761-77799 

0330,  0333 

Breast-conserving  surgery 

In  the  9  mo  after 
diagnosis 

8520-3,  8525 

19110,  19120,  19125, 
19160, 19162 

Mastectomy 

In  the  9  mo  after 
diagnosis 

8541-8 

19180,  19182,  19200, 
19220,  19240 

Axillary  LN  dissection 

3  mo  before  to 

1  y  after 

4023,  4051,  8543, 
8547 

38740,  38745,  19162, 
19200,  19220,  19240 

Chemotherapy  (any 
agent) 

In  the  9  mo  after 
diagnosis 

V581,  V662,  V672 

9925 

96400-96549, J9000-9, 
Q0083-5 

0331,  0332, 
0335 

Doxorubicin 

J9000-1, J9180 

Paclitaxel 

J9170, J9265 

(continued) 
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Predictor  variable 

Time  period 
searched 

Source  file 

ICD-9  diagnosis 

ICD-9  procedure 

CPT 

Revenue 

center 

Preventive  care  and 

interaction  with  health 

care  system 

No.  physician  visits 

In  the  y  before 
diagnosis 

All 

Screening  ammography 

In  the  y  before 
diagnosis 

V7611,  V7612 

8737,  8736 

77055-6,  77058-9, 
76090-2,  G0202, 
G0204,  G0206 

0401,  0403 

Influenza  vaccine 

In  the  y  before 
diagnosis 

V0481 

90658,  G0008 

General  health  status 

Any  hospital  admission 

In  the  y  after 
diagnosis 

Inpatient 

Charlson  comorbidity 
score 

In  the  y  before 
diagnosis  up 
to  1  mo  after 
diagnosis 

Per  prior  studies* 

*Klabunde  CN,  Potosky  AL,  Legler  JM,  Warren  JL.  Development  of  a  comorbidity  index  using  physician  claims  data.  J  Clin  Epidemiol  2000;53:1258-1267. 
CPT,  Common  Procedural  Terminology;  LN,  lymph  node;  PET,  positron  emission  tomography. 
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ABSTRACT 


Purpose 

Brachytherapy  is  a  method  for  delivering  partial-breast  irradiation  after  breast-conserving  surgery 
(BCS).  It  is  currently  used  in  the  community  setting,  although  its  efficacy  has  yet  to  be  validated 
in  prospective  comparative  trials.  Frequency  and  factors  influencing  use  have  not  been  previ¬ 
ously  identified. 

Methods 

In  a  nationwide  database  of  6,882  Medicare  beneficiaries  (age  a  65  years)  with  private 
supplemental  insurance  (MarketScan  Medicare  Supplemental),  claims  codes  identified  patients 
treated  with  brachytherapy  versus  external-beam  radiation  after  BCS  for  incident  breast  cancer 
(diagnosed  from  2001  to  2006).  Logistic  regression  modeled  predictors  of  brachytherapy  use. 

Results 

Frequency  of  brachytherapy  use  as  an  alternative  to  external-beam  radiation  after  BCS 
increased  over  time  (<  1  %  in  2001,  2%  in  2002,  3%  in  2003,  5%  in  2004,  8%  in  2005,  10% 
in  2006;  P  <  .001 ).  Increased  use  correlated  temporally  with  US  Food  and  Drug  Administration 
approval  and  Medicare  reimbursement  of  brachytherapy  technology.  Brachytherapy  use  was 
more  likely  in  women  with  lymph  node-negative  disease  (odds  ratio  [OR],  2.19;  95%  Cl,  1 .17 
to  4.1 1)  or  axillary  surgery  (OR,  1 .74;  95%  Cl,  1 .23  to  2.44).  Brachytherapy  use  was  also  more 
likely  in  women  with  non-health  maintenance  organization  insurance  (OR,  1.81;  95%  Cl,  1.24 
to  2.64)  and  in  areas  with  higher  median  income  (OR,  1.58;  95%  Cl,  1.05  to  2.38),  lower 
density  of  radiation  oncologists  (OR,  1 .78;  95%  Cl,  1.11  to  2.86),  or  higher  density  of  surgeons 
(OR,  1.57;  95%  Cl,  1.07  to  2.31). 

Conclusion 

Despite  ongoing  questions  regarding  efficacy,  breast  brachytherapy  was  rapidly  incorporated  into 
the  care  of  older,  insured  patients.  In  our  era  of  frequently  emerging  novel  technologies  yet 
growing  demands  to  optimize  costs  and  outcomes,  results  provide  insight  into  how  clinical,  policy, 
and  socioeconomic  factors  influence  new  technology  diffusion  into  conventional  care. 
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INTRODUCTION 


Brachytherapy  is  a  radiation  treatment  administered 
through  local  implantation  of  a  radioactive  source. 
In  recent  years,  brachytherapy  has  emerged  as  an 
important  technique  for  the  delivery  of  partial- 
breast  irradiation  after  breast-conserving  surgery 
(BCS),  particularly  in  selected  patients  with  breast 
cancer  with  low-risk  features  such  as  smaller  tumor 
size  and  node-negative  disease.  In  these  patients, 
either  multiple  interstitial  catheters  or  a  single  bal¬ 
loon  catheter  are  placed,  allowing  the  radiation 
treatment  to  be  delivered  locally  to  tissue  at  the  high¬ 
est  risk  for  recurrence.  Compared  with  a  standard 
course  of  external-beam  radiation  therapy  (EBRT) 
to  the  whole  breast,  a  brachytherapy  course  is  typi¬ 


cally  several  weeks  shorter,  thereby  reducing  incon¬ 
venience  associated  with  radiation  treatment. 

However,  the  use  of  breast  brachytherapy  as  the 
sole  radiation  treatment  after  BCS  remains  contro¬ 
versial.  EBRT  to  the  whole  breast  after  BCS  is  widely 
accepted  as  the  definitive  standard  therapy  for  early- 
stage  disease  based  on  extensive  randomized  studies. 
Data  have  accumulated  demonstrating  local  control 
and  survival  benefits  attributable  to  EBRT  over 
follow-up  as  long  as  20  years. 1-4  In  contrast,  phase  III 
trials  directly  comparing  breast  brachytherapy  with 
standard  EBRT  after  BCS  have  yet  to  mature.5 
Moreover,  existing  phase  II  studies  of  breast  brachy¬ 
therapy  have  generally  included  relatively  small 
sample  sizes  and  have  limited  median  follow-up 
times,  typically  around  3  to  5  years.6"12 
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Although  it  seems  that,  in  the  community  (nonclinical  trial) 
setting,  there  is  ongoing  use  of  breast  brachytherapy  after  BCS,  the 
actual  frequency  of  utilization  is  unknown.  Additionally,  the  factors 
influencing  its  use  have  not  been  previously  studied.  The  lack  of 
studies  on  breast  brachytherapy  utilization  patterns  is  surprising, 
given  that  substantial  controversy  exists  over  its  current  use  in  the 
community  setting.  Specifically,  some  experts  consider  this  treatment 
strategy  still  unproven  because  of  the  existing  gaps  in  the  scientific 
literature  regarding  its  efficacy.13  Furthermore,  others  considered  this 
strategy  one  of  the  costlier  options  for  patients  with  candidate  low-risk 
tumors  treated  in  the  era  of  its  initial  diffusion.14 

Documenting  the  pattern  by  which  this  novel — but  also  poten¬ 
tially  unproven  and  costly — treatment  has  diffused  into  the  care  of 
patients  with  breast  cancer  across  the  United  States  offers  a  unique 
opportunity  to  help  clinicians  and  policy  makers  better  understand 
how  clinical  factors,  policy  factors,  and  socioeconomic  factors  influ¬ 
ence  the  dissemination  of  new  technologies  into  the  health  care  sys¬ 
tem.  Therefore,  in  a  national  cohort  of  older,  insured  patients 
diagnosed  with  invasive  breast  cancer,  we  sought  to  quantify  the  fre¬ 
quency  of  breast  brachytherapy  utilization,  identify  time  trends  in 
breast  brachytherapy  use,  and  identify  clinical  and  nonclinical  factors 
associated  with  brachytherapy  use. 


METHODS 


Data  Set  and  Study  Sample 

The  MarketScan  Medicare  Supplemental  database  is  a  large,  nationwide, 
employment-based  claims  database  that  includes  Medicare  beneficiaries  with 
private  supplemental  insurance  (Appendix  Table  Al,  online  only).15  During 
the  study  interval,  approximately  two  thirds  of  Medicare  beneficiaries  had 
private  supplemental  insurance.  The  University  of  Texas  M.  D.  Anderson 
Cancer  Center’s  Institutional  Review  Board  exempted  this  study  for  approval 
because  all  patient  data  were  de-identified. 

We  developed  an  algorithm  to  identify  patients  with  incident  breast 
cancer  treated  with  BCS.  Our  method  was  based  on  a  prior  validated  algorithm 
for  claims  data.16,17  Our  initial  study  population  consisted  of  23,576  women 
who  had  any  diagnosis  claims  code  indicating  invasive  breast  cancer  (defined 
as  International  Classification  of  Diseases,  Ninth  Revision  diagnosis  code  of 
174)  plus  a  procedure  code  related  to  breast  cancer  (Appendix,  online  only) 
during  the  study  period  (2001  to  2006).  The  earliest  date  with  a  diagnosis  or 
procedure  code  indicative  of  breast  cancer  was  considered  the  index  date  (a 
proxy  measure  for  the  diagnosis  date).  After  applying  criteria  to  exclude 
prevalent  cases,  patients  younger  than  age  65  years  at  diagnosis,  men,  and 
patients  with  incomplete  treatment  course  information  (see  Appendix  for 
expanded  algorithm),  our  final  sample  included  6,882  women. 

Radiation  Therapy  and  Other  Covariates 

Patients  were  classified  as  treated  with  any  radiation  therapy  if  a  diagnosis 
or  procedure  claim  was  recorded  within  6  months  of  the  index  claim  date.  Of 
those,  patients  were  further  classified  as  treated  with  EBRT,  brachytherapy,  or 
both  (EBRT  plus  brachytherapy  boost),  as  indicated  by  claims  codes.  Patients 
treated  with  brachytherapy  were  further  classified  as  having  received  balloon- 
based  treatment  if  a  procedure  code  was  found  specifying  that  procedure 
(Appendix  1). 

Other  clinical  variables  were  defined  using  claims-based  criteria  and  were 
used  as  covariates.  A  complete  list  of  variables  is  provided  in  the  Appendix. 
Definitions  of  the  clinical  variables  were  derived  from  claims-based  definitions 
used  in  prior  studies  of  patients  with  breast  cancer.18  Published  clinical  guide¬ 
lines  on  the  use  of  breast  brachytherapy  offer  recommendations  based  on  age, 
tumor  size,  extent,  histology,  multicentricity,  receptor  status,  and  surgical 
history.19'21  Of  these  variables,  only  extent  (presence  of  lymph  node  involve¬ 
ment  or  metastatic  disease)  and  surgical  history  were  available  in  our  data  set  to 


be  evaluated.  Other  covariates  included  in  our  analysis  were  demographic 
variables  and  socioeconomic  variables,  listed  in  detail  in  the  Appendix.  For 
univariate  and  multivariate  analyses,  continuous  variables  were  categorized 
into  quartiles  a  priori. 

Statistical  Analysis 

All  analyses  were  conducted  using  SAS  software  version  9.1.3  (SAS  Insti¬ 
tute,  Cary,  NC),  and  all  statistical  tests  assumed  a  two-tailed  a  =  .05.  To 
address  our  first  objective,  we  calculated  percent  use  of  brachytherapy,  EBRT, 
and  EBRT  plus  brachytherapy  boost  for  the  entire  sample  and  by  year  of 
diagnosis.  To  address  our  second  objective,  we  evaluated  for  a  time  trend  using 
the  Mantel-Haenszel  y2  and  Cochran-Armitage  tests  for  trend.22  We  also 
benchmarked  the  time  trends  against  the  following  two  major  policy  events 
that  occurred  during  the  study  period:  US  Food  and  Drug  Administration 
(FDA)  approval  of  the  first  balloon-based  brachytherapy  device  for  breast 
cancer  in  June  2002  and  Medicare  reimbursement  of  breast  brachytherapy  in 
April  2004. 

Because  our  third  objective  was  to  understand  the  clinical  and  nonclini¬ 
cal  factors  influencing  selection  of  brachytherapy  alone  versus  conventional 
EBRT  and  because  the  number  of  patients  who  used  EBRT  plus  brachytherapy 
boost  was  small,  we  excluded  28  patients  who  received  EBRT  plus  brachyther¬ 
apy  boost  from  the  denominator  for  subsequent  analyses.  In  addition,  because 
few  patients  received  brachytherapy  alone  before  2003,  1,505  patients  treated 
before  2003  were  excluded  from  analyses  designed  to  address  our  third  objec¬ 
tive,  for  a  total  of  5,349  included  patients.  Multivariate  analyses  on  this  group 
included  5,031  patients  with  complete  covariate  information. 

To  address  our  third  objective,  we  tested  univariate  and  multivariate 
associations,  with  the  outcome  variable  dichotomized  as  receipt  of  brachy¬ 
therapy  alone  versus  EBRT.  The  univariate  associations  between  this  outcome 
and  clinical,  demographic,  and  socioeconomic  covariates  were  assessed  using 
the  Pearson  y2  test  for  categorical  variables  and  the  Wilcoxon  rank  sum  test  for 
continuous  variables.  Multivariate  logistic  models  tested  adjusted  associa¬ 
tions  (Appendix).22,23 


RESULTS 


Patient  Characteristics 

Of  6,882  women  treated  with  BCS  and  radiation  treatment,  the 
mean  age  was  75  years  (standard  deviation,  6  years).  Eight  percent  of 
women  had  axillary  involvement,  and  4%  had  metastatic  disease  at 
diagnosis.  Seventy-eight  percent  of  women  underwent  any  axillary 
surgery,  10%  received  chemotherapy,  and  65%  received  endocrine 
therapy  as  part  of  the  initial  treatment  course. 

Frequency  of  Brachytherapy  and  EBRT  Use  and 
Temporal  Trends 

Of  the  entire  sample,  333  (5%)  of 6,882  women  received  brachy¬ 
therapy  alone  (multicatheter  or  balloon-based),  6,521  (95%)  of  6,882 
women  received  EBRT,  and  28  (<  1%)  of  6,865  women  received 
EBRT  plus  brachytherapy  boost.  Treatments  with  brachytherapy 
alone  significantly  increased  over  time,  from  less  than  1%  in  2001, 
2%  in  2002,  3%  in  2003,  5%  in  2004,  8%  in  2005,  and  10%  in  2006 
( P  <  .001;  Table  1).  The  most  notable  increases  could  be  bench- 
marked  against  two  major  policy  events.  First,  an  increase  in  use  was 
noted  after  July  2002,  correlating  with  FDA  approval  of  the  balloon- 
based  breast  brachytherapy  device  (June  2002).  Second,  a  further 
increase  was  noted  after  July  2004,  correlating  with  Medicare  reim¬ 
bursement  of  treatment  (April  2004;  Fig  1).  Of  patients  treated  with 
any  form  of  brachytherapy  alone,  the  proportion  who  received 
balloon-based  treatment  also  increased  dramatically  over  time,  with 
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Table  1.  Temporal  Trends:  Patients  Receiving  EBRT  Alone,  EBRT  Plus  Brachytherapy  Boost,  or  Brachytherapy  Alone  During  the  Study  Period  (P<  .001) 

Time  Period 

Total  No. 
of  Patients 

Patients  Receiving 

EBRT  Alone 

Patients  Receiving 

EBRT  +  Boost 

Patients  Receiving 
Brachytherapy  Alone 

No. 

% 

No. 

% 

No. 

% 

1/1/2001  to  6/30/2001 

183 

181 

98.91 

0 

0.00 

2 

1.09 

7/1/2001  to  12/30/2001 

363 

359 

98.90 

2 

0.55 

2 

0.55 

1/1/2002  to  6/30/2002 

469 

466 

99.36 

1 

0.21 

2 

0.43 

7/1/2002  to  12/30/2002 

494 

478 

96.76 

1 

0.20 

15 

3.04 

1/1/2003  to  6/30/2003 

711 

691 

97.19 

1 

0.14 

19 

2.67 

7/1/2003  to  12/30/2003 

828 

796 

96.14 

4 

0.48 

28 

3.38 

1/1/2004  to  6/30/2004 

874 

827 

94.62 

3 

0.34 

44 

5.03 

7/1/2004  to  12/30/2004 

851 

804 

94.48 

2 

0.24 

45 

5.29 

1/1/2005  to  6/30/2005 

942 

865 

91.83 

8 

0.85 

69 

7.32 

7/1/2005  to  12/30/2005 

534 

489 

91.57 

2 

0.37 

43 

8.05 

1/1/2006  to  6/30/2006 

633 

565 

89.26 

4 

0.63 

64 

10.11 

Total 

6,882 

6,521 

94.75 

28 

0.41 

333 

4.83 

Abbreviation:  EBRT,  external-beam  radiation  therapy. 


89%  receiving  balloon-based  treatment  by  2006  (Fig  2) .  In  multivari¬ 
ate  analysis,  the  temporal  trend  indicating  a  steady  increase  in  the  use 
of  brachytherapy  remained  significant  (P  <  .001;  Table  2). 

Clinical  Predictors  of  Treatment  With 
Brachytherapy  Alone 

In  multivariate  analysis,  treatment  with  brachytherapy  alone  was 
also  associated  with  multiple  clinical  factors.  Brachytherapy  alone  was 
more  likely  in  women  who  had  lymph  node-negative  disease  (odds 
ratio  [OR],  2.19;  95%  Cl,  1.1 7  to  4. 11),  did  not  receive  chemotherapy 
(OR,  1.68;  95%  Cl,  1.01  to  2.80),  and  received  any  axillary  surgery 
(OR,  1.74;95%CI,  1.23  to  2.44;  Tables2and3). Reassuringly, only  2% 
of  patients  with  metastatic  disease  and  3%  of  patients  with  axillary 
nodal  involvement  received  brachytherapy  alone. 

Nonclinical  Predictors  of  Treatment  With 
Brachytherapy  Alone 

Several  nonclinical  factors  were  significant  predictors  of  treat¬ 
ment.  Geographic  region  was  a  strong  predictor  of  treatment  with 


Time  Period 


Fig  1.  Temporal  trends:  percentage  of  patients  treated  with  brachytherapy 
(multicatheter  or  balloon  based)  as  the  sole  modality  of  radiotherapy  after 
breast-conserving  therapy  (P  <  .001).  The  “a”  after  year  refers  to  January 
through  June,  and  “b”  refers  to  July  through  December.  FDA,  US  Food  and 
Drug  Administration. 


brachytherapy  alone.  Patients  residing  in  the  Northeast  were  less  likely 
to  receive  brachytherapy  alone,  particularly  compared  with  patients 
residing  in  the  West  and  South  (Tables  2  and  3).  Second,  health  care 
access  as  reflected  by  type  of  private  supplemental  insurance  was  also  a 
significant  predictor.  Brachytherapy  alone  was  more  likely  in  women 
whose  supplemental  insurance  was  non-HMO  (OR,  1.81;  95%  Cl, 
1.24  to  2.64). 

Finally,  several  neighborhood  socioeconomic  factors  marginally 
influenced  treatment.  Brachytherapy  alone  was  more  likely  for  pa¬ 
tients  living  in  areas  with  a  higher  median  income  (OR,  1.58;  95%  Cl, 
1 .05  to  2.38  for  the  second  quartile  compared  with  the  lowest  quartile) 
and  in  areas  with  a  lower  density  of  radiation  oncologists  (OR,  1.78; 
95%  Cl,  1. 1 1  to  2.86  for  the  second  quartile  compared  with  the  lowest 
quartile).  However,  these  comparisons  were  not  significant  for  higher 
quartiles.  In  addition,  in  the  provider-based  model  (socioeconomic 
variables  were  defined  based  on  provider  county),  brachytherapy 
alone  was  more  likely  for  patients  who  received  their  BCS  in  areas  with 
the  highest  quartile  density  of  surgeons  (OR,  1 .60;  95%  Cl,  1 .09  to  2.34 
for  the  highest  quartile  density  compared  with  the  lowest  quartile 


Fig  2.  Temporal  trends:  percentage  of  patients  treated  with  brachytherapy  as 
the  sole  modality  of  radiotherapy  after  breast-conserving  therapy  who  received 
balloon-based  brachytherapy  (P  <  .001).  The  “a”  after  year  refers  to  January 
through  June,  and  "b”  refers  to  July  through  December. 
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Table  2.  Multivariate  Logistic  Model 

Adjusted  OR  for  Predictors  of  Use  of  Brachytherapy  Alone  Versus  EBRT  Alone  in  the  Study  Period  From  2003  to  2006 

(N  =  5,031) 

Patient-Based  Model* 

Predictor  Variable 

OR 

95%  Cl 

P 

Age,  years 

<  70 

I.OOt 

70-79 

0.77 

0.57  to  1.03 

.08 

a  80 

1.08 

0.76  to  1.54 

.65 

Employee  relationship 

Employee 

1.00 

Spouse 

1.12 

0.85  to  1.46 

.42 

Geographic  region 

Northeast  (reference  group) 

1.00 

West 

2.69 

1 .46  to  4.93 

.0010 

Midwest 

1.64 

0.93  to  2.91 

.09 

South 

2.39 

1 .33  to  4.30 

.004 

Charlson  comorbidity  score 

0 

1.00 

>  1 

1.14 

0.68  to  1.90 

.62 

Any  inpatient  admission 

Yes 

1.00 

No 

0.87 

0.44  to  1 .69 

.68 

No.  of  outpatient  physician  visits 

1.01 

0.98  to  1.04 

.62 

Screening  mammography 

Yes 

1.00 

No 

0.70 

0.44  to  1.11 

.13 

Axillary  lymph  node  involvement 

Yes 

1.00 

No 

2.18 

1 . 1 7  to  4.10 

.02 

Metastatic  disease 

Yes 

1.00 

No 

1.30 

0.60  to  2.84 

.51 

Any  axillary  surgery 

No 

1.00 

Yes 

1.73 

1 .23  to  2.44 

.002 

Chemotherapy 

Yes 

1.00 

No 

1.67 

1.01  to  2.79 

.05 

Endocrine  therapy 

Yes 

1.00 

No 

0.93 

0.72  to  1.21 

.60 

Insurance  type 

HMO 

1.00 

Non-HMO 

1.80 

1.23  to  2.62 

.002 

Surgeon  density** 

Q1  (s  8) 

1.00 

Q2  (9-11) 

0.97 

0.65  to  1.46 

.90 

Q3  (12-15) 

0.91 

0.57  to  1.45 

.68 

Q4  (>  15) 

1.39 

0.83  to  2.33 

.21 

Radiation  oncologist  density 

Q4  (>  19) 

1.00 

Q1  (s  8) 

1.78 

1 .10  to  2.86 

.02 

Q2  (9-13) 

1.28 

0.83  to  1.96 

.26 

Q3  (14-19) 

0.91 

0.61  to  1.35 

.65 

Metropolitan  area 

No  (reference  group) 

1.00 

Yes 

1.30 

0.80  to  2.13 

.29 

Median  income, $ 

Q1  (s  40,305) 

1.00 

Q2  (40,306^15,290) 

1.51 

1.00  to  2.28 

.05 

Q3  (45,291-52,936) 

1.39 

0.84  to  2.31 

.20 

Q4  (>  52,936) 

1.69 

0.86  to  3.32 

.13 

(continued  on  following  page) 
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Table  2.  Multivariate  Logistic  Model:  Adjusted  OR  for  Predictors  of  Use  of  Brachytherapy  Alone  Versus  EBRT  Alone  in  the  Study  Period  From  2003  to  2006 

(N  =  5,031)  (continued) 

Predictor  Variable 

Patient-Based  Model* 

OR 

95%  Cl 

P 

Percentage  living  in  poverty 

Q1  (==  9.3) 

1.00 

Q2  (9.4-1 1 .8) 

0.96 

0.64  to  1 .44 

.83 

Q3  (11.9-15.2) 

1.45 

0.89  to  2.35 

.14 

Q4  (>  15.2) 

1.32 

0.75  to  2.32 

.34 

Percentage  with  college  education 

or  more 

Q1  (£  17.2) 

Q2  (17.3-24.9) 

0.86 

0.56  to  1 .30 

.47 

Q3  (25.0-29.5) 

1.46 

0.93  to  2.30 

.10 

Q4  (>  29.5) 

1.12 

0.64  to  1 .98 

.69 

Time  period  of  treatment 

Before  7/1/2004  (reference  group) 

1.00 

7/1/2004  to  12/31/2004 

1.58 

1.08  to  2.30 

.02 

1/1/2005  to  6/30/2005 

2.02 

1.43  to  2.86 

<  .001 

7/1/2005  to  12/31/2005 

2.01 

1.32  to  3.04 

.001 

1/1/2006  to  6/30/2006 

2.60 

1.81  to  3.76 

<  .001 

Abbreviations:  OR,  odds  ratio;  EBRT,  external-beam  radiation  therapy;  HMO,  health  maintenance  organization;  Q  quartile. 

*By  patient  county  of  residence. 

+An  OR  of  1 .00  indicates  the  reference  group. 
tPer  1000,000  persons. 

density),  although  this  difference  was  not  significant  for  other  quartile 
comparisons  (Table  2). 


DISCUSSION 


To  our  knowledge,  our  analysis  is  the  first  to  detail  frequency  and 
patterns  of  breast  brachytherapy  use  in  a  nonclinical  trial  setting  in  the 
United  States  and,  further,  to  identify  potential  clinical,  policy,  and 
socioeconomic  influences  on  treatment  utilization  in  older,  well- 
insured  women.  In  our  cohort  of  older  insured  women  diagnosed 
with  breast  cancer,  we  found  that  5%  of  patients  received  brachyther¬ 
apy  as  the  sole  radiotherapy  modality  after  BCS  during  the  study 
period.  More  importantly,  our  results  demonstrated  that  in  recent 
years,  the  use  of  brachytherapy  alone  after  BCS  increased  significantly. 
By  2006,  a  substantial  percentage  of  our  cohort — as  much  as  10% — 
was  treated  with  this  modality,  with  the  vast  majority  of  this  group 
using  a  balloon-based  delivery  system.  Temporal  trends  found  in 
our  study  could  be  benchmarked  both  against  FDA  approval  of  the 
balloon-based  brachytherapy  device  for  breast  cancer  and  Medi¬ 
care  reimbursement  of  breast  brachytherapy.  Interestingly,  FDA 
approval  was  directed  at  the  brachytherapy  balloon  device  only  and 
thus  did  not  evaluate  clinical  efficacy  data  (unlike  FDA  approval 
policy  for  pharmaceutical  agents).  However,  the  temporal  correla¬ 
tion  found  in  our  study  still  suggests  a  potential  influence  of  device 
approval  policy,  independent  of  accrual  of  clinical  evidence.  How¬ 
ever,  given  that  this  correlation  does  not  necessarily  imply  causa¬ 
tion,  the  effect  of  approval  policy  on  treatment  selection  may  merit 
further  investigation. 

To  our  knowledge,  no  prior  studies  of  breast  brachytherapy  have 
reported  the  frequency  of  its  use  in  the  community  setting.  The  ma¬ 
jority  of  published  studies  have  occurred  in  the  clinical  trial  setting, 


generally  limited  to  select  patients  who  meet  criteria  for  lower  risk 
tumors.6  In  contrast,  our  study  provides  new  data  to  suggest  that  the 
use  of  breast  brachytherapy,  and  particularly  balloon-based  tech¬ 
niques,  is  being  incorporated  both  rapidly  and  extensively  into  treat¬ 
ment  after  BCS.  In  a  recent  study,  an  estimated  13%  to  27%  of  patients 
treated  with  BCS  were  considered  potential  candidates  for  breast 
brachytherapy.24  Thus,  our  results  suggest  that  a  substantial  propor¬ 
tion  of  available  candidate  patients  are  currently  being  treated  with 
this  new  technology.  This  increasing  utilization  may  reflect,  in  part,  a 
growing  interest  by  both  patients  and  providers  in  the  convenience  of 
accelerated  partial-breast  irradiation  techniques.  For  breast  brachy¬ 
therapy,  the  typical  treatment  course  is  shortened  by  4  to  6  weeks 
compared  with  standard  EBRT  to  the  whole  breast,  thus  reducing  the 
inconvenience  of  a  long  course  of  therapy.  Alternatively,  increasing 
utilization  could  have  reflected  a  response  to  financial  incentives  for 
providers,  with  potentially  higher  reimbursements  for  this  technology 
during  the  study  period. 

The  relatively  high  frequency  of  use  and  the  temporal  trends 
indicating  a  rapid  diffusion  of  this  technology  reported  in  our  study 
are  particularly  noteworthy  in  the  context  of  unanswered  questions 
that  remain  regarding  the  most  appropriate  use  of  breast  brachy¬ 
therapy.  In  particular,  the  long-term  efficacy  of  breast  brachyther¬ 
apy  alone  after  BCS,  specifically  in  comparison  with  standard 
treatment  with  EBRT  to  the  whole  breast,  remains  unclear.  Only  a 
small  total  number  of  patients  have  been  accrued  in  randomized 
clinical  studies  evaluating  the  efficacy  of  all  accelerated  partial- 
breast  irradiation  techniques  including  brachytherapy,  yielding 
approximately  1 ,500  patient-years  of  follow-up.6'25  In  comparison, 
cumulative  data  from  randomized  trials  evaluating  the  added  ben¬ 
efit  of  EBRT  to  the  whole  breast  after  BCS  have  accumulated 
approximately  60,000  to  100,000  patient-years  of  follow-up  and 
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Table  3.  Univariate  Associations: 

Predictors  of  Use  of  Brachytherapy  Alone  Versus  EBRT  Alone  in  the  Study  Period  From  2003  to  2006 

Brachytherapy  Alone 

(n  =  294) 

EBRT  Alone  (n  = 

4,737) 

No.  of 

No.  of 

Variable 

Patients 

% 

Patients 

% 

P 

Demographic 

Age,  years 

Mean 

74 

75 

.25 

SD 

6 

6 

<  70 

79 

6.7 

1,098 

93.4 

.03 

70-79 

142 

5.1 

2,631 

94.9 

IV 

00 

o 

73 

6.9 

1,008 

93.1 

Employee  relationship 

.17 

Employee 

180 

5.5 

3,091 

94.5 

Spouse 

114 

6.0 

1,646 

93.5 

Geographic  region 

.03 

Northeast 

15 

3.3 

437 

96.7 

West 

103 

6.1 

1,595 

93.9 

Midwest 

97 

5.5 

1,676 

94.5 

South 

79 

7.1 

1,029 

92.9 

Clinical 

Charlson  comorbidity  score 

.43 

0 

275 

5.8 

4,479 

94.2 

>  1 

19 

6.9 

258 

93.1 

Any  inpatient  admission 

.48 

No 

283 

5.8 

4,593 

94.2 

Yes 

11 

7.1 

144 

92.9 

No.  of  outpatient  physician  visits 

.21 

Mean 

6 

6 

SD 

4 

4 

Screening  mammography 

.06 

No 

21 

4.0 

503 

96.0 

Yes 

273 

6.1 

4,234 

93.9 

Disease  stage  and  treatment 

Axillary  lymph  node  involvement 

.01 

No 

283 

6.1 

4,349 

93.9 

Yes 

11 

2.8 

388 

97.2 

Metastatic  disease 

.61 

No 

287 

5.9 

4,584 

94.1 

Yes 

7 

4.4 

153 

95.6 

Any  axillary  surgery 

.01 

No 

43 

4.1 

1,002 

96.0 

Yes 

251 

6.3 

3,735 

93.7 

Chemotherapy 

.03 

No 

276 

6.1 

4,262 

93.9 

Yes 

18 

3.7 

475 

96.4 

Endocrine  therapy 

.32 

No 

99 

5.4 

1,735 

94.6 

Yes 

195 

6.1 

3,002 

93.9 

Health  care  access* 

Insurance  type 

.005 

HMO 

59 

4.3 

1,309 

95.7 

Non-HMO 

235 

6.4 

3,428 

93.6 

Surgeon  density*! 

.24 

Q1  (s  8) 

79 

6.3 

1,182 

93.7 

Q2  (9-11) 

69 

5.3 

1,222 

94.7 

Q3  (12-15) 

60 

5.0 

1,136 

95.0 

Q4  (>  15) 

86 

6.7 

1,197 

93.3 

Radiation  oncologist  density*! 

.34 

Q1  (s  8) 

84 

6.6 

1,182 

93.4 

Q2  (9-13) 

71 

5.8 

1,146 

94.2 

Q3  (14-19) 

63 

5.0 

1,211 

95.1 

Q4  (>  19) 

76 

6.0 

1,198 

94.0 

(continued  on  following  page) 
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Table  3.  Univariate  Associations: 

Predictors  of  Use  of  Brachytherapy  Alone  Versus  EBRT  Alone 

in  the  Study  Period  From  2003  to  2006  (continued) 

Variable 

Brachytherapy  Alone 
(n  =  294) 

No.  of 

Patients 

% 

EBRT  Alone  (n 

No.  of 

Patients 

=  4,737) 

% 

P 

Socioeconomic  status* 

Metropolitan  area 

No 

28 

5.1 

520 

94.9 

.50 

Yes 

266 

5.9 

4,217 

94.1 

Median  income,  $ 

Q1  (<  40,305) 

54 

4.4 

1,188 

95.7 

.06 

Q2  (40,306-45,290) 

80 

6.4 

1,175 

93.6 

Q3  (45,291-52,936) 

72 

5.9 

1,160 

94.1 

Q4  (>  52,936) 

88 

6.7 

1,214 

93.2 

Percentage  living  in  poverty 

Q1  (<  9.3) 

83 

6.4 

1,216 

93.6 

.25 

Q2  (9.4-1 1 .8) 

61 

4.9 

1,185 

95.1 

Q3  (11.9-15.2) 

85 

6.5 

1,216 

93.5 

Q4  (>  15.2) 

65 

5.5 

1,120 

94.5 

Percentage  with  >  college  education 
Q1  (<  17.2) 

68 

5.4 

1,191 

94.6 

.08 

Q2  (17.3-24.9) 

56 

4.7 

1,128 

95.3 

Q3  (25.0-29.5) 

93 

7.1 

1,228 

93.0 

Q4  (>  29.5) 

77 

6.1 

1,190 

93.9 

Abbreviations:  EBRT,  external-beam  radiation  therapy;  SD,  standard  deviation;  HMO,  health  maintenance  organization;  Q,  quartile. 
*By  patient  county  of  residence. 
tPer  1000,000  persons. 


thus  conclusively  demonstrate  a  local  recurrence,  breast  cancer- 
specific  survival,  and  overall  survival  benefit.1,3 

Although  initial  studies  of  breast  brachytherapy  suggest  reason¬ 
able  short-term  locoregional  control  rates  (for  example,  a  3-year  ac¬ 
tuarial  recurrence  rate  of  2%  in  the  American  Society  of  Breast 
Surgeons  MammoSite  Breast  Brachytherapy  Trial7),  longer  term  re¬ 
sults  are  highly  anticipated,  given  the  long  natural  history  of  breast 
cancer  with  an  ongoing  risk  of  recurrence  even  after  several  decades. 
Current  published  guidelines  on  the  use  of  breast  brachytherapy  con¬ 
sider  the  limitations  of  available  clinical  data  and  are  generally  conser¬ 
vative  in  identifying  patients  at  the  lowest  risk  of  recurrence  as  the 
most  appropriate  candidates  for  this  therapy.19"21  Consistent  with 
guidelines,  nearly  all  patients  (94%)  in  our  study  treated  with  brachy¬ 
therapy  were  free  of  axillary  involvement  or  metastatic  disease  at  the 
time  of  treatment,  although  future  studies  with  explicit  tumor  staging 
information  may  seek  to  evaluate  in  more  detail  the  frequency  of 
guideline  concordance.  At  the  same  time,  less  consistent  with  available 
data,  our  results  demonstrated  the  dramatically  increasing  popularity 
of  the  balloon-based  technique  compared  with  the  interstitial  tech¬ 
nique  of  catheter  placement,  ostensibly  as  a  result  of  ease  of  balloon- 
based  catheter  placement.  The  majority  of  historical  data  on  breast 
brachytherapy  was  based  on  studies  using  interstitial  catheters,  and 
only  recently  have  larger  patient  cohorts  been  assembled  to  assess 
long-term  efficacy  of  balloon-based  techniques,7  which  may  help  to 
address  the  concern  over  potential  differences  in  treatment  volumes 
and  dosimetry  associated  with  the  balloon-based  system. 

Nonclinical  factors  also  played  an  important  role  in  influencing 
the  utilization  of  brachytherapy  alone,  with  our  results  suggesting  that 
the  use  of  brachytherapy  was  nonuniform  across  patient  subgroups 
and  geographic  regions.  Neighborhood  socioeconomic  factors  may 


have  influenced  patients’  access  or  physicians’  incentives  to  use  this 
novel  therapy.  Various  aspects  of  standard  breast  cancer  treatment, 
including  the  receipt  of  chemotherapy,  BCS,  and  radiotherapy  after 
BCS,  have  been  found  to  vary  by  sociodemographic  variables,  and  in 
general,  it  has  been  reported  that  socioeconomic  disadvantage  is  a 
barrier  to  appropriate  receipt  of  standard  therapy.18,26"28  Our  study  is 
hypothesis  generating,  suggesting  that  factors  associated  with  provider 
characteristics,  provider-patient  interaction,  insurance  environment, 
geographic  practice  patterns,  and  socioeconomic  status29"31  may  also 
affect  the  diffusion  of  emerging  treatments  in  breast  cancer  care.  The 
significance  of  low  density  of  radiation  oncologists  in  predicting  use  of 
brachytherapy  may  suggest  that  the  convenience  factor  associated 
with  brachytherapy  courses  may  indeed  have  affected  decisions  to 
implement  this  treatment  option.32  The  significance  of  high  density  of 
surgeons  in  predicting  brachytherapy  may  have  further  reflected  pa¬ 
tient  access  to  new  technology  or  may  have  potentially  reflected  pro¬ 
vider  behavior — for  example,  a  response  to  competitive  pressures  or 
financial  rewards  in  high-density  markets  or  availability  of  collabora¬ 
tive  knowledge  and  experience  using  new  treatments  associated  with 
higher  patient  volumes.  The  associations  with  socioeconomic  vari¬ 
ables,  which  require  further  validation,  could  suggest  a  unique  phe¬ 
nomenon  of  a  reverse  disparity,  in  which  more  advantaged,  better 
insured  patients  would  be  more  likely  to  receive  a  new,  although 
nonstandard,  treatment. 

Our  sample  focused  on  older  patients  with  breast  cancer  who 
were  well  insured  with  continuous  Medicare  plus  private  supplemen¬ 
tal  insurance  coverage.  Therefore,  future  studies  will  be  required  to 
validate  the  magnitude  and  significance  of  predictors  of  brachyther¬ 
apy  use  among  younger  patients  and  patients  with  less  comprehensive 
insurance  coverage,  in  whom  breast  brachytherapy  may  have  been 
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adopted  less  readily.3 1  Because  our  definition  of  radiotherapy  type  was 
claims  based,  misclassification  may  have  occurred.  However,  we  chose 
to  err  on  the  side  of  increased  specificity  for  the  claim  of  brachy- 
therapy  versus  EBRT,  with  patients  who  had  nonspecific  radiotherapy 
claims  codes  classified  as  receiving  EBRT.  Some  patient,  disease,  and 
treatment  details  were  unavailable  in  our  claims-based  database.  Fu¬ 
ture  studies  may  compare  breast  brachytherapy  utilization  with  other 
partial-breast  radiation  treatment  strategies  (eg,  three-dimensional 
conformal  therapy).  Additionally,  future  studies  may  focus  on  treat¬ 
ment  adherence  to  other  guideline  criteria,  such  as  younger  age,  tumor 
size,  histology,  multicentricity,  and  receptor  status.  Notably,  the  pa¬ 
tient  sample  size  (denominator)  decreased  in  our  analytic  sample  in 
the  last  two  time  periods.  This  potentially  reflects  the  influence  of  the 
implementation  of  Medicare  Part  D  on  January  1,  2006.  Residual 
confounding  attributable  to  this  factor  may  exist,  and  thus,  future 
studies  are  particularly  required  to  validate  whether  increased  brachy¬ 
therapy  use  during  this  time  period  existed  in  patients  with  various 
types  of  insurance  coverage.  Finally,  additional  studies  are  required 
to  identify  whether  there  were  any  underlying  changes  in  patient, 
provider,  or  institutional  factors  during  the  study  period  that  con¬ 
tributed  to  the  apparent  temporal  and  geographic  variations  in 
treatment  utilization. 

Despite  ongoing  questions  regarding  the  long-term  outcomes 
associated  with  breast  brachytherapy  and  controversy  over  its  routine 
use  outside  of  the  clinical  trial  setting,  breast  brachytherapy  has  been 
rapidly  incorporated  into  treatment  of  older,  well-insured  patients 
with  breast  cancer  in  recent  years,  with  the  frequency  of  use  as  high  as 
10%  in  older  insured  patients.  Our  results  suggest  that  the  availability 


of  clinical  evidence  demonstrating  treatment  efficacy  is  unlikely  to 
have  been  the  major  force  determining  diffusion  of  this  new  technol¬ 
ogy.  Instead,  nonclinical  factors,  including  policy  and  socioeconomic 
influences,  seem  to  play  an  important  role.  In  an  era  when  new  tech¬ 
nologies  and  therapies  are  advancing  rapidly  yet,  simultaneously, 
demands  are  growing  to  contain  costs  and  establish  treatment  effec¬ 
tiveness,  it  is  important  to  analyze  available  data  to  understand  how 
decisions  are  made  to  use  developing  treatments  such  as  breast 
brachytherapy.  Insights  gained  may  help  to  improve  the  rationale  by 
which  future  therapies  are  promoted  and  adopted  into  care. 
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2011  Gastrointestinal  Cancers  Symposium:  Register  Today! 


This  3-day  multidisciplinary  symposium  (January  20-22,  2011,  San  Francisco,  California)  is  the  premier  Gl  oncology  event, 
bringing  together  leading  experts  in  diagnostic  imaging,  gastroenterology,  medical  oncology,  radiation  oncology,  and 
surgical  oncology.  Tailored  to  be  a  discussion-based  meeting,  the  Symposium  is  intended  to  foster  dialogue  among 
oncologists  and  other  members  of  the  cancer  care  community.  Sessions  will  feature  invited  abstract  presentations  on 
the  latest  science  and  its  applicability  to  optimizing  treatment  of  patients  with  gastrointestinal  cancers. 

Symposium  highlights  include: 


•  Meet  the  Professor  Sessions— breakfast  and  evening  events  in  four  concurrent  sessions* 

•  Translational  Research  Sessions— daily  lunch  event 

•  Fellows,  Residents,  and  Junior  Faculty  Networking  Luncheon* 

•  NEW:  Posters  will  be  displayed  for  all-day  viewing 

•  NEW:  Q&A,  text,  tweet,  or  email  your  questions  in  selected  sessions 

•  NEW:  Abstract  supplement  in  Journal  of  Clinical  Oncology 
*Ticket  required 
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For  more  information,  please  visit  www.gicasym.org/2011. 
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